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of  Kearny  Point  in  Newark  Bay.  The  proposed  diversion  tunnel  and  surface  works 
are  generally  designed  to  protect  against  the  100-year  flood  event,  during  which 
it  is  anticipated  that  a  flood  volume  equal  to  that  of  Newark  Bay  will  be 
discharged  over  a  two -day  period. 

The  Corps  is  studying  the  impact  of  discharging  water  from  a  diversion 
tunnel  into  Newark  Bay  and  its  tributaries,  the  Passaic  and  Hackensack  Rivers. 
While  a  diversion  tunnel  discharge  itself  would  not  cause  significant  flooding, 
it  is  the  combination  of  a  storm  surge  generated  by  hurricanes  and  northeasters 
coincident  with  the  diversion  tunnel  operation  that  the  Corps  is  evaluating.  The 
main  purpose  of  this  study  is  to  determine  the  stage -frequency  curves  for  various 
locations  and  conditions  throughout  Newark  Bay  and  the  tidally  influenced  portion 
of  the  Passaic  and  Hackensack  Rivers.  These  curves  can  be  used  to  determine  the 
heights  of  storms  in  the  vicinity  of  Kearny  Point.  Another  objective  of  this 
study  is  to  examine  the  likelihood  that  a  significant  diversion  tunnel  discharge 
could  be  coincident  with  a  major  storm  surge,  and  to  simulate  such  a  combination, 
thereby  determining  the  potential  for  additional  flooding. 

During  the  1970's  and  1980's,  a  two-dimensional  hurricane  surge  model  was 
used  for  the  Newark  Bay/Passaic/Hackensack  River  system  as  boundary  forcing  for 
an  interior  Dynamic  Estuary  Model  (DEM)  to  estimate  flood  frequencies.  The 
current  study  builds  on  previous  models  and  revalidates  and  uses  a  modified  DEM 
to  accurately  simulate  tidal  elevations  in  the  system  using  recent  data  on  system 
bathymetry,  river  cross  sections,  topography,  and  historical  hurricane, 
northeaster,  and  rain  flow  data  through  1991,  in  combination  with  a  refined  grid. 
The  revalidated  DEM  is  a  valuable  engineering  tool  more  accurately  simulating 
tidal  and  river  flows  than  previous  studies. 

The  coincidence  analysis  confirmed  that  large  river  flows  are  unlikely  to 
coincide  with  large  surges  from  New  York  Bight.  Time  history  plots  of  several 
individual  events  show  that  there  is  a  weak  relationship  between  storm  tides  and 
river  flows,  and  that  the  river  flows  lag  the  tides  by  several  days;  large  storm 
tides  with  large  Passaic  River  flows  (lagged  by  0-2  days)  are  extremely  rare 
events.  An  increase  in  elevations  due  to  sea  level  rise  has  been  accounted  for 
in  the  stage -frequency  distribution  curves.  The  bootstrap  method,  a  modification 
of  the  joint  probability  method,  was  used  to  generate  stage -frequency 
distributions  with  33  synthetic  storms  approximating  the  entire  population  of 
historical  storms.  Results  of  coincidence  analysis  show  that  large  storm  tides 
together  with  large  river  flows  is  an  extremely  unlikely  combination,  and  they 
can  be  treated  as  essentially  independent  events.  However,  the  analysis 
indicates  that  the  likelihood  of  elevated  river  flows  following  storm  tides  (by 
0-2  days)  is  quite  strong. 

The  effect  of  tunnel  discharges  to  Newark  Bay  will  have  little  effect  (less 
than  0.1  foot)  on  tidal  elevations  in  the  lower  Passaic  River,  Hackensack  River, 
and  Newark  Bay  where  flood-prone  areas  are  located.  The  inclusion  of  future 
project  conditions,  such  as  control  levees  in  the  vicinity  of  Kearny  Point, 
changes  the  future  unimproved  stage -frequency  distributions  by  less  than  0.2  ft 
in  the  lower  parts  of  the  Passaic  and  Hackensack  Rivers  and  in  Newark  Bay. 
However,  the  effect  of  tunnel  diversion  into  Newark  Bay  will  be  to  reduce  river 
flooding  in  the  upper  tidal  Passaic  River  by  several  feet. 


PREFACE 


The  investigation  described  in  this  report  was  authorized  by  the  US  Army 
Corps  of  Engineers,  New  York  District,  Passaic  River  Division  (CENAN-PR-C) . 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

cubic  feet 

cubic  feet  per  second 
feet 

feet  per  second 

feet  per  second 
per  second 

inches 

miles  (U.S.  statute) 
nautical  miles 
knots  ( internat iona 1 ) 
square  feet 

square  feet  per  second 


_ Ex. _ 

0.028317 

0.028317 

0.3048 

0.3048 

0.3048 

25.4 

1.609347 

1.853187 

0.514444 

0.09290304 

0.09290304 


_ To  Obtain _ 

cubic  meters 

cubic  meters  per  second 

meters 

meters  per  second 

meters  per  second 
per  second 

millimeters 

kilometers 

kilometers 

meters  per  second 

square  meters 

square  meters  per  second 
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PASSAIC  RIVER  BASIH  FLOOD  PROTECTION  PROJECT 


STORM  SURGE  ANALYSIS 

PART  I :  INTRODUCTION 


Passaic  River  Flood  Protection  Project 

1.  The  935-square-mile  Passaic  River  Basin  in  New  Jersey  has  seen 
serious  flooding  during  the  last  century.  This  is  due  to  the  fact  that  it 
lies  in  the  East  Coast  storm  belt,  and  over  the  years  has  become  more 
urbanized,  with  about  30,000  structures  in  the  flood-prone  areas  (Fact  Sheet 
1992,  Appendix  D) .  Damages  average  about  $100  million  per  year,  and  the 
system  has  been  declared  a  national  disaster  area  on  a  number  of  occasions. 

2.  The  U.S.  Army  Corps  of  Engineers,  after  evaluating  more  than  150 
alternatives,  proposed  in  December  1987  to  reduce  the  upstream  flooding  threat 
by  constructing  a  dual  inlet  diversion  tunnel  system  as  the  key  to  a  Passaic 
River  Flood  Protection  Project.  The  Hater  Resources  Development  Act  of  1990 
authorized  funds  for  construction.  The  major  elements  of  the  project  are  two 
underground  tunnels,  a  20.1-mile-long  main  tunnel  about  40  feet  in  diameter 
and  a  1.2-mile  spur  tunnel  about  22  feet  in  diameter  (Figure  1)  that  convey 
Central  Basin  flood  waters  to  an  outlet  in  the  vicinity  of  Kearny  Point  in 
Newark  Bay.  The  overall  Passaic  River  Flood  Protection  Project  includes  the 
purchase  of  natural  water  storage  areas,  including  wetlands,  continued 
protection  for  flood-threatened  (but  protected)  areas,  and  the  preservation  of 
fish  and  wildlife.  The  proposed  diversion  tunnel  and  surface  works  are 
generally  designed  to  protect  against  the  100-year-flood  event.  During  such 
an  event,  it  is  anticipated  that  a  flood  volume  equal  to  that  of  Newark  Bay 
will  be  discharged  over  a  two-day  period.  The  Fact  Sheet  (Appendix  D) 
provides  an  overview  of  the  project. 

Study  Objectives 

3.  As  stated  in  the  Fact  Sheet,  the  U.S.  Army  Corps  is  studying  the 
impacts  of  discharging  water  from  a  diversion  tunnel  into  Newark  Bay  and  its 
tributaries,  the  Passaic  and  Hackensack  Rivers.  In  particular,  the  Hackensack 
River  has  been  the  scene  of  previous  flooding  due  to  storm  surges  entering  the 


9 


system  from  New  York  Bight.  As  recently  as  October  1991,  the  "Halloween" 
Northeaster  created  flooding  in  the  low-lying  areas  along  the  banks  of  the 
Hackensack,  including  the  Meadowlands.  At  this  time,  no  damages  were  reported 
for  this  event. 

4.  While  a  diversion  tunnel  discharge  by  itself  would  not  cause 
significant  flooding  in  these  areas,  it  is  the  combination  of  a  storm  surge 
(hurricane  or  northeaster  generated)  coincident  with  the  diversion  tunnel 
operation  that  the  Corps  is  evaluating.  The  main  purpose  of  this  study  is  to 
determine  the  stage- frequency  curves  for  various  locations  and  conditions 
throughout  Newark  Bay  and  the  tidally-inf luenced  portion  of  the  Passaic  and 
Hackensack  Rivers.  These  curves  can  be  used  to  determine  the  heights  of  storm 
tide  levees  in  the  vicinity  of  Kearny  Point.  Another  objective  of  this  study 
is  to  examine  the  likelihood  that  a  significant  diversion  tunnel  discharge 
could  be  coincident  with  a  major  storm  surge,  and  to  simulate  such 
combinations,  thereby  determining  the  potential  for  additional  flooding. 

The  Study  Area 

5.  The  study  area  for  the  storm  surge  analysis  includes  Newark  Bay 
and  the  tidal  portions  of  its  tributaries,  the  Passaic  and  Hackensack  Rivers. 
The  Passaic  River  is  tidal  as  far  upstream  as  Dundee  Dam,  a  distance  of  17.7 
river  miles.  Upstream  of  Dundee  Dam,  the  Passaic  River  Basin  has  a  drainage 
area  of  809.9  square  miles.  The  Hackensack  River  is  tidal  as  far  as  the 
Oradell  Reservoir,  a  distance  of  23  river  miles.  The  tributary  area  above  the 
New  Milford  gage  (at  Oradell  Dam)  is  113  square  miles.  Beyond  Newark  Bay, 
water  enters  the  system  through  Kill  Van  Kull  and  Arthur  Kill,  on  each  side  of 
Staten  Island.  In  recent  years  these  areas  have  been  dredged  to  deepen  and 
maintain  the  navigation  channel. 

Previous  Studies 


6.  During  the  late  1970s,  Camp  Dresser  &  McKee  Inc.  (CDM)  performed  a 
storm  surge  study  for  FEMA  (CDM  1983a).  The  FEMA  study  used  a  two- 
dimensional,  hurricane  surge  model  to  generate  boundary  conditions  at  four 
interior  points,  which  in  turn  were  used  as  boundary  forcing  for  an  interior 
model  of  the  Newark  Bay/Passaic  River/Hackensack  River  system.  The  interior 


10 


model  was  an  earlier  version  of  the  Dynamic  Estuary  Model  (DEM).  That 
version,  with  some  minor  modifications,  was  used  for  this  study.  The  CDM 
(1983a)  study  was  based  on  the  Joint  Probability  Method,  in  which  295 
synthetic  storms  (both  hurricanes  and  northeasters)  were  generated  assuming 
independence  of  individual  storm  parameters.  Assigning  probabilities  of 
occurrence  to  each  storm  parameter,  frequency  distributions  of  surge 
elevations  were  developed  for  key  locations  in  the  study  area. 

7.  The  study  was  extended  by  the  Corps  to  look  at  flood  frequencies 
on  the  Passaic  River  (CDM  1983b) .  The  DEM  grid  was  modified  to  provide 
improved  resolution  in  Newark  Bay,  the  Hackensack  River,  and  particularly  the 
Passaic  River. 

8.  In  1988,  another  study  examined  the  cost  and  benefit  of 
constructing  a  storm  surge  barrier  on  the  Hackensack  River  to  protect  flood- 
prone  areas  in  the  tidal  reach.  In  that  study  (CDM  1988),  the  DEM  grid  was 
refined  along  the  Hackensack  River,  and  a  subset  of  24  synthetic  storm  events 
was  used  to  simulate  surge  elevations  with  varying  proposed  tidal  barrier 
heights. 


Approach 

9.  This  current  study  builds  on  previous  models  described  above,  and 
now  simulates  the  interaction  of  storm  surges  and  coincident  diversion  tunnel 
discharges.  Using  recent  data  on  system  bathymetry  and  topography,  the 
Hackensack  River  version  of  the  DEM  grid  (CDM  1988)  was  modified  to: 

(1)  incorporate  the  most  recent  information  on  bathymetry  and  river  cross 
sections,  and  (2)  provide  increased  resolution  in  parts  of  Kill  Van  Kull, 
Arthur  Kill,  Newark  Bay,  and  the  lower  parts  of  the  Passaic  and  Hackensack 
Rivers.  This  information  was  used  to  create  three  grids:  (1)  a  pre-1980 
condition,  (2)  an  October  1991  condition,  and  (3)  a  future  (2050)  condition. 

10.  The  "pre-1980"  grid  was  used  to  revalidate  the  DEM  model  using 

the  historical  conditions  from  the  Passaic  study  (CDM  1983b).  This  consisted 
of:  (1)  an  astronomical  tide  of  11-13  October  1976;  (2)  a  hurricane  of  21 

September  1938;  (3)  Hurricane  Donna  of  12  September  1960;  (4)  a  northeaster  of 
24-26  November  1950;  and  (5)  a  northeaster  of  6-8  November  1953.  The  "October 
1991"  grid  was  then  used  to  validate  the  model  to  the  "Halloween"  northeaster 
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of  October  1991.  The  updated  1991  grid  included  partial  dredging  that  has 
occurred  in  Kill  Van  Kull  and  Arthur  Kill. 

11.  During  the  Hackensack  River  study  (CDM  1988),  the  coincidence  of 
storm  surges  and  river  flows  was  investigated  by  examining  flow  and  tide  gage 
records  around  the  time  of  observed  peak  surge  elevations.  In  this  study,  we 
reviewed  the  analysis,  and  extended  its  scope  to  determine  whether  a 
relationship  exists  between  storm  surges  and  proposed  diversion  tunnel  flows. 

12.  The  Hackensack  River  study  (CDM  1988)  used  a  modification  of  the 
joint  probability  method  to  identify  a  limited  number  of  storm  surge  events 
capable  of  providing  stage  frequency  distributions  that  would  approximate  the 
distribution  of  the  entire  population  of  synthetic  events.  Twenty-four  storms 
were  selected  for  the  CDM  (1988)  study.  The  modified  "bootstrap"  method  was 
further  examined  to  determine  whether  the  24  events  were  adequate  to  develop 
accurate  stage  frequency  distributions  for  the  purposes  of  this  study. 

13.  Several  project  alternatives  were  identified  for  detailed 
investigation  using  the  newly  revalidated  storm  surge  model.  Options  included 
various  tunnel  discharge  locations  and  the  inclusion  or  exclusion  of  other 
flood  protection  works  (i.e.,  storm  surge  levee  system).  These  various  design 
conditions,  and  the  diversion  tunnel  discharging  to  Newark  Bay  without  a 
coincident  storm  surge,  were  simulated  with  the  suite  of  storm  surges  and 
stage  frequency  distributions  developed.  Simulations  were  performed  both  with 
and  without  project  conditions  to  examine  the  project  impacts  on  storm  surges 
and  frequency  distributions. 

Findings  and  Recommendations 


14.  The  revalidated  DEM  model  more  accurately  simulated  tidal 
elevations  in  the  Passaic/Hackensack/Newark  Bay  system  than  previous  studies. 
This  was  due  to  inclusion  of  more  recent  system  geometry  data.  The  model, 
therefore,  is  a  valuable  tool  for  simulating  tidal  and  river  flow  scenarios  in 
the  system. 

15.  The  coincidence  analysis  confirmed  that  large  river  flows  are 
unlikely  to  coincide  with  large  storm  surges  from  New  York  Bight.  Time- 
history  plots  of  several  individual  events  show  that  there  is  a  weak 
relationship  between  storm  tides  and  river  flows,  and  that  the  river  flows  lag 
the  tides  by  several  days.  In  general,  for  a  relatively  small  surge  (say. 
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less  than  four  feet),  a  wide  range  of  river  flows  is  possible.  For  a 
relatively  small  (less  than  one-year)  river  flow,  a  wide  range  of  storm  surges 
is  possible.  However,  large  storm  tides  with  large  Passaic  River  flows 
(lagged  by  0-2  days)  are  extremely  rare  events.  The  data  were  further 
analyzed  to  determine  the  average  surge  condition  that  might  occur  within  72- 
hours  of  specific  significant  river  flows.  The  results  indicated  that  over 
the  range  of  river  flows  observed,  an  average  storm  surge  is  on  the  order  of 
1.5-2  feet.  In  a  similar  fashion,  the  data  were  also  analyzed  to  determine 
the  average  Passaic  River  flow  that  might  occur  within  72  hours  of  storm 
surges  of  different  elevations.  The  results  showed  that  for  all  storm  surge 
elevations  observed,  an  average  Passaic  River  flow  of  1000-2000  cfs  would  be 
found,  which  is  approximately  the  average  daily  river  flow  (1150  cfs). 

16.  The  National  Oceanic  Service  (NOS)  concluded  that  sea  level  is 
rising  at  an  approximate  rate  of  2.7  mm/yr  (0.0089  ft/yr).  Based  on  the 
previous  model  studies,  we  estimated  a  sea  level  rise  of  0.89  feet  in  the  100 
years  from  1950  to  2050.  Using  this  value  of  sea  level  rise,  a  sensitivity 
analysis  performed  for  future  (year  2050)  conditions,  for  three  tide  events 
and  four  river  flow  conditions,  demonstrated  that  the  distribution  and 
magnitude  of  river  and  tunnel  discharges  has  a  very  small  effect  (less  than 
0.1  ft)  on  tidal  elevations  in  the  lower  part  of  the  system,  including  the 
lower  tidal  parts  of  the  Passaic  and  Hackensack  Rivers  and  Newark  Bay.  The 
placement  of  flood  protection  levees  around  Kearny  Point  also  has  a  small 
effect  (less  than  0.2  ft  in  the  vicinity  of  Kearny  Point  and  0.4  ft  in  the 
vicinity  of  Newark  Airport)  on  tidal  elevations.  Upstream  in  the  Passaic 
River,  toward  Dundee  Dam,  the  effect  of  the  Diversion  Tunnel  is  to  lower  water 
elevations  during  major  river  flow  events.  For  the  100-year  flow,  for 
example,  even  though  the  peak  fresh  water  discharge  at  the  mouth  of  the 
Passaic  may  increase  (from  33,645  cfs  to  48,544  cfs),  flood  elevations  will  be 
lowered  by  several  feet  as  a  result  of  the  tunnel  bypassing  the  upper  part  of 
the  tidal  river.  The  flood  characteristics  of  the  Hackensack  should  remain 
essentially  the  same,  due  to  comparatively  larger  cross-sectional  areas. 

17.  Stage- frequency  distribution  curves  for  various  locations  in  the 
system  were  developed  using  33  synthetic  storm  events  and  the  one-year  river 
flows  in  the  Passaic  and  Hackensack.  The  new  distribution  curves  are  similar 
to  the  distributions  developed  previously  (CDM  1983a,  1988),  except  that  an 
increase  in  elevations  due  to  sea  level  rise  has  been  accounted  for  in  these 
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new  curves.  These  distributions  permit  the  selection  of  appropriate  flood 
protection  levees  around  Kearny  Point  to  protect  against  a  given  return  period 
storm  surge  event. 

18.  From  the  analyses  in  this  study,  the  following  conclusions  are 

drawn : 

1.  The  inclusion  of  future  project  conditions,  including  flood 
control  levees  in  the  vicinity  of  Kearny  Point,  changes  the  future 
unimproved  stage- frequency  distributions  by  less  than  0.2  foot  in 
the  lower  parts  of  the  Passaic  and  Hackensack  Rivers  and  in  Newark 
Bay. 

2.  The  stage-frequency  distributions  developed  for  the  vicinity  of 
Kearny  Point,  will  allow  the  selection  of  flood  protection  levee 
elevations  to  protect  against  storm  tides  with  a  given  return 
period.  The  distributions  in  the  upper  tidal  Passaic  and 
Hackensack  Rivers,  do  not  include  the  effect  of  significant  river 
flows,  and  therefore  reflect  only  tidal  conditions. 

3.  The  coincidence  of  large  storm  tides  and  large  river  flows  is  an 
extremely  unlikely  combination,  and  can  be  treated  as  essentially 
independent  events.  However,  the  likelihood  that  elevated  river 
flows  occur  following  storm  tides  (by  0-2  days)  is  quite  strong. 

4.  Sea  level  rise  is  estimated  to  be  on  the  order  of  2.7  mm/yr 
(0.0089  ft/yr),  or  0.89  feet  from  1950  to  2050. 

5.  The  magnitude  and  distribution  of  significant  river  and  tunnel 
discharges  has  a  very  small  effect  (less  than  0.1  ft)  on  tidal 
elevations  in  the  lower  Passaic  River,  Hackensack  River,  and 
Newark  Bay. 

6.  The  effect  of  tunnel  discharges  to  Newark  bay  will  have  little 
effect  (less  than  0.1  ft)  on  the  existing  characteristics  of 
flood-prone  areas  along  the  Hackensack  River,  such  as  the 
Keadowlands . 

7.  The  effect  of  tunnel  diversions  into  Newark  Bay  will  be  to  reduce 
river  flooding  in  the  upper  tidal  Passaic  River  by  several  feet. 

8.  The  tidal  boundary  condition  for  a  river  flow  model  of  the  Passaic 
River  system  should  be  the  superposition  of  an  average  tide  and  a 
storm  surge  of  1.5-2  feet. 
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9. 


The  coincidence  analysis  indicated  that  a  Passaic  River  flow  of 
1000-2000  cfs  may  occur  during  a  typical  storm  tide.  The  add 
"conservatism"  to  the  storm  tide  analysis  in  this  report,  we 
specified  the  one-year  river  flows  in  each  of  the  major 
tributaries.  This  results  in  a  flow  in  excess  of  5000  cfs  in  the 
Passaic  River,  or  nearly  three  times  the  average  flow  expected. 
However,  this  flow  will  occur  relatively  infrequently  with 
significant  storm  tides. 

10.  The  coincidence  analysis  was  used  to  develop  a  storm  tide  stage- 
frequency  distributions  of  events  in  which  the  Passaic  River  flow 
exceeded  6000  cfs.  This  analysis  can  be  used  to  determine 
appropriate  tidal  heights  for  design  of  tunnel  outlet  works. 

11.  Dredging  operations  in  the  Kills  and  Newark  Bay  will  have  little 
effect  on  circulation  and  water  levels  in  the  study  area.  Water 
levels  may  change  by  less  than  0.15  feet,  and  the  tidal  phase  may 
be  altered  by  less  than  one-half  hour. 
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Figure  1.  Pompton  River /Passaic  River  dual  inlet  tunnel 
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PART  II:  MODEL  MODIFICATIONS 


Description  of  Model 

19.  The  study  area  is  characterized  by  a  large  number  of  bays, 
embayments ,  interconnecting  channels,  waterways,  and  rivers.  The  offshore 
surge  model  developed  for  the  FEMA  study  (CDM  1983a)  used  a  spatial  resolution 
of  five  nautical  miles  for  the  computational  grid.  Due  to  the  grid  size  and 
the  lack  of  geometric  flexibility  of  its  finite-difference  scheme,  it  became 
necessary  to  select  a  different  model  for  the  nearshore  portion  of  the  study 
area.  The  Dynamic  Estuary  Model  (DEM)  was  selected  for  this  purpose.  DEM 
(WRE  1977)  is  a  two-dimensional,  vertically  integrated,  hydrodynamic  model 
based  on  a  link-node  computational  method.  The  fundamental  feature  of  the  DEM 
model  is  its  very  efficient  integration  of  the  one-dimensional  and  two- 
dimensional  domains  of  flow.  This  model  is  ideally  suited  to  modeling  the 
one-dimensional  Passaic  River  and  its  complex  connection  with  the  New  York 
Bight  through  Newark  Bay,  and  Arthur  Kill  and  Kill  Van  Kull  through  the  Lower 
New  York  Harbor.  The  model  was  modified  during  the  FEMA  study  (CDM  1983a)  to 
run  directly  using  a  binary  sequential  output  file  from  the  surge  generation 
model,  which  transfers  the  surge  heights  at  four  boundary  nodes  and  also  storm 
information. 

20.  The  hydrodynamic  behavior  of  a  water  body  influenced  by  external 
forces  is  governed  by  two  fundamental  equations;  namely,  the  continuity 
equation  and  the  momentum  equation.  A  complete  description  of  the  movement  of 
water  requires  a  simultaneous  solution  of  the  two  equations.  In  order  to 
solve  these  two  coupled  equations,  the  DEM  program  uses  a  link  and  node 
discretization  of  the  prototype  water  body  (Figure  2). 

21.  A  reach  (or  link)  is  defined  as  a  channel  where  water  flows  from 
one  end  to  the  other.  The  ends  of  the  reach  are  further  defined  as  nodes  (or 
junctions)  in  which  water  volumes  are  stored.  If  the  heads  at  the  end  nodes 
at  any  given  time  are  known,  the  hydraulic  gradient  along  the  channel  and  thus 
the  flow  rate  may  be  calculated.  By  the  transfer  of  water  from  one  end  to  the 
other  according  to  the  computed  flow,  the  head  at  the  node  will  subsequently 
be  altered,  defining  a  stepwise  hydraulic  computation  on  the  discretized 
system  representing  the  water  body. 
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22.  The  nodes  are  defined  by  the  surface  area,  depth,  and  a  starting 
water  surface  elevation.  The  nodal  surface  area  (which  can  vary  with  depth) 
is  defined  as  the  area  enclosed  by  the  perpendicular  bisectors  (Figure  2)  of 
its  connecting  links  and  the  coastline  for  nodes  near  the  shoreline.  To 
assure  mass  conservation,  the  sum  of  the  surface  areas  for  all  nodes  must 
equal  the  total  surface  area  of  the  water  body.  For  areas  in  which  the 
surface  area  varies  greatly  with  water  elevation,  the  program  has  an  option  to 
vary  the  surface  area  with  depth.  This  option  becomes  very  important  when 
modeling  the  storm  surge  elevations  in  flat  areas  such  as  the  Meadowlands. 

The  depth  of  the  node  is  defined  as  the  average  depth  of  the  area  that  the 
node  represents.  Depths  on  charts  are  referenced  to  some  datum,  usually  Mean 
Low  Water  (MLW) .  A  correction  factor  is  included  to  correct  this  to  other 
datum  planes.  Zn  this  study,  all  elevations  are  referenced  to  the  National 
Geodetic  Vertical  Datum  of  1929  (NGVD)  since  this  is  the  one  standard  datum 
that  does  not  change  with  time  or  space.  The  nodes  are  also  identified  by  an 
X  and  Y  coordinate.  Although  not  strictly  needed  by  the  computational  scheme, 
this  information  is  used  to  compute  link  lengths  if  desired,  and  to  compute 
the  northeaster  and  hurricane  wind  and  pressure  effects. 

23.  The  hydraulic  computation  requires  the  following  information  for 
links:  the  width,  length,  depth  or  hydraulic  radius  and  friction  coefficient 
(Manning's  "n").  Values  for  these  parameters  are  selected  to  represent  the 
conveyance  characteristics  of  the  water  body  between  the  nodes.  For  a  two- 
dimensional  water  body,  the  appropriate  width  is  often  difficult  to  assess. 

For  consistency,  the  width  is  defined  by  the  length  of  the  perpendicular 
bisector  of  the  channels  between  their  intersections  with  adjacent  bisectors. 
The  depth  of  the  link  represents  the  controlling  depth  for  the  link  or 
channel.  Often  this  is  the  average  of  the  depths  of  the  two  end  nodes  and  DEM 
will  compute  it  if  unspecified.  The  length  of  the  link  is  also  computed  from 
the  X  and  Y  coordinates  of  the  end  nodes  if  it  is  unspecified.  The  Manning's 
"n"  friction  coefficient  is  the  only  calibration  parameter  of  the  model 
assuming  a  correct  spatial  geometric-discretization  of  the  water  body, 
although  experience  has  found  that  a  correct  geometric  description  is  the  most 
important  factor  for  accurate  model  simulations.  The  range  of  hydraulically 
plausible  values  for  Manning's  "n"  is  usually  between  0.015  and  0.08,  although 
a  range  of  0.025  -  0.04  is  more  typical  of  open  water.  Values  greater  than 
0.04  may  be  found  in  very  shallow  or  overbank  vegetated  areas.  Finally, 
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identification  of  the  nodes  at  the  ends  of  the  link  are  also  needed  in  order 
to  define  the  connectivity  of  the  network. 

24.  The  one-dimensional  momentum  equation  is  written  for  each  link 

as: 


BV 

dt 


-  g—  -  gS,  *  gS. 


(1) 


where : 

V  =  Horizontal  velocity, 
t  =  Time, 

x  =  Distance, 

H  =  Water  surface  elevation  from  the  datum  plane, 

g  =  Gravitational  acceleration, 

S,  =  Energy  gradient,  and 
Sw  =  Wind  stress. 

The  energy  gradient,  Sf,  of  turbulent  flow  is  proportional  to  the  square  of 
the  mean  velocity  according  to  Manning's  equation: 


n 2  V  \V\ 
2.2  R*° 


(2) 


where : 

n  =  Friction  coefficient  (Manning's  "n"),  and 
R  =  Hydraulic  radius  of  the  channel  (or  link). 
The  wind  stress  term  is  given  by: 

S  =  U2  cos 

d  «Pw 


(3) 


where : 

K  =  Dimensionless  coefficient  varying  with  U 
=  (0.8  +  0.0475  U)  x  2.63  10* 

d  =  Depth  of  flow, 
p,  =  Air  density, 
pw  =  Water  density, 

U  =  Wind  velocity  in  mph,  and 
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♦  =  Angle  between  the  wind  direction  and  the  axis  of  the  channel. 

This  wind  stress  formulation  was  modified  from  the  original  program  (WRE  1977) 
to  provide  a  more  accurate  determination  of  the  water  surface  wind  stress 
under  gale  force  winds.  It  is  used  to  compute  the  surface  stress  from  the 
simulated  hurricane  winds.  For  northeasters,  a  planetary  boundary  layer  model 
was  used  to  compute  the  surface  stress  terms. 

25.  The  second  equation  is  the  continuity  equation,  which  states  that 
the  rate  of  water  flowing  into  a  node  through  links  or  specified  inflows 
affects  the  rate  of  change  in  the  water  surface  elevation  at  the  node: 


k 


(4) 


where: 

A„j  =  Surface  area  of  node  j , 

Qi  =  Flow  of  a  connecting  channel, 

Qj  =  Water  inflow  rate  to  node  j,  and 
k  *  Number  of  links  entering  node  j. 

26.  The  numerical  solution  of  Equations  1  and  4  entails  rewriting  both 
equations  in  finite  difference  form.  Two  types  of  integration  are  necessary 
for  their  solution;  namely,  a  space  integration  and  a  time  integration.  A 
stepwise  procedure  is  used  for  the  spatial  integration  of  Equations  1  and  4, 
whereby  the  momentum  equation  is  used  to  solve  for  link  flow  and  the 
continuity  equation  is  used  to  solve  for  the  nodal  head.  Both  equations  are 
explicitly  solved  for  the  unknown  variables.  The  hydraulic  computations 
within  a  given  time  interval  proceed  as  follows: 

1.  Compute  the  flow  rate  in  each  link  according  to  the  hydraulic 
gradient  and  other  hydraulic  conditions  existing  at  the  beginning 
of  a  time  interval. 

2.  Compute  the  rise  or  fall  of  the  water  surface  (head)  at  each  node 
based  on  the  link  flow  and  the  importation  or  withdrawal  of  water 
at  the  node. 

3.  Update  the  geometric  and  hydraulic  conditions  for  the  computation 
of  the  next  time  interval. 

27.  Time  integration  is  performed  numerically  by  using  a  modified  Euler 
technique.  Each  computational  time-step  is  subdivided  into  two  parts,  whereby 
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the  time  rate  of  change  of  a  variable  calculated  at  the  half  time-step  is  used 
to  project  the  variable  from  one  time-step  to  the  next  full  time-step.  The 
intermediate  time-step  computations  improve  the  stability  and  accuracy  of  the 
model.  Stability  of  the  scheme  for  each  link  is  governed  by  the  Courant- 
Friederichs-Levy  condition: 


A  t  < 


Ax 

fgh 


(5) 


where : 

At  =  Computation  time-step. 

Ax  =  Length  of  link, 

g  =  Gravitational  acceleration, 

h  =  Depth  of  water,  and 

Vgh  =  Celerity  of  the  free  gravity  wave. 

An  efficient  discretization  has  an  equal  maximum  allowable  (stable)  time-step 
for  all  links.  Therefore,  shallow  links  should  have  shorter  lengths  and  deep 
areas  should  have  longer  links.  The  stability  of  the  entire  model  is  governed 
by  the  time-step  of  the  most  critical  link.  The  above  stability  condition 
provides  a  guide  to  the  link  lengths  to  be  used  but  may  often  be  exceeded 
without  producing  stability  problems  in  the  DEM  model. 

28.  Two  types  of  boundary  conditions  are  used  in  the  DEM  model,  namely 
flow  (or  discharge)  boundary  conditions  and  tidal  (or  water  surface  elevation) 
boundary  conditions.  Flow  boundary  conditions  are  used  for  the  headwaters  of 
rivers,  whereas  tidal  boundary  conditions  are  used  at  the  mouth  of  estuaries. 
The  tidal  boundary  conditions  can  be  fitted  astronomic  tides,  measured  tidal 
elevations,  or  direct  input  from  the  offshore  surge  generation  model. 

29.  The  version  of  the  Dynamic  Estuary  Model  (DEM)  applied  for  the  FEMA 
purge  study  of  New  York  Bight  (CDM  1983a),  and  the  Corps  study  of  the  Passaic 
River  (CDM  1983b)  included  only  open  channels  (using  a  double  trapezoid)  as 
flow  conveyances.  In  the  Hackensack  River  study  (CDM  1988)  there  were  two 
types  of  physical  structures  that  needed  to  be  included  in  the  model,  culverts 
and  tidal  barriers. 

30.  The  Hackensack  River  basin  is  crossed  by  a  number  of  roads  and 
railways.  A  number  of  culverts  have  been  built  to  allow  the  uninterrupted 
f lows  of  many  of  the  creeks  in  the  area.  From  experience  with  culvert 
dynamics,  we  know  that  they  can  cause  dramatic  changes  in  local  hydrodynamics. 
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particularly  across  relatively  small  structures,  as  in  the  Bolsa  Chica  study 
(Walton  1988).  Therefore,  it  became  necessary  to  include  culvert  dynamics  in 
the  Hackensack  River  study  version  of  DEM  (CDM  1988). 

31.  The  main  impact  of  a  culvert  is  to  change  the  cross-sectional  area, 
with  friction  and  other  losses  being  relatively  minor.  There  are  a  number  of 
possible  water  surface  profiles  through  a  culvert  (Chow  1959;  Bodhaine  1968) 
but  for  relatively  short  and  flat  culverts,  the  following  generally  applies: 


H,,,  =  min 

<HC,  H„J 

(6) 

H*  =  min  (He, 

max  [H„,  h*,]) 

(7) 

where : 

H,  =  depth  at  upstream  end  of  culvert, 

H*  =  depth  at  downstream  end  of  culvert. 

He  =  height  of  culvert, 

=  depth  in  waterbody  upstream  of  culvert, 

H^,  *  depth  in  waterbody  downstream  of  culvert,  and 
=  critical  depth  at  outlet  of  culvert. 

Knowing  the  average  depth  of  flow  in  the  culvert  (0.5*  [H^+H,,,]),  the  cross- 
sectional  area,  wetted  perimeter,  and  the  hydraulic  radius  can  be  calculated 
for  either  a  box  or  circular  culvert. 

32.  It  was  necessary  to  include  the  capability  of  simulating  barrier 
overtopping  in  the  Hackensack  River  study  version  of  DEM  (CDM  1988).  This  was 
done  by  including  a  sharp-crested  weir  equation: 

Q«  =  3.3  B  H*2  -  Qp  (8) 

where: 

Q.  =  net  flow  over  the  barrier, 

B  =  width  of  the  barrier, 

B  =  height  of  surge  above  the  barrier,  and 

Q,  *  the  river  flow  pumped  around  the  barrier. 

The  last  term,  Q,,,  was  included  to  simulate  che  capability  of  the  tidal 
structure  to  pass  the  river  flow  around  the  barrier  using  a  pump  station.  The 
volume  of  water  flowing  over  the  barrier  during  each  time  step  is  then  Q„ 
multiplied  by  the  time  step,  with  the  flow  assumed  to  remain  constant  during 
each  time  step. 
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Data  Sources 


33.  The  original  DEM  link-node  model  (CDM  1983a)  was  developed  using 
NOAA  charts  and  USGS  quads  of  the  New  York  Bight  area.  For  the  original 
extension  to  the  Passaic  River  (CDM  1983b),  these  data  were  supplemented  with 
surveyed  cross-sections  of  the  Passaic  River  for  a  fluvial  flood  investigation 
using  the  computer  model  HEC2.  For  the  Hackensack  River  Surge  Barrier  study 
(CDM  1988),  refinements  were  made  to  the  overbank  areas  of  the  Hackensack 
River  flood  plain  using  overbank  topographic  maps  and  "as-built"  drawings  of 
the  major  bridges  in  the  area. 

34.  Since  the  Hackensack  River  study,  additional  information  has  been 
collected  in  the  study  area,  and  was  provided  by  the  Passaic  River  Division 
(PRD)  of  the  New  York  Corps  of  Engineers  (NAN).  Additional  information 
consisted  of: 

1.  Kill  Van  Kull  and  Newark  Bay  Channels,  Navigation  Improvement 
Project.  U.S.  Army  Engineer  District,  Corps  of  Engineers,  New 
York,  NY  (no  date).  Additional  dredging  information  was  obtained 
from  the  U.S.  Corps  of  Engineers  Waterways  Experiment  Station 
(WES),  Vicksburg,  MS. 

2.  Survey  of  Hackensack  River,  Condition  Survey  -  Echo  Soundings. 

File  No.  185,  Surveyed  by  Tallamy,  Van  Kuren,  Gertis  &  Associates, 
Orchard  Park,  New  York,  for  New  York  Department  of  the  Army,  New 
York  District  Corps  of  Engineers,  New  York,  NY  (April  1989). 

3.  Survey  of  Newark  Bay,  NJ,  Hackensack  River,  Echo  Soundings 
Condition  Survey.  File  No.  186,  Surveyed  by  Ocean  Surveys  Inc., 
Old  Saybrook,  Connecticut,  for  Department  of  the  Army,  New  York 
District  Corps  of  Engineers,  New  York,  NY  (1990). 

4.  Survey  of  Passaic  River,  Condition  Survey  -  Echo  Soundings.  File 
No.  797,  Surveyed  by  Tallamy,  Van  Kuren,  Gertis  &  Associates, 
Orchard  Park,  New  York,  for  New  York  Department  of  the  Army,  New 
York  District  Corps  of  Engineers,  New  York,  NY  (June  1989). 

5.  Survey  of  Passaic  River,  Condition  Survey  -  Echo  Soundings.  File 
No.  798,  Surveyed  by  Tallamy,  Van  Kuren,  Gertis  &  Associates, 
Orchard  Park,  New  York,  for  New  York  Department  of  the  Army,  New 
York  District  Corps  of  Engineers,  New  York,  NY  (November  1989). 
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6.  Hackensack  River  (Overbank)  from  Newark  Bay  to  Oradell  Reservoir. 
Prepared  by  Aerial  Data  Reduction  Associates,  Inc.,  Pennsauken, 

NJ,  for  the  Department  of  the  Army,  New  York  District  Corps  of 
Engineers,  New  York,  NY.  Survey  photography  (March  1985). 

7.  Passaic  River  Flood  Control  Project  (DACW54-90-C-003)  with  UNET 
Cross  Section  Locations.  Prepared  by  Aerial  Data  Reduction 
Associates,  Inc.,  Pennsauken,  NJ,  for  the  Department  of  the  Army, 
New  York  District  Corps  of  Engineers,  New  York,  NY.  Survey 
photography  (January  1990). 

Pre-1980  Grid  Refinements 

35.  The  additional  data  sources  were  used  to  refine  the  model  grid  used 
for  previous  simulations  of  water  level  elevations  in  the  Kills,  Newark  Bay, 
Passaic  River,  and  Hackensack  River.  Refinements  included  adding  nodes  and 
links  in  some  areas  and  modifying  some  node  and  link  parameters.  Figure  3 
shows  the  modified  link-node  grid  used  for  this  study. 

36.  Additional  nodes  were  added  in  Kill  Van  Kull  and  Arthur  Kill  in 
order  to  fit  the  geographical  shape  more  closely  and  increase  resolution. 

Five  nodes  were  inserted  between  the  eastern  end  of  Kill  Van  Kull  and  the 
confluence  of  the  Rahway  River  with  Arthur  Kill  (Figure  4).  The  link-node 
description  of  Newark  Bay  (Figure  4)  was  significantly  refined  with  the  new 
model  using  about  twice  as  many  nodes  in  that  area  as  in  the  previous  studies. 
Nodes  were  added  to  distinguish  shallow  areas  from  the  main  channel,  and 
provided  pathways  for  lateral  water  exchanges  to  occur.  Nodes  were  also  added 
near  Kearny  Point  to  provide  additional  resolution  for  improved  condition 
simulations  that  included  the  proposed  diversion  tunnel  outflows.  The 
increased  resolution  allowed  for  a  finer  level  of  detail  in  the  examination  of 
effects  on  Newark  Bay  water  levels  of  possible  future  modifications  to  the 
existing  waterways.  One  additional  node  was  added  in  both  the  Passaic  and 
Hackensack  Rivers  slightly  north  of  Kearny  Point.  These  nodes  were  placed  in 
the  existing  sharp  bends,  which  both  rivers  make,  before  resuming  general 
northward  directions.  Finally,  two  nodes  were  added  to  include  the  Saddle 
River  in  the  model  system.  The  Saddle  River  flows  into  the  Passaic  River 
approximately  one-half  mile  south  of  Garfield. 


24 


37.  Soundings  and  contour  maps  of  the  Kills,  Newark  Bay,  and  Passaic 
and  Hackensack  Rivers  provided  the  source  for  refining  the  physical 
parameterizations  of  nodes  and  links.  Physical  parameters  for  the  newly-added 
nodes  and  links  were  taken  from  the  maps. 

1991  Grid  Refinements 


38.  Refinements  to  the  grid  to  represent  October  1991  conditions,  for 
simulation  of  the  "Halloween  Storm" ,  consisted  of  revising  node  and  link 
physical  parameters.  These  changes  were  made  to  the  pre-1980  grid.  Revisions 
reflected  depth  and  area  changes  due  to  the  ongoing  dredging  project  in  the 
Kills  and  Newark  Bay.  Generally  the  main  channels,  except  for  Arthur  Kill, 
were  deepened,  and  in  some  cases  slightly  widened.  However,  the  depth  of 
Arthur  Kill  remained  at  35  feet.  The  authorized  project  depth  for  other  areas 
of  the  system  was  40  feet,  although  actual  depths  were  up  to  one  foot  deeper. 
Extraction  of  the  rock  substrate  was  not  included  at  this  time,  since  the 
dredging  removed  only  the  covering  sediment  layer  using  a  suction  dredge.  The 
major  change  in  the  Newark  Bay  area  was  the  removal  of  a  large  shoal  off  the 
Port  Elizabeth  Marine  Terminal.  The  appropriate  nodes  and  links  were 
modified.  Other  physical  changes  resulted  in  model  modifications  at  corners 
of  connecting  waterways.  Only  slight  changes  to  channel  widths  resulted. 

39.  There  were  no  new  nodes  added  to  the  1991  grid.  The  additional 
resolution  and  enhanced  physical  representation  were  purposefully  added  to  the 
pre-1980  grid  in  anticipation  of  the  needs  of  subsequent  simulations, 
including  model  validation. 


Future  Project  Grid  Refinements 


40.  Prior  to  the  October  1991  storm,  mainly  sediment  overburden  had 
been  removed  from  the  waterways,  and  removal  of  the  rock  substrate  had  not 
begun.  Underlying  rock  substrate  removal  began  after  October  1991.  According 
to  dredging  plans,  rock  substrate  removal  was  to  be  performed  in  very 
localized  areas  rather  than  uniformly  over  the  waterways.  Removal  thicknesses 
were  generally  small.  Because  model  predictions  are  insensitive  to  small 
changes  in  waterway  depth  and  volume,  and  because  revisions  of  the  1991  grid 
to  account  for  rock  substrate  removal  would  not  necessarily  be  accurate  enough 
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to  reflect  the  actual  changes,  watering  geometry  and  parameterization  for  the 
future  project  grid  were  taken  to  be  the  same  as  the  1991  grid.  At  present, 
no  other  plans  are  known  for  future  projects  to  modify  the  waterways. 

However,  some  simulations  of  future  conditions  included  modification  of 
overbank  links  and  nodes,  thereby  reflecting  the  proposed  flood-protection 
levees  at  Kearny  Point.  Waterway  parameters  were  not  affected. 
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Figure  2.  Pseudo-two-dimensional  geometric  representation 
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Long  Island  Sound 


PART  III:  MODEL  VALIDATION 


Introduction 

41.  Once  the  model  refinements  were  completed  (see  Figures  3  and  4), 
the  model  was  revalidated  against  six  events: 

1.  The  astronomic  tide  of  11-13  October  1976 

2.  The  hurricane  of  21  September  1938 

3.  The  hurricane  (Donna)  of  12  September  1960 

4.  The  northeaster  of  24-26  November  1950 

5.  The  northeaster  of  6-7  November  1953 

6.  The  "Halloween  Storm"  of  29-31  October  1991 

42.  In  a  previous  study  of  the  area  (CDM  1983b)  the  two  hurricanes,  and 
the  1950  and  1953  northeasters,  were  simulated  first  using  the  offshore  model, 
and  then  transferred  to  the  interior  link-node  model  at  exterior  boundary 
nodes  1-4.  Nodes  1-3  are  located  in  New  York  Bight,  approximately  two  miles 
east  of  Sandy  Hook.  Node  4  is  located  in  Long  Island  Sound  in  the  vicinity  of 
Bridgeport,  Connecticut.  However,  in  this  study  we  chose  to  force  the  link- 
node  model  with  direct  observations  at  Sandy  Hook  (node  12)  and  The  Battery 
(node  62)  in  Figure  3.  Therefore,  we  developed  a  smaller  version  of  the  link- 
node  grid  containing  176  nodes  and  263  links.  The  two  hurricanes  and  three 
northeaster  events  were  simulated  using  this  grid,  and  the  astronomical  tide 
was  simulated  using  the  entire  grid,  with  forcing  at  the  exterior  boundary 
nodes  1-4  and  at  the  interior  nodes  at  Sandy  Hook  and  The  Battery. 

Astronomic  Tide  of  October  1976 


43.  The  astronomic  tide  of  11-13  October  1976  was  used  to  calibrate  the 
original  model  of  the  Passaic  River  (CDM  1983b).  This  tide  was  also  used  to 
validate  the  changes  made  to  the  grid  for  the  Hackensack  River  Surge  Barrier 
study  (CDM  1988).  The  full  modified  grid  was  used  to  simulate  this  event, 
with  interior  forcing  from  observations  at  Sandy  Hook  (node  12)  and  The 
Battery  (node  62),  and  forcing  at  boundary  nodes  1-4  developed  from 
extrapolated  readings  from  tide  gages  in  the  area.  Since  the  interior 
boundary  conditions  were  also  specified,  forcing  at  these  nodes  dominates  the 
interior  solution  in  Newark  Bay  and  its  tributaries. 
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44.  Figure  5  shows  the  locations  of  observation  points.  Figures  6  and 
7  show  the  observed  and  predicted  tides  at  Sandy  Hook  and  The  Battery. 

Figures  8  through  23  each  show:  a)  observed  tidal  heights;  b)  the  results  of 
simulations  using  the  modified  grid  for  this  study;  and  c)  the  results  from 
the  Hackensack  River  Surge  Barrier  Study  (COM  1988).  These  figures  show  that 
the  model  agrees  well  with  observations,  and  that  at  most  locations  the 
agreement  is  at  least  as  good  as,  if  not  better  than,  the  results  of  the 
Hackensack  River  Surge  Barrier  study.  From  this  comparison,  we  conclude  that 
the  model  is  accurately  simulating  tidal  conditions  in  the  Newark  Bay  system. 
Sensitivity  of  forcing  from  observations  at  Sandy  Hook  (node  12)  and  The 
Battery  (node  62)  was  tested  by  rerunning  the  larger  grid  of  the  system  with 
forcing  at  only  the  exterior  nodes  1-4.  The  results  of  this  simulation  are 
shown  in  Figures  6  and  7,  respectively,  and  demonstrate  that  the  interior 
region  can  be  accurately  forced  using  boundary  conditions  specified  at  either 
sets  of  locations. 

45.  One  additional  simulation  was  performed  in  which  friction  values 
along  the  Hackensack  River  were  reduced  by  about  25  percent  to  determine 
whether  significant  changes  in  predicted  tidal  elevations  along  that  river 
would  be  seen.  The  results  indicated  that  little  change  (less  than  0.1  ft) 
was  found,  confirming  that  system  geometry  often  plays  a  more  important  role 
in  defining  the  model  than  does  bottom  friction. 

Historical  Hurricanes 


46.  In  the  Passaic  River  study  (CDM  1983b)  two  hurricanes  were 
simulated  to  validate  the  model.  These  events  were  not  simulated  in  the 
Hackensack  River  Surge  Barrier  study,  as  it  was  a  reconnaissance  study  of 
limited  scope.  For  the  current  study,  we  resimulated  these  two  events  using 
the  smaller  grid  of  the  system,  and  direct  forcing  at  nodes  12,  24  and  25 
using  observations  at  Sandy  Hook,  and  at  nodes  62  and  63  using  observations  at 
The  Battery. 

47.  The  hurricane  of  21  September  1938  was  simulated  using  50  hours  of 
observations  from  0000  hours,  20  September  1938.  Figure  24  and  25  show  the 
observed  storm  tides  at  Sandy  Hook  at  The  Battery,  respectively.  Figures  26 
and  27  show  comparisons  with  observed  high  and  low  tides  at  Elm  Park  and  East 
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Newark,  respectively.  These  figures  show  that  the  model  predicts  this  event 
accurately. 

48.  The  hurricane  of  12  September  1960  was  simulated  using  48  hours  of 
observations  from  0000  hours,  11  September  1960.  Figures  28  and  29  show  the 
observed  tides  at  Sandy  Hook  and  The  Battery,  respectively.  Figures  30 
through  32  show  comparisons  with  observed  tides  at  Elm  Park,  Fort  Hamilton  and 
Carteret,  respectively.  Figure  33  shows  predicted  and  observed  tidal  heights 
at  Kearny  Point.  These  figures  demonstrate  that  the  model  predicts  storm  tide 
heights  accurately  during  this  event. 

Historical  Northeasters 

49.  The  northeaster  of  24-26  November  1950  was  simulated  using  the 
smaller  grid,  and  direct  forcing  at  the  Sandy  Hook  and  The  Battery  boundary 
nodes.  The  event  was  simulated  from  0000  hours,  24  November  1950  for  52 
hours.  Figures  34  and  35  show  the  observed  tides  at  Sandy  Hook  and  The 
Battery,  respectively.  Figure  36  shows  a  comparison  with  observed  tidal 
heights  at  Elm  Park,  and  the  agreement  is  generally  very  good.  However,  it 
appears  from  Figure  36  that  the  gage  at  Elm  Park  topped  out  and  missed  two 
consecutive  peak  tides. 

50.  The  northeaster  of  6-7  November  1953  was  simulated  from  0000 
hours,  6  November  1953  for  48  hours.  Figure  37  and  38  show  the  observed  storm 
tide  at  Sandy  Hook  and  The  Battery,  respectively.  Figures  39  through  41  show 
comparisons  with  observations  at  Elm  Park,  East  Newark  and  Fort  Hamilton, 
respectively.  These  figures  further  validate  that  the  model  accurately 
simulates  this  northeaster  event. 

MHalloweenM  Northeaster  of  1991 

51.  During  30-31  October  1991,  the  "Halloween"  storm  hit  the  Newark  Bay 
area.  The  New  York  District  of  the  Corps  of  Engineers  recorded,  as  accurately 
as  possible  in  the  circumstances,  high  water  levels  in  the  afternoon  of 
Wednesday,  30  October  and  the  morning  of  Thursday,  31  October  1991.  We 
simulated  this  event  for  71  hours  from  0000  hours,  29  October  1991,  using  a 
grid  that  reflected  the  partial  dredging  in  the  Kills  and  Newark  Bay  to  this 
date. 
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52.  Figures  42  and  43  show  the  observed  tides  at  Sandy  Hook  and  The 
Battery,  respectively.  Figure  44  shows  predicted  and  observed  tidal  heights 
at  Ridgefield  Park  (Bergen  Generating  Station).  Table  1  compares  observations 
throughout  the  Newark  Bay/Passaic/Hackensack  River  system  with  simulation 
results.  The  table  shows  generally  agreement  throughout  the  system.  However, 
there  are  two  locations  along  the  Hackensack  River  where  observed  heads  are 
more  than  a  foot  lower  than  predicted.  It  is  interesting  to  note  the  close 
agreement  at  either  end  of  the  Hackensack  River,  and  both  higher  predicted  and 
observed  tidal  elevations  than  at  these  two  more  central  locations. 

Therefore,  either:  a)  there  is  some  systematic  error  in  the  readings  at  these 
observation  points,  b)  the  readings  represent  overbank  observations  where 
obstructions  in  the  flood  plain  reduce  water  levels  from  those  seen  in  the 
main  river ,  or  c)  the  model  does  not  predict  the  storm  tide  correctly.  Other 
simulations  suggest  that  the  model  accurately  predicts  storm  surge  events,  and 
it  is  unlikely  that  flood  plain  obstructions  close  to  the  main  channel  would 
produce  the  differences  in  water  surface  "observed".  We  therefore  suspect 
that  the  differences  are  produced  by  some  systematic  error  in  the  readings  at 
locations  near  nodes  206  and  210. 

53.  Finally,  to  determine  the  effect  of  dredging  on  tidal  heights  in 
the  Newark  Bay  system,  the  "Halloween"  Storm  was  simulated  using  the  "pre- 
1980"  grid.  The  results  showed  that,  at  most,  predicted  tidal  heights  for 
this  extreme  event  changed  by  less  than  0.02  feet.  In  addition,  we  simulated 
the  observed  hurricanes  of  September  1938  and  1960,  and  the  northeasters  of 
1950  and  1953,  using  the  "1991"  grid  without  sea  level  rise.  The  results 
showed  that  peak  surge  elevations  compared  within  0.15  feet,  and  that  the 
effect  of  the  dredging  was  to  slightly  alter  the  phasing  of  the  peak 
elevations,  but  by  less  than  1/2  hour.  Therefore,  we  conclude  that  the 
dredging  operations  in  the  Kills  and  Newark  Bay  has,  and  will  have,  little 
effect  on  circulation  and  water  levels  in  the  study  area. 
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Table  1 


Comparison  of  Observed  and  Simulated  Tidal  Elevations 
For  the  "Halloween"  Storm  of  October  1991 


I  Location 

Node 

Observed 

Height' 

Observed 

Time1 2 

Predicted 

Height' 

Predicted 

Time 

New  York  Harbor 

60 

6.2 

W  (2pm) 

6.02 

W  (2  pm) 

Caven  Point ,  New 
York  Harbor 

Passaic  River 

BASF  Plant, 

Kearny 

262 

7.62-7.94 

T  (debris) 

7.23 

T  (3:45  am) 

4th  St . ,  Newark 

171 

6.83 

7.43 

W  (3  pm) 

T  (debris) 

6.66 

7.29 

W  (3  pm) 

T  (4  am) 

Fire  Pier, 

Newark 

172 

7.75 

T  (debris) 

7.31 

T  (4  am) 

Page  St . , 
Lyndhurst 

216 

5.72-6.22 

6.25-7.6 

W  (3:30  pm) 

T  (debris) 

6.78 

7.37 

W  (3:30  pm) 

T  (3:45  am) 

Old  Boat  House, 
East  Rutherford 

219 

6.34 

T  (debris) 

7.39 

T  (3:45  am) 

Hackensack  River 

Daily  News, 

Kearny 

264 

>  7.1 

T  (77) 

7.22 

T  (3:30  am) 

Fence  Outlet, 
Ridgefield  Park 

207 

5.34 

5.94 

W  (3  pm) 

T  (7-8  am) 

W  (3  pm) 

T  (7:30  am) 

Gas  Station, 

South  Hackensack 

206 

4.95 

T  (7  am) 

6.06 

T  (7  am) 

Marina,  South 
Hackensack 

210 

5.79 

T  (debris) 

6.91 

T  (5:15  am) 

Hess,  Bogota 

211 

7.2 

T  (debris) 

6.94 

T  (5:15  am) 

1  Feet  NGVD 

2  W  =  Wednesday  afternoon,  30  October  1991 
T  =  Thursday  morning,  31  October  1991 


Figure  5. 

Gage  Locations  for 
October  1976  Tidal 
Elevation  Study 

(Note  Numbers  Identify 
Nodal  Points) 

Cartography  by 
EBAS00  Environmental 
Hay  1992 


Figure  6.  Tidal  heights  at  Sandy  Hook 
for  tide  of  11  “13  October  1976 
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PREDICTED  .  OBSERVED 


Figure  7.  Tidal  heights  at  the  Battery 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED 


Figure  8.  Tidal  heights  at  Point  No 
Point  for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (1908)  STUDY 


Figure  9.  Tidal  heights  at  Belleville 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  - CDM  (1988)  STUDY 


Figure  10.  Tidal  heights  at  East 
Rutherford  for  tide  of  11-13  October 

1976 
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PREDICTED  .  OBSERVED  -  CDM  (1988)  STUDY 


Figure  11.  Tidal  heights  at  Garfield 
for  tide  of  11  —  13  October  1976 
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PREDICTED  .  OBSERVED  CDM  (1988)  STUDY 
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Figure  13.  Tidal  heights  at  Belleville 
Tpke  for  tide  of  11-13  October  1976 
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Figure  14.  Tidal  heights  at  Carlstadt 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  - CDM  (1988)  STUDY 


Figure  15.  Tidal  heights  at  Mill  Creek 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  CDM  (1988)  STUDY 


Figure  16.  Tidal  heights  at  Cromakill 
Creek  for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (1988)  STUDY 


Figure  17.  Tidal  heights  at  Ridgefield 
Park  for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  CDM  (1988)  STUDY 


Figure  18.  Tidal  heights  at  Hackensack 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (1988)  STUDY 


Figure  19.  Tidal  heights  at  Port  Ivory 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (1988)  STUDY 


Figure  20.  Tidal  heights  at  South  Amboy 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (198B)  STUDY 


Figure  21.  Tidal  heights  at  Sayreville 
for  tide  of  11  —  13  October  1976 


51 


PREDICTED  .  OBSERVED  - CDM  (1988)  STUDY 


Figure  22.  Tidal  heights  at  Keasbey 
for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (1988)  STUDY 


Figure  23.  Tidal  heights  at  New 
Brunswick  for  tide  of  11-13  October  1976 
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PREDICTED  .  OBSERVED  -  CDM  (1988)  STUDY 


Figure  24.  Observed  tidal  heights  at 
Sandy  Hook  for  hurricane  of  21 
September  1938 
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TIME  IN  HOURS 


Figure  25.  Observed  tidal  heights  at 
The  Battery  for  hurricane  of  21 
September  1938 
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TIME  IN  HOURS 


Figure  26.  Tidal  heights  at  Elm  Park 
for  hurricane  of  21  September  1938 
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PREDICTED  ‘  OBSERVED 


Figure  27.  Tidal  heights  at  East  Newark 
for  hurricane  of  21  September  1938 


57 


PREDICTE 


Figure  28.  Observed  tidal  heights  at 
Sandy  Hook  for  hurricane  of  12 
September  1960 
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TIME  IN  HOURS 


Figure  29.  Observed  tidal  heights  at 
The  Battery  for  hurricane  of  12 
September  1960 
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PREDICTED  .  OBSERVED 


Figure  31.  Tidal  heights  at  Fort 
Hamilton  for  hurricane  of  12 
September  1960 
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PREDICTED  .  OBSERVED 


Figure  32.  Tidal  heights  at  Carteret 
for  hurricane  of  12  September  1960 
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PREDICTED  .  OBSERVED 


Figure  33.  Tidal  heights  at  Kearny 
Point  for  hurricane  of  September  1960 
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Figure  34.  Observed  tidal  heights  at 
Sandy  Hook  for  northeaster  of  24-26 
November  1950 
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Figure  35.  Observed  tidal  heights  at 
The  Battery  for  northeaster  of  24-26 
November  1950 
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Figure  36.  Tidal  heights  at  Elm  Park 
>r  northeaster  of  24-26  November  1950 
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PREDICTED  .  OBSERVED 


Figure  37.  Observed  tidal  heights  at 
Sandy  Hock  for  northeaster  of  6-8 
November  1953 
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Figure  38.  Observed  tidal  heights  at 
The  Battery  for  northeaster  of  6-8 
November  1953 
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PREDICTED  .  OBSERVED 


Figure  40.  Tidal  heights  at  East  Newark 
for  northeaster  of  6-8  November  1953 
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PREDICTED  .  OBSERVED 


Figure  41.  Tidal  heights  at  Fort 
Hamilton  for  northeaster  of  6-8 
November  1953 
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Figure  42.  Observed  tidal  heights  at 
Sandy  Hook  for  the  "Halloween"  storm 
of  October  1991 
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TIME  IN  HOURS 


Figure  43.  Observed  tidal  heights  at 
The  Battery  for  the  "Halloween”  storm 
of  October  1991 
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TIME  IN  HOURS 


Figure  44.  Tidal  heights  at  Ridgefield 
Park  for  the  "Halloween”  storm  of 
October  1991 
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PARS  XV:  COINCIDENCE  ANALYSIS 


Introduction 

54.  Astronomic  tides,  meteorological  conditions,  and  river  flows  cause 
fluctuations  of  water  level  in  the  modeled  system.  Astronomical  tides  are 
independent  of  the  other  two,  although  tidally- induced  water  surface 
elevations  can  combine  with  either  or  both  meteorological  and  fluvial  events 
to  cause  very  high  water  elevations.  Since  the  greatest  effect  would  occur 
when  meteorological  and  fluvial  events  reinforce  each  other,  information 
regarding  their  joint  relationship  is  useful  in  examining  the  potential 
occurrence  of  high  water  levels. 

Approach  and  Data  Preparation 

55.  Several  approaches  were  taken  to  investigate  the  relation  between 
river  flows  and  storm  surges  (storm  tide  minus  astronomic  tide).  The 
approaches  were  based  on  observed  water  surface  elevations  and  predicted  tides 
at  The  Battery  in  New  York  Harbor  used  in  conjunction  with  Passaic  River  flows 
measured  at  Little  Falls,  New  Jersey  (785  square  miles).  Approaches  included 
cross  correlations  and  linear  regressions  between  the  data  serieB,  examination 
of  the  frequency  of  joint  occurrence  of  large  flow  volume  and  surge  height 
events,  and  scatter  plots  of  corresponding  flow  and  surge  events. 

56.  Time  series  of  hourly  observed  water  level  at  The  Battery  were 
obtained  through  WES  from  the  National  Ocean  Service  (NOS)  office  in 
Rockville,  Maryland.  A  standard  tidal  prediction  model  and  ancillary  data 
files  were  obtained  directly  from  the  same  office.  Daily  river  flow  data 
measurements  were  obtained  from  the  New  Jersey  USGS  office.  River  flow  data 
spanned  the  interval  of  1898  to  1990,  and  observed  water  level  data  extended 
from  1920  to  1991. 

57.  Daily  values  cf  water  level  elevations  caused  by  influences  other 
than  tides  (which  will  be  referred  to  as  surge  heights)  were  obtained  in  the 
following  way.  The  tide  prediction  package  was  used  to  produce  time  series  of 
hourly  predicted  tide  heights.  Predicted  tides  took  into  account  the  increase 
in  mean  sea  level  from  1920  to  the  present  using  Epoch  values  provided  by  NOS. 
Predicted  tide  values  were  subtracted  from  the  corresponding  hourly  observed 
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values  to  obtain  hourly  surge  heights,  and  the  maximum  daily  surge  heights 
were  extracted  from  the  respective  24-hour  values  for  each  day  of  the  time 
series. 

58.  Time  series  of  yearly  maximum  values  of  river  flow  and  surge  height 
were  extracted  from  the  respective  time  series  of  daily  maximum  values.  The 
smallest  value  in  each  of  the  time  series  of  annual  maximum  values  was  taken 
as  an  indication  of  the  magnitude  of  the  respective  one-year  event.  The 
values  were  1.85  feet  for  surge  height  and  3040  cubic  feet  per  second  (cfs) 
for  river  flow.  Since  1.85  feet  was  0.5  feet  less  than  the  next-larger  yearly 
maximum  of  2.35  feet,  and  because  there  were  nearly  2,600  values  larger  than 
1.85  feet,  the  value  of  2.35  feet  was  taken  as  an  estimate  of  the  one-year 
surge  event. 

59.  The  time  series  of  observed  water  levels  at  the  Battery  contained 
several  gaps,  which  produced  corresponding  gaps  in  the  surge  height  time 
series.  Surge  height  time  series  gaps  of  five  days  or  less  were  bridged  by 
inserting  a  linear  fit  between  values  at  the  ends  of  the  gaps.  No  attempt  was 
made  to  bridge  larger  gaps.  Gaps  were  from  1921  to  1927,  January  1958  to  May 
1958,  May  1973  to  May  1977,  and  September  1981  to  October  1981.  Six  time 
series  resulted  after  considering  gaps  and  start /end  dates  for  tidal  epochs: 
1927-1952,  1953-1958,  1958-1961,  1962-1973,  1977-1980,  and  1981-1991.  This 
resulted  in  21,171  daily  observations  for  the  analysis. 

Comparison  of  Water  Surface  Elevations  at  Kearnv  Point  and  The  Battery 

60.  The  time  series  of  observed  water  elevations  is  from  measurements  at 
The  Battery,  which  is  close  to  the  area  of  interest,  namely  Newark  Bay  and  the 
Passaic  and  Hackensack  rivers.  In  order  to  establish  that  both  the  timing  and 
magnitude  of  water  surface  elevation  at  The  Battery  were  representative  of 
conditions  at  Kearny  Point,  water  surface  elevations  at  these  two  locations 
were  plotted  against  the  same  time  axis  for  three  historical  cases:  (1)  the 
northeaster  of  24-26  November  1950  (Figure  45),  (2)  the  hurricane  of  12 
September  1960  (Figure  46),  and  (3)  tides  during  11-13  October  1976  (Figure 
47). 

61.  For  these  three  types  of  events,  the  three  figures  show  close 
agreement  between  timing  and  magnitude  of  water  surface  elevation  at  The 
Battery  and  Kearny  Point.  Therefore,  using  water  surface  elevations  obtained 
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at  The  Battery  for  examining  coincidence  between  surges  and  Passaic  River 
flows  will  yield  results  applicable  to  the  vicinity  of  Kearny  Point. 

Suroe  Height  and  River  Plow  Cross  Correlation 

62.  Lagged  cross  correlation  between  surge  height  and  river  flow  showed 
that  peak  correlation  (*0.22)  occurred  between  surge  height  on  a  given  day  and 
river  flow  measured  two  days  later  (Figure  48).  That  is,  changes  in  river 
flow  tended  to  lag  surge-caused  changes  in  water  surface  elevation  by  about 
two  days.  Positive  correlation  means  that  increases  in  surge  height  are 
accompanied  by  increases  (rather  than  decreases)  in  river  flow.  The  cross 
correlation  coefficient  of  0.22  represents  a  weak  correlation  (where  a  value 
of  1.0  implies  complete  dependence  and  a  value  of  0.0  implies  complete 
independence) . 

63.  The  degree  of  correlation  begins  to  increase  when  river  flow  leads 
surge  height  by  one  day.  The  correlation  peaks  when  surge  height  leads  by  two 
days,  and  shows  a  slow  decline  thereafter.  This  indicates  that  river  flow 
events  tend  to  occur  later  than  surge  events  and  to  decrease  more  slowly. 

64.  Linear  regression  was  calculated  between  surge  height  values  and 
river  flows  two  days  later,  which  represents  the  timing  that  would  yield  the 
best  fit.  The  resulting  regression  relation  could  explain  only  4.4  percent  of 
the  variance  of  the  river  flow  when  using  the  surge  height  time  series  as  the 
independent  variable.  Surge  height  is  therefore  a  poor  predictor  of  river 
flow,  which  is  to  be  expected  in  view  of  the  relatively  weak  correlation 
between  the  two. 

Coincidence  Between  Suroe  Height  and  River  Flow 

65.  Coincidence  was  examined  by  comparing  daily  values  of  surge  height 
and  river  flow  with  respective  threshold  values  of  2.35  feet  and  3040  cfs. 
Counted  were  the  number  of  times  (1)  both  daily  values  exceeded  their 
thresholds,  (2)  neither  exceeded  its  threshold,  (3)  surge  height  exceeded  its 
threshold,  but  river  flow  did  not,  and  (4)  river  flow  exceeded  its  threshold, 
but  surge  height  did  not.  Comparisons  were  made  for  seven  cases  that  ranged 
from  surge  height  leading  river  flow  by  three  days  to  river  flow  leading  surge 
height  by  three  days. 
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66.  The  results  (Table  2)  indicate  greater  coincidence  when  surge  leads 
river  flow,  which  reflects  the  findings  of  the  correlation  analysis. 
Coincidence  with  surge  height  leading  was  approximately  double  that  with  river 
flow  leading.  The  largest  count  of  both  series  exceeding  thresholds,  192, 
occurred  with  surge  height  leading  by  two  days.  A  count  of  192  represents 
only  0.9  percent.  In  contrast,  both  series  were  jointly  less  than  threshold 
approximately  89  percent  of  the  time.  The  analysis  shows  that,  given 
relatively  moderate  event  threshold  values,  surge  height  and  river  flow  events 
coincide  only  a  small  fraction  of  the  time. 

67.  The  relation  between  water  level  and  river  flow  was  further  examined 
for  five  historical  storm  events  by  plotting  predicted  tide  height  and 
observed  water  level  at  The  Battery,  along  with  Passaic  River  flow,  for  the 
time  interval  from  three  days  preceding  the  event  to  seven  days  after.  The 
selected  storm  events  were  21  September  1938  (Figure  49),  25  November  1950 
(Figure  50),  12  September  1960  (Figure  51),  10  November  1977  (Figure  52),  and 
7  April  1984  (Figure  53).  Maximum  surge  heights  range  from  approximately  four 
to  eight  feet,  and  maximum  river  flows  from  approximately  4,000  to  17,000  cfs. 
The  plots  show  that,  for  these  selected  events,  the  river  flow  peaks  within 
about  one  to  three  days  after  a  peak  surge.  However,  the  plots  also  show  that 
the  magnitudes  of  surge  and  flow  peaks  are  not  necessarily  directly  related. 
That  is,  a  large  surge  peak  may  be  accompanied  by  either  large  or  small  flow 
peaks,  and  vice-versa.  For  example,  a  surge  of  about  seven  feet  and  flow  of 
4,000  cfs  occurred  during  the  1950  storm,  while  a  surge  of  four  feet  and  flow 
of  17,000  cfs  occurred  during  the  1984  storm.  The  indication  is  that,  for 
major  events,  peaks  in  surge  and  flow  tend  to  accompany  each  other  within 
approximately  72  hours,  but  that  there  is  no  consistent  relation  between  their 
magnitudes . 


Scatter  Diagram  of  Surge  Height  and  River  Flow  Values 

68.  In  order  to  investigate  the  relation  between  surge  height  and  river 
flow  visually,  a  scatter  diagram  was  plotted  with  surge  height  values  on  the 
horizontal  axis  and  river  flow  values  on  the  vertical  axis.  Only  those  points 
for  which  surge  height  exceeded  0.5  feet  and  river  flow  exceeded  1,000  cfs 
jointly  were  plotted.  Note  that  for  each  identified  surge  height  value,  only 
the  single  maximum  river  flow  within  the  succeeding  72  hours  was  selected. 
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69.  The  scatter  diagram  (Figure  54)  shows  the  absence  of  any  strong 
connection  between  surge  height  and  river  flow.  Points  tend  to  be  distributed 
primarily  along  the  vertical  axis  of  the  diagram.  This  pattern  indicates 
that,  for  a  given  surge  height,  there  is  a  wide  range  of  possible  river  flows, 
with  most  of  them  being  relatively  small.  To  a  lesser  extent,  there  is  a 
range  of  surge  heights  for  a  given  river  flow.  Close  relationship  between 
surge  height  and  river  flow  would  have  been  indicated  by  points  being 
distributed  about  some  diagonal  from  the  origin. 

70.  The  scatter  plot  reinforces  the  results  of  correlation  and 
regression  calculations  and  the  findings  of  the  coincidence  analysis.  When 
only  moderate  thresholds  are  specified,  surge  height  events  and  river  flow 
events  are  only  weakly  related,  and  the  occurrence  of  one  is  not  a  conclusive 
predictor  of  either  the  occurrence  or  the  magnitude  of  the  other.  However, 
there  is  some  indication  (based  upon  plots  of  the  five  selected  storm  events) 
of  greater  likelihood  of  joint  occurrence  of  surge  and  flow  increases  during 
major  storm  events,  even  though  their  magnitudes  may  not  be  directly  related. 

Stage  Frequency  Statistics  for  Surge  Heights 

71.  The  distribution  of  water  surface  elevations  at  The  Battery  was 
examined  for  peak  values  occurring  within  72  hours  prior  to  Passaic  River 
flows  exceeding  6,000  cfs.  This  is  the  threshold  flow  when  the  proposed 
diversion  tunnel  begins  to  function.  The  analysis  was  performed  to  examine 
what  magnitude  of  storm  tide  might  occur  when  the  diversion  tunnel  was 
operating. 

72.  According  to  the  record  at  the  Little  Falls  gage  on  the  Passaic 
River,  230  of  the  21,174  days  examined  had  flows  that  exceeded  6000  cfs.  For 
each  of  these  230  records,  the  maximum  tidal  height  in  the  previous  72  hours 
was  identified  from  the  tide  gage  record  at  The  Battery.  These  maximum  tidal 
heights  were  then  ranked  and  plotted  in  Figure  55  as  a  cumulative  frequency 
distribution.  The  plot  is  a  conditional  cumulative  frequency  (or  exceedance) 
distribution.  The  plot  presents  both  the  percent  daily  exceedance  and  the 
return  period  for  storm  tides  that  are  coincident  with  river  flows  exceeding 
6000  cf 8 .  Also  plotted  on  Figure  55  is  the  stage  frequency  distribution  for 
all  river  flows.  Comparison  between  these  two  distributions  shows  that  a 
given  tidal  elevation  coincident  with  Passaic  River  flow  greater  than  6000  cfs 
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occurs  much  less  frequently  than  when  no  river  flow  condition  is  included. 

For  example,  a  tidal  elevation  of  six  feet  independent  of  river  flow  has  a 
return  period  of  2.6  years.  However,  a  tidal  elevation  of  six  feet  coincident 
with  flow  greater  than  6000  cfs  has  a  return  period  of  19.3  years,  which 
corresponds  to  only  13.5  percent  as  often.  Examination  of  Figure  55  shows 
that  the  elevation  of  6.15  feet  (the  maximum),  in  combination  with  flow 
greater  than  6000  cfs,  is  equalled  or  exceeded  once  in  the  21,174  days  of 
record,  or  0.005  percent  of  days  (return  period  of  58  years).  An  elevation  of 
2.15  feet  (a  typical  high  tide  elevation)  is  exceeded  only  1.1  percent  of  the 
days  (return  period  of  0.25  years).  Tidal  heights  of  4.2  feet  and  5.9  feet 
are  coincident  with  Passaic  River  flows  exceeding  6000  cfs  for  return 
intervals  of  one  year  and  10  years,  respectively.  It  is  possible  to  estimate 
longer  return  period  events  (say,  100  years)  by  extrapolating  the  conditional 
distribution  in  Figure  55. 


Condition t 


Low  Models 


73.  The  weak  coincidence  between  river  flow  and  storm  surges  may  at 
first  suggest  that  the  two  "events"  be  considered  to  be  independent.  However, 
the  results  of  the  analysis  also  shows  that,  while  the  coincidence  is  weak, 
the  events  are  not  completely  independent.  One  question  that  is  asked  is  what 
type  of  tidal  boundary  condition  should  be  specified  when  using  a  river  model 
to  simulate  significant  flow  events. 

74.  If  the  probability  of  a  specific  river  flow  can  be  statistically 
determined,  then  the  overall  statistic  is  not  altered  by  selecting  the  average 
downstream  tidal  boundary  condition  for  that  flow  condition.  Consequently, 
the  daily-maximum  surge  elevations  within  72  hours  of  observed  river  flows 
were  averaged  for  specified  ranges  of  river  flows,  and  the  results  summarized 
in  Table  3.  Table  3  lists  the  average  surge  (without  tide)  elevation;  (1)  on 
the  same  day  as  the  river  flow,  (2)  one  day  prior,  (3)  two  days  prior,  and 

(4)  the  maximum  of  the  72  hours  prior  to  the  river  flow. 

75.  Table  3  shows  that  a  storm  surge  of  1.5-  to  2-feet  will  occur  on 
average  within  72  hours  of  various  flow  ranges  in  the  Passaic  River.  However, 
it  is  very  difficult  to  analyze  the  shape  of  the  storm  (duration  and 
"peakiness")  that  would  occur.  Overall,  we  recommend  that  a  river  discharge 
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model  of  the  Passaic  River  specify  a  tidal  boundary  that  is  the  combination  of 
an  average  tide  plus  a  1. 5-2.0  feet  surge. 

River  flows  for  Storm  Tide_simulatioos 

76.  The  weak  coincidence  between  river  flows  and  storm  surges  also  leads 
to  the  question  of  what  river  flow  to  specify  when  simulating  a  storm  tide 
event.  This  question  was  addressed  in  a  similar  fashion  as  before  by 
investigating  the  river  flows  that  occurred  within  the  72  hours  following  an 
observed  storm  surge.  The  data  were  analyzed  to  determine  the  average  flow 
observed  for  a  surge  of  a  specified  height.  The  results  are  presented  in 
Table  4a.  In  Table  4b,  Passaic  River  flows  associated  with  historical  storm 
tide  events  are  tabulated  to  illustrate  flows  observed  during  these  events. 

77.  Table  4  lists  the  average  flows;  (1)  on  the  same  day  as  the  storm 
surge,  (2)  one  day  following,  (3)  two  days  following,  and  (4)  the  maximum  of 
the  72  hours  following  the  surge  event.  The  results  show  little  variation  in 
average  flows  between  surge  heights,  and  indicate  that  when  a  surge  event 
occurs,  a  river  flow  in  the  range  1000-2000  cfs  might  be  expected. 

78.  To  add  "conservatism"  to  the  subsequent  storm  tide  analyses,  we 
decided  to  the  specify  the  one-year  flows  in  each  of  the  major  tributaries. 

On  the  Passaic  River,  this  results  in  a  flow  in  excess  of  5000  cfs,  or  nearly 
three  times  the  average  flow  that  might  be  expected.  Additional  coincidence 
analyses  have  shown  that  this  flow  will  occur  relatively  infrequently  with 
storm  tide  events. 
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Table  2 


Tallies  of  coincidence  between 
surge  height  and  river  flow  events 


***********************  COINCIDENCE  ANALYSIS  ********************** 

THE  TOTAL  NUMBER  OF  POINTS  USED  WAS  21171 

TIME  SERIES  SURGE. DLY  :  MEAN  «  .95  STD.  DEV.  *=  .72 

TIME  SERIES  RIVER. DLY  :  MEAN  *  1096.90  STD.  DEV.  «*  1330.02 

SURGE  HEIGHT  THRESHOLD  VALUE  IS  2.35  FEET 
RIVER  FLOW  THRESHOLD  VALUE  IS  3040.00  CFS 

THRESHOLD  VALUES  WERE  SPECIFIED  BY  USER 

WHEN  COMPARING  SURGE  HEIGHT  WITH  RIVER  FLOW  3  DAYS  LATER: 

BOTH  SERIES  EXCEED  THRESHOLDS  179  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  640  TIMES 

SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1432  TIMES 

BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18917  TIMES 

WHEN  COMPARING  SURGE  HEIGHT  WITH  RIVER  FLOW  2  DAYS  LATER: 

BOTH  SERIES  EXCEED  THRESHOLDS  192  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  627  TIMES 

SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1419  TIMES 

BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18930  TIMES 

WHEN  COMPARING  SURGE  HEIGHT  WITH  RIVER  FLOW  1  DAY  LATER: 

BOTH  SERIES  EXCEED  THRESHOLDS  191  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  628  TIMES 

SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1420  TIMES 

BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18929  TIMES 

WHEN  COMPARING  SURGE  HEIGHT  WITH  RIVER  FLOW  ON  THE  SAME  DAY: 

BOTH  SERIES  EXCEED  THRESHOLDS  137  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  682  TIMES 
SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1474  TIMES 
BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18875  TIMES 

WHEN  COMPARING  RIVER  FLOW  WITH  SURGE  HEIGHT  1  DAY  LATER: 

BOTH  SERIES  EXCEED  THRESHOLDS  91  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  728  TIMES 
SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1520  TIMES 
BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18829  TIMES 

WHEN  COMPARING  RIVER  FLOW  WITH  SURGE  HEIGHT  2  DAYS  LATER: 

BOTH  SERIES  EXCEED  THRESHOLDS  94  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  725  TIMES 
SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1517  TIMES 
BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18832  TIMES 

WHEN  COMPARING  RIVER  FLOW  WITH  SURGE  HEIGHT  3  DAYS  LATER: 

BOTH  SERIES  EXCEED  THRESHOLDS  88  TIMES 

SURGE  EXCEEDS  THRESHOLD  AND  FLOW  IS  LESS  THAN  THRESHOLD  731  TIMES 
SURGE  IS  LESS  THAN  THRESHOLD  AND  FLOW  EXCEEDS  THRESHOLD  1523  TIMES 
BOTH  SERIES  ARE  LESS  THAN  THRESHOLDS  18826  TIMES 
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Table  3 


Average  Surges1  Coincident  with  Given  Passaic  River  Flows 


Surae 

( feet » 

Flow  Ranae 

2  days 

1  day 

Same 

Maximum 

(cfs) 

Prior 

Prior 

Day 

72  Hour 8 

0.- 

1000. 

.90 

.90 

.88 

1.19 

1000.- 

2000. 

.99 

1.00 

1.02 

1.38 

2000.- 

3000. 

1.05 

1.04 

1.06 

1.46 

3000.- 

4000. 

1.04 

1.02 

1.12 

1.46 

4000.- 

5000. 

1.01 

.98 

1.16 

1.47 

5000.- 

6000. 

1.07 

1.16 

1.21 

1.55 

6000.- 

7000. 

1.28 

1.21 

1.25 

1.61 

7000.- 

8000. 

1.15 

1.03 

1.16 

1.44 

8000.- 

9000. 

1.13 

1.27 

1.40 

1.64 

9000.- 

10000. 

1.18 

1.33 

1.42 

1.65 

10000.- 

11000. 

1.65 

1.61 

1.51 

1.85 

11000.- 

12000. 

1.36 

1.63 

1.64 

1.91 

12000.- 

13000. 

1.42 

1.32 

1.27 

1.42 

13000.- 

14000. 

1.60 

1.40 

1.64 

1.75 

14000.- 

15000. 

1.91 

1.32 

1.49 

2.04 

15000.- 

16000. 

— 

no  flows 

in  record 

16000.- 

17000. 

2.16 

2.23 

1.91 

2.23 

17000.- 

18000. 

— - - 

no  flows 

in  record 

— 

18000.- 

19000. 

1.32 

1.11 

1.05 

1.32 

>  19000.  -  no  flows  in  record  - 


Note :  'Surge  in  storm  tide  minus  astronomic  tide. 
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Table  4a 


Average  Passaic  River  Flows  for  Given  Surge  Heights 


Surge1 


Passaic  River  Flow  (cfsl 


height 

2  days 

1  day 

Same  Maximum 

(ft) 

Later 

Later 

Dav  72  Hours 

-2.50 

101. 

116. 

108. 

116. 

-2.00 

1522. 

1751. 

1836. 

1858. 

-1.50 

1335. 

1351. 

990. 

1459. 

-1.00 

1213. 

1247. 

1095. 

1400. 

-.50 

1126. 

1207. 

1174. 

1339. 

.00 

888. 

942. 

938. 

1053. 

.50 

889. 

931. 

946. 

1053. 

1.00 

1193. 

1195. 

1194. 

1373. 

1.50 

1145. 

1125. 

1123. 

1307. 

2.00 

1356. 

1293. 

1276. 

1527. 

2.50 

1673. 

1507. 

1507. 

1855. 

3.00 

1876. 

1518. 

1449. 

2013. 

3.50 

1717. 

1324. 

1286. 

1837. 

4.00 

1715. 

1187. 

1150. 

1805. 

4.50 

901. 

919. 

910. 

1002. 

5.00 

1782. 

1561. 

1534. 

1818. 

5.50 

2602. 

1002. 

1018. 

2629. 

6.00 

1356. 

363. 

387. 

1356. 

6.50 

943. 

148. 

161. 

943. 

7.00 

- no 

surges 

in  record 

— 

7.50 

— -  no 

surges 

in  record 

- — 

8.00 

- no 

surges 

in  record 

— 

8.50 

1000. 

296. 

348. 

1000. 

Note: 

Surge  in  storm  tide 

minus 

astronomic 

tide. 

Table 

4b. 

Passaic  River  Flows  for  Selected  Storm  Events 


Passaic  River  Flow  (cfs 


Storm 

2  days 

Event 

Later 

21  Sept  1938 

9,140 

25  Nov  1950 

3,920 

12  Sept  1960 

5,040 

10  Nov  1977 

5,310 

7  Apr  1984 

13,200 

1  day 

Same 

Maximum 

Later 

Dav 

72  Hours 

7,020 

4,590 

9,140 

4,130 

1,000 

4,130 

4,910 

2,190 

5,040 

7,130 

8,960 

8,960 

16,000 

18,000 

18,000 
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Figure  45.  Tidal  heights  at  Kearny 
Point  and  The  Battery  for  northeaster 
of  24-26  November  1950 
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-  PREDICTED  .  OBSERVED 

(KEARNY  POINT)  (THE  BATTERY) 


Figure  46.  Tidal  heights  at  Kearny 
Point  and  The  Battery  for  hurricane 
of  12  September  1960 
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-  PREDICTED  .  OBSERVED 

(KEARNY  POINT)  (THE  BATTERY) 


Figure  47.  Tidal  heights  at  Kearny 
Point  and  The  Battery  for  tide  of 
11-13  October  1976 


87 


-  PREDICTED  .  OBSERVED  -  OBSERVED 

(KEARNY  POINT)  (KEARNY  POINT)  (THE  BATTERY) 


CROSS  CORRELATION 


ZP 


TIME  DIFFERENCE  (DATS) 


Figure  48. 


Cross  correlation  between  surge  height  at 
The  Battery  and  Passaic  River  flow  measured 
at  the  Little  Falls,  New  Jersey  gage. 


CD 


Figure  49.  Predicted  tide  height  and 
observed  water  level  at  The  Battery  and 
flow  of  the  Passaic  River  measured  at 
the  Little  Falls,  New  Jersey  gage  during 
the  21  September  1938  storm  event 
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Figure  50.  Predicted  tide  height  and 
observed  water  level  at  The  Battery  and 
flow  of  the  Passaic  River  measured  at 
the  Little  Falls,  New  Jersey  gage  during 
the  25  November  1950  storm  event 


os  —  >  wo: 
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PREDICTED  TIDE  HT. 


Figure  51.  Predicted  tide  height  and 
observed  water  level  at  The  Battery  and 
flow  of  the  Passaic  River  measured  at 
the  Little  Falls,  New  Jersey  gage  during 
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Figure  52.  Predicted  tide  height  and 
observed  water  level  at  The  Battery  and 
flow  of  the  Passaic  River  measured  at 
the  Little  Falls,  New  Jersey  gage  during 
the  10  November  1977  storm  event 
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PREDICTED  TIDE  HT. 


Figure  53.  Predicted  tide  height  and 
observed  water  level  at  The  Battery  and 
flow  of  the  Passaic  River  measured  at 
the  Little  Falls,  New  Jersey  gage  during 
the  7  April  1984  storm  event 
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PREDICTED  TIDE  HT. 


RIVER  FLOW  (CFS) 


20000 


15000 


10000 


5000 


0 


0  2  4  6  8  10 

SURGE  HEIGHT  (FEET) 


Figure  54.  Values  of  surge  height  greater  than  0.5  feet 

and  maximum  river  flow  greater  than  1,000  cfs 
within  the  following  72  hours. 
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Figure  55.  Stage-frequency  curves  an 
return  periods  for  water  surface 
elevation  at  The  Battery  and  Passaic 
River  flows  at  the  Little  Falls  gage 
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Part  V:  STATISTICAL  PROCEDURE 


Joint  Probability  Analysis 

79.  The  Joint  Probability  Method  (JPM)  was  used  to  develop  hurricane  and 
northeaster  surge  frequencies  for  previous  studies  of  the  area  (CDM  1983a, 
1983b,  and  1988).  The  method  involves  the  parameterization  of  the  storm 
characteristics  relevant  to  surge  formation,  which  includes:  (1)  central 
pressure  deficit,  (2)  radius  to  maximum  winds,  (3)  forward  speed,  (4)  travel 
direction,  and  (5)  latitude  at  which  the  storm  meets  the  coast.  The 
probability  distribution  of  each  of  these  parameters  is  derived  from  the 
observations  of  historical  hurricanes  and  northeasters.  Other  factors  which 
might  influence  peak  surge  elevation  are  the  tidal  amplitude,  and  the  temporal 
relationship  (lag)  between  tidal  periodicity  and  storm  occurrence.  Assuming 
that  the  storm  characteristics  are  statistically  independent,  and  that  the 
observed  distributions  are  accurate  indicators  of  the  underlying  parameter 
frequency  distributions,  the  frequency  of  occurrence  of  a  given  storm  is  the 
product  of  the  probability  of  occurrence  of  each  of  its  defining  parameter 
values. 

80.  Three  models  were  used  to  estimate  the  surge  elevation  along  the 
length  of  the  Passaic  and  Hackensack  Rivers,  which  would  occur  in  response  to 
a  synthetic  storm  defined  by  some  combination  of  parameter  values.  These 
models  are: 

1.  Wind  Field  Models  -  simulation  of  hurricane  or  northeaster  forcing 
on  the  ocean  surface. 

2.  Offshore  Surge  Model  -  two-dimensional,  hydrodynamic  model  of 
surge  generation  over  the  continental  shelf. 

3.  Inshore  Surge  Model  (DEM)  -  hydrodynamic  model  of  surge 
propagation  within  bays,  harbors  and  waterways. 

81.  The  JPM  was  used  to  develop  252  synthetic  hurricanes  and  43 
synthetic  northeasters.  These  295  total  synthetic  storm  events  were  then 
convolved  with  five  tide  types  and  50  half-hour  lags.  This  produced  over 
70,000  combinations,  each  with  a  probability  of  occurrence,  which  were  used  to 
develop  stage- frequency  curves  at  various  locations  in  the  study  area. 

Details  of  the  method  can  be  found  in  CDM  (1983a  and  1983b). 
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Modified  Bootstrap  Method 


82.  Rather  than  simulate  the  295  synthetic  storm  events,  and  then 
convolve  the  results  with  the  five  tidal  types  and  50  half-hour  lags  (to 
produce  295x5x50  =  73,750  "events"),  the  approach  used  for  the  Hackensack 
River  Surge  Barrier  Study  (COM  1988)  was  selected.  The  Hackensack  study  used 
a  smaller  subset  of  the  combinations  which  approximated  the  original  stage- 
frequency  distributions  (Prater  and  Butler  1986).  A  modified  version  of  the 
"Bootstrap"  method  (Efron  1982)  was  used  to  estimate  the  error  introduced  by 
simulating  the  more  than  70,000  events,  which  were  used  to  generate  the  stage- 
frequency  distributions  in  the  Passaic  River  Study  (CDM  1983b). 

83.  The  modification  to  the  "Bootstrap"  method  was  to  "select"  a  few 
events  over  the  range  of  return  periods  of  interest,  rather  than  to  choose 
many  at  random.  The  range  of  interest  in  CDM  (1988)  was  return  periods  of  10- 
100  years.  The  method  involves  the  selection  of  a  set  of  events 
representative  of  the  population  frequency  distribution  (Prater  et  al.  1988). 
The  percent  probability  of  occurrence  of  a  given  peak  surge  elevation  may  be 
obtained  from  a  stage-frequency  distribution  plot.  If  a  set  of  elevations  is 
selected  from  a  stage-frequency  distribution  curve,  the  horizontal  axis  (the 
probability  of  exceedance)  may  be  divided  into  ranges  that  equal  the 
probability  of  occurrence  for  each  specified  elevation. 

Storm  Selection 


84.  In  the  Hackensack  Surge  Barrier  Study  (CDM  1988),  24  synthetic  storm 
events  were  selected  from  the  stage- frequency  curve  at  node  44  at  the  mouth  of 
Newark  Bay,  with  the  majority  selected  for  return  periods  of  10-100  years. 

The  adequacy  of  this  number  of  events  was  tested  by  choosing,  from  the  24 
events  originally  selected,  three  sets  of  eight  events  and  two  sets  of  12 
events,  and  examining  how  the  sample  number  matched  the  stage- frequency 
distributions  at  various  locations  in  the  system  developed  from  the  full  set 
of  events.  The  results  (an  example  is  shown  in  Figure  56)  showed  that 
convergence  to  the  complete  distribution  was  obtained  as  more  events  were 
included  in  the  sample.  When  all  24  selected  events  were  used,  the  agreement 
between  the  full  distribution  and  the  bootstrap  distribution  was  acceptable. 
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85.  For  the  Passaic  River  Basin  Flood  Protection  study,  however,  we 
decided  to  use  more  events  to  fill  in  the  distribution  for  return  periods  of 
less  than  10  years  and  for  return  periods  of  greater  than  100  years.  The 
additional  events  were  selected  by  specifying  storm  tide  elevations  from  the 
stage-frequency  distribution  for  node  44  at  the  mouth  of  Newark  Bay,  and 
identifying  the  surge/tide/lag  combinations  from  the  more  than  70,000  events 
that  produced  total  maximum  tidal  surge  elevations  within  0.001  feet  of  the 
specified  value.  From  the  number  of  events  available  at  each  specified 
elevation,  a  single  event  combination  was  selected  for  the  study.  The  events 
were  selected  to  provide  a  mix  of  hurricanes  and  northeasters,  and  to 
distribute  these  two  types  of  storm  events  over  the  range  of  surge  elevations 
specified.  Table  5  lists  the  storm  events  selected.  The  nine  new  storms 
selected  for  this  study  are  indicated  with  an  asterisk.  Figure  57  shows  the 
stage-frequency  distributions  at  node  44  at  the  mouth  of  Newark  Bay,  and  the 
33  event  combinations  selected. 

86.  The  observed  range  of  elevation  differences,  from  the  full  stage- 
frequency  distribution,  for  a  given  return  period,  may  be  used  as  an  estimate 
of  error.  If  so,  one  might  expect  that  the  range  would  be  a  function  of  the 
form  (Leon  Borgman,  personal  communication,  1988): 

R  -  k  nc  (9) 

where:  R  =  the  range, 

n  =  sample  size,  and 
k  and  c  =  function  parameters. 

Figure  58  shows  the  range  plotted  as  a  function  of  sample  size  for  the  24 
events  selected  for  the  Surge  Barrier  study  (COM  1988),  whose  analysis  yielded 
an  error  of  approximately  0.24  feet.  With  the  inclusion  of  nine  additional 
storm  events  from  this  study  (i.e.,  a  total  of  33  storm  events)  reduces  the 
estimated  error  to  0.15  feet.  This  error  value  is  well  within  the  accuracy  of 
numerical  models  used  herein  and  in  similar  flooding  projects. 
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Tal'xe  5 


Selected  Storm  Events 


Incremental 

Lag 

Tide 

Peak  Surge 

Probability 

Cumulative 

Storm1 

/Hours) 

No. 

/Ft.  NGVD) 

of  Occurrence 

Probabilitv 

N017 

4.5 

3 

12.25 

0.0013 

0.0013 

H242* 

11.5 

2 

11.75 

0.0007 

0.0020 

H135 

15.0 

1 

11.25 

0.0006 

0.0026 

H203* 

8.5 

3 

11.00 

0.0006 

0.0032 

H144 

16.0 

5 

10.75 

0.0007 

0.0039 

H196* 

14.5 

1 

10.5 

0.0008 

0.0047 

H117 

14.5 

3 

10.25 

0.0009 

0.0056 

N014* 

4.5 

4 

10.0 

0.0010 

0.0066 

H192 

20.5 

4 

9.75 

0.0017 

0.0083 

H200 

10.5 

4 

9.50 

0.0021 

0.0104 

H090 

18.5 

4 

9.25 

0.0024 

0.0128 

N01S 

10.5 

4 

9.08 

0.0027 

0.0155 

N001 

4.0 

3 

8.92 

0.0027 

0.0182 

H184 

20.5 

5 

8.75 

0.0041 

0.0223 

H225 

2.0 

2 

8.58 

0.0053 

0.0276 

H151 

13.0 

1 

8.42 

0.0053 

0.0329 

N020 

3.5 

3 

8.25 

0.0075 

0.0404 

H027 

13.0 

2 

8.08 

0.0093 

0.0497 

H191 

10.5 

5 

7.92 

0.0093 

0.0590 

H23S 

15.0 

5 

7.75 

0.0120 

0.0710 

H025 

5.0 

5 

7.58 

0.0140 

0.0850 

H041 

10.5 

5 

7.42 

0.0140 

0.0990 

H250 

7.0 

3 

7.25 

0.0210 

0.1200 

N041* 

17.5 

1 

7.125 

0.0200 

0.1400 

H078 

21.5 

2 

7.00 

0.0300 

0.1700 

H008* 

22.5 

4 

6.875 

0.0200 

0.1900 

N028 

2.5 

3 

6.75 

0.0600 

0.2500 

H066* 

21.5 

3 

6.25 

0.0700 

0.3200 

N009 

6.5 

1 

6.25 

0.1400 

0.4600 

H166* 

8.5 

2 

6 

0.1300 

0.5900 

H249 

1.5 

4 

5.75 

0.2700 

0.8600 

N018* 

16.0 

2 

5.5 

0.1310 

0.9910 

H001 

1.5 

5 

5.25 

0.0090 

1 . 0000 

Note: 

Hxxx  = 

refers  to 

a  hurricane 

Nxxx  = 

refers  to 

a  northeaster 

refer b  to  new 

storm  event  for 

this  study 
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<D  >  CO  +*  —  O  C 


Return  Period  for  Peak  Surge 

Node  172 


10  50  100  1000 


Return  Period  (years) 

#  Storms  in  Sample 

-  Return  Period  -£*  24  Storms 

12  Storms  — 8  Storms 

Figure  56.  Return  period  for  peak  surge  at  Node  172. 
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•  Selected  events  (CDM  1988) 
O  Additional  events  for  this  study 


GA9N  133d  NOI1VA313 


Stage  frequency  curve  at  entrance  to  Newark  Bay 
(source:  CDM  1988) 


£E  CO  C  0)0) 


Figure  58.  Range  of  elevation  vs  size  of  sample  for  return  periods  of  10  to  100  years 


PART  VI:  SEA  LEVEL  RISE 


87.  Marine  water  level  records  in  the  vicinity  of  New  York  indicate  a 
long-term  rise  in  mean  sea  level  since  at  least  the  1920s.  Because  sea  level 
influences  surge  elevations  in  the  modeled  system,  sea  level  rise  must  be 
taken  into  account  in  order  to  provide  accurate  simulations  of  surge 
elevations  under  future  conditions. 

88.  Figure  59  shows  mean  annual  and  19-year  averages  of  sea  level  at  The 
Battery.  Mean  sea  level  fluctuates  substantially  from  year  to  year,  but  the 
long-term  trend  becomes  apparent  when  a  19-year  moving  average  is  applied  to 
the  annual  means.  Nineteen  years  is  approximately  the  period  of  one  of  the 
lowest- frequency  tidal  components,  so  that  averaging  over  19  years  smooths  all 
signals  with  periods  less  than  19  years.  According  to  calculations  made  by 
the  National  Ocean  Survey  (NOAA  1988),  mean  sea  level  rise  at  The  Battery 
since  1920  has  been  2.7  millimeters  (0.00886  feet)  per  year.  This  estimate  is 
based  upon  historical  observations. 

89.  Previous  simulations  (COM  1983b)  assumed  a  value  of  0.5  feet  for 
mean  sea  level  above  NGVD.  This  value  is  approximately  correct  for  the  year 
1950,  when  the  long-term  average  was  about  5.48  feet  (0.53  feet  NGVD),  on  the 
tide  staff.  From  Figure  59,  long-term  mean  sea  level  in  1950  was 
approximately  0.3  feet  higher  than  that  in  1929. 

90.  Simulations  of  surge  elevations  under  future  conditions  in  the  year 
2050  incorporated  an  estimated  sea  level  rise  of  0.89  feet  at  The  Battery. 

This  estimate  assumes  that  sea  level  rise  continues  at  the  constant  historical 
rate  of  0.00886  feet  per  year  over  the  next  100  years  from  1950. 

91.  Figure  59  shows  that  there  can  be  large  interannual  variability  in 
annual  mean  sea  level.  One  large  fluctuation  was  in  the  early  1980s  when  sea 
level  appeared  to  increase  by  approximately  0.3  feet  in  one  year.  Such 
irregular  changes  cannot  be  anticipated  in  simulations  of  surge  elevations 
under  future  conditions. 
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Figure  59.  Mean  sea  level  1921-1990  at  the  Battery,  New  York 


PART  VZXs  SENSITIVITY  ANALYSIS 
Events  Selected 


92.  The  purpose  of  the  sensitivity  analysis  was  to  complement  the 
coincidence  analysis  by  determining  the  order  of  magnitude  of  coincident  tidal 
surge  and  river  flow  events,  both  with  and  without  the  proposed  improvements. 
In  this  analysis,  we  postulated  that  the  effect  of  large  river  flows  in 
increasing  tidal  elevation  near  Kearny  Point  will  be  small.  This  postulate 
was  based  on  the  previous  analysis  of  CDM  (1983b),  in  which  storm  surges  and 
river  flows  were  assumed  to  be  independent  events.  The  superimposed  stage- 
frequency  distributions  in  that  study  showed  that  the  distribution  produced  by 
river  flows  alone  was  small  for  locations  well  downstream  in  the  Passaic  and 
Hackensack  Rivers,  and  in  Newark  Bay.  Based  on  this  analysis,  we  modeled  the 
river  flows  in  the  majority  of  the  sensitivity  simulations  as  being  constant 
because,  if  the  effects  are  small  in  the  region  of  interest,  then  the  effects 
produced  by  an  actual  hydrograph  will  be  smaller  still. 

93.  From  the  available  storm  and  non-storm  events,  we  selected  three 
that  cover  the  range  of  surge  elevations  expected  from  the  complete  set  of 
storms.  The  three  events  selected  were: 

1.  The  tide  of  11-13  October  1976  -  this  represents  a  non-surge 
condition. 

2.  Synthetic  hurricane  H249  -  this  represents  a  small  hurricane/tide 
combination  that  produced  a  small  maximum  tidal  elevation  of  5.75 
feet  at  the  mouth  of  Newark  Bay. 

3.  Synthetic  hurricane  H135  -  this  represents  a  large  hurricane/ tide 
combination  that  produced  a  large  maximum  tidal  elevation  of  11.25 
feet  at  the  mouth  of  Newark  Bay. 

94.  Each  of  these  three  events  was  simulated,  with  various  coincident 
river  flows,  for  year  2050  conditions.  Consequently,  a  datum  plane  shift  of 
0.89  feet  was  used  to  include  the  effect  of  sea  level  rise  from  1950-2050. 


liver  Flow  Conditions 


95.  The  New  York  Corps  has  analyzed  the  hydrology  of  the  Passaic  River 
Basin  above  the  Dundee  Dam,  and  has  developed  flow  hydrographs  for  the  100- 
year  river  flood  event.  Figure  60  shows  the  Phase  1  General  Design  Memorandum 


105 


100-year  flow  hydrographs  for  the  year  2050,  with  and  without  the  effect  of 
the  diversion  tunnel.  Without  the  diversion  tunnel,  the  100-year  hydrograph 
is  double  peaked,  with  a  maximum  flow  of  33,645  cfs.  The  effect  of  the 
diversion  tunnel  is  to  decrease  the  time  of  travel  of  the  flood  wave,  and 
gives  a  maximum  flow  of  48,544  cfs,  with  approximately  30,000  cfs  through  the 
diversion  tunnel  and  20,000  cfs  in  the  Passaic  River. 

96.  In  the  sensitivity  analysis,  the  three  events  were  simulated  for  the 
following  four  conditions: 

1.  Future  unimproved  conditions  with  low  river  flows  -  for  the  tide 
of  11-13  October  1976,  the  observed  river  flows  were  used;  for  the 
two  synthetic  hurricanes,  the  one-year  flow  of  5,285  cfs  in  the 
Passaic  River  and  1,630  cfs  in  the  Hackensack  River  were  used. 

2.  Future  unimproved  conditions  with  high  river  flows  -  the  100-year 
flows  of  33,645  cfs  in  the  Passaic  River,  and  7,135  cfs  in  the 
Hackensack  River,  were  used. 

3.  Future  improved  conditions  with  no  flood  levees  -  the  100-year 
flows  of  30,000  cfs  through  the  diversion  tunnel,  20,000  cfs  in 
the  upper  Passaic  River,  and  7,135  cfs  in  the  Hackensack  River 
were  used. 

4.  Future  improved  conditions  with  flood  levees  -  the  same  as 
condition  3,  except  that  overbank  links  in  the  vicinity  of 
proposed  flood  protection  levees  were  removed,  to  simulate  the 
elimination  of  overbank  flooding  due  to  the  construction  of  levees 
of  sufficient  height. 

Analysis  of  Simulations 

97.  The  four  conditions  were  simulated  for  each  of  the  three  events,  and 
the  results  plotted  at  various  locations  along  the  Passaic  and  Hackensack 
Rivers,  and  in  Newark  Bay.  The  results  are  plotted  in  Figures  61  through  84 
for  several  sites  selected  in  consultation  with  the  PRD.  The  peak  tidal  surge 
elevations  at  these  sites  are  compared  in  Tables  6  through  8,  for  the  three 
different  events. 

98.  The  results  show  that  for  those  stations  in  the  lower  Passaic  and 
Hackensack  Rivers,  and  in  Newark  Bay,  the  magnitude  of  the  river  flow  (low 
versus  100-year  event),  and  the  distribution  of  flows  (future  unimproved 
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versus  future  improved)  have  little  effect  on  the  resulting  tidal  surge 
elevations.  Prom  Tables  6  through  8,  it  can  be  seen  that  the  difference  is 
only  a  few  tenths  of  a  foot,  and  we  conclude  that  inflow  distributions  and 
magnitudes  and  flood  protection  levees  have  a  negligible  effect  on  tidal 
elevations  in  the  lower  part  of  the  system.  Certainly  the  variability  is 
within  the  accuracy  of  the  model,  and  within  the  factor  of  safety  normally 
designed  into  this  type  of  project. 

99.  For  those  stations  farther  up  the  Passaic  River,  the  effect  of  the 
diversion  tunnel  is  to  reduce  the  magnitude  of  the  100-year  flood  by  several 
feet.  Thus  in  the  tidal  portion  of  the  system,  below  the  river  dams,  the 
overall  effect  of  the  tunnel  diversion,  with  or  without  flood  protection 
levees,  is  to  reduce  the  river  flood  elevations  in  the  upper  portion  of  the 
tidal  Passaic  River,  while  not  significantly  increasing  flood  elevations  in 
the  remaining  tidal  system.  The  benefit  of  the  diversion  tunnel  may  by  itself 
warrant  considering  a  diversion  tunnel  for  the  project  area. 
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Table  6 


Maximum  Predicted  Water  Elevations  (feet)  at  Selected  Nodes 
October  1976  Tides  and  Four  Project  Scenarios 


Node 


Year  2050  Conditions 


Low  Flows 
1976 

Unimproved 


100-Year 

Flows 

Unimproved 


100-Year 
Flows  With 
Tunnel 


100-Year 
Flows  With 
Tunnel 
and  Levees 


Newark  Bav 

49  (Newark  Airport) 
262  (Kearny  Point) 

Passaic  River 

173  (East  Newark) 
217  (Third  River) 

226  (Garfield) 

Hackensack  River 

182  (Snake  Hill) 

190  (Berry's  Creek) 
213  (New  Milford) 


3.75 

3.78 


3.86 

3.95 

4.06 


3.82 

3.82 

4.03 


3.80 

3.83 


7.06 

12.60 

18.61 


3.92 

4.02 

10.20 


3.81 

3.89 


4.99 

8.05 

13.93 


3.95 

4.04 

10.21 


3.81 

3.89 


4.99 

8.05 

13.93 


3.95 

4.04 

10.21 


108 


Table  7 


Maximum  Predicted  Hater  Elevations  (feet)  at  Selected  Nodes  for 
Synthetic  Hurricane  135  and  Four  Project  Scenarios 


Node  _ Year  2050  Conditions 


Low  Flows 

100-Year 

100-Year 

100-Year 
Flows  With 

1976 

Flows 

Flows  With 

Tunnel 

UnimDroved 

UnimDroved 

Tunnel 

and  Levees 

Newark  Bav 


49  (Newark  Airport) 

13.14 

13.41 

13.44 

13.82 

262  (Kearny  Point) 

13.41 

13.86 

13.89 

14.09 

Passaic  River 

173  (East  Newark) 

13.58 

14.86 

14.89 

14.99 

217  (Third  River) 

13.97 

16.34 

15.80 

15.99 

226  (Garfield) 

14.89 

19.75 

17.52 

17.69 

Hackensack  River 

162  (Belleville) 

12.87 

13.23 

13.24 

13.29 

190  (Berry's  Creek) 

11.36 

11.72 

11.74 

11.78 

213  (New  Milford) 

8.43 

11.63 

11.66 

11.65 
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Table  8 


Maximum  Predicted  Water  Elevations  (feet)  at  Selected  Nodes  for 
Synthetic  Hurricane  249  and  Four  Project  Scenarios 


Node  _ Year  2050  Conditions 


Low  Flows 

100-Year 

100-Year 

100-Year 
Flows  With 

1976 

Flows 

Flows  With 

Tunnel 

UnimDroved 

UnimDroved 

Tunnel 

and  Levees 

Newark  Bav 


49  ( Newark  Airport ) 

5.65 

5.70 

5.73 

5.73 

262  (Kearny  Point) 

5.69 

5.75 

5.80 

5.80 

Passaic  River 

173  (East  Newark) 

5.90 

8.46 

6.67 

6.67 

217  (Third  River) 

6.19 

13.14 

9.17 

9.17 

226  (Garfield) 

7.08 

18.74 

14.39 

14.39 

Hackensack  River 

182  (Belleville) 

5.66 

5.75 

5.78 

5.78 

190  (Berry's  Creek) 

5.54 

5.68 

5.71 

5.71 

213  (New  Milford) 

5.89 

10.62 

10.63 

10.63 
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Figure  60.  Passaic  River  Hydrographs 


Figure  61.  Predicted  water  level 
elevations  at  Newark  Airport  for  year 
2050  conditions  and  October  1976  tides 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

itt  nw  FLOW 


Figure  62.  Predicted  water  level 
elevations  at  Kearny  Point  for  year 
2050  conditions  and  October  1976  tides 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  63.  Predicted  water  level 
elevations  at  East  Newark  for  year 
2050  conditions  and  October  1976  tides 


WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FLOW 


Figure  64.  Predicted  water  level 
elevations  at  Third  River  for  year 
2050  conditions  and  October  1976  tides 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  65.  Predicted  water  level 
elevations  at  Garfield  for  year  2050 
conditions  and  October  1976  tides 
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Figure  66.  Predicted  water  level 
elevations  at  Snake  Hill  for  year  2050 
conditions  and  October  1976  tides 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FL0W 


Figure  67.  Predicted  water  level 
elevations  at  Berry’s  Creek  for  year 
2050  conditions  and  October  1976  tides 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FLOW 


Figure  68.  Predicted  water  level 
elevations  at  New  Milford  for  year 
2050  conditions  and  October  1976  tides 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  69.  Predicted  water  level 
elevations  at  Newark  Airport  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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NO  TUNNEL,  NO  LEVEES,  LOW  FLOW -  NO  TUNNEL,  NO  LEVEES ,  HIGH  FLC 

WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FL0W 


Figure  70.  Predicted  water  level 
elevations  at  Kearny  Point  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  71.  Predicted  water  level 
elevations  at  East  Newark  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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NO  TUNNEL,  NO  LEVEES,  LOW  FLOW -  NO  TUNNEL,  NO  LEVEES,  HIGH  FLOW 

WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  72.  Predicted  water  level 
elevations  at  Third  River  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FLOW 


Figure  73.  Predicted  water  level 
elevations  at  Garfield  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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NO  TUNNEL,  NO  LEVEES,  LOW  FLOW -  NO  TUNNEL,  NO  LEVEES ,  HIGH  FLOW 

WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

ri  n  w  FLOW 


Figure  74.  Predicted  water  level 
elevations  at  Snake  Hill  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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WITH  TUNNEL.  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FL0W 


Figure  75.  Predicted  water  level 
elevations  at  Berry’s  Creek  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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NO  TUNNEL,  NO  LEVEES,  LOW  FLOW -  NO  TUNNEL,  NO  LEVEES ,  HIGH  FL( 

WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FL0W 


Figure  76.  Predicted  water  level 
elevations  at  New  Milford  for  year 
2050  conditions  and  synthetic 
hurricane  135 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FLOW 


Figure  77.  Predicted  water  level 
elevations  at  Newark  Airport  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FLOW 


Figure  78.  Predicted  water  level 
elevations  at  Kearny  Point  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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WITH  TUNNEL.  NO  LEVEES,  HIGH  .  WITH  TUNNEL.  WITH  LEVEES,  HIGH 

irr  nw  FLOW 


Figure  79.  Predicted  water  level 
elevations  at  East  Newark  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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Figure  80.  Predicted  water  level 
elevations  at  Third  River  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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WITH  TUNNEL,  NO  LEVEES.  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  81.  Predicted  water  1 
elevations  at  Garfield  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FLOW 


Figure  82.  Predicted  water  level 
elevations  at  Snake  Hill  for  year 
2050  conditions  and  synthetic 
hurricane  249 


133 


WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 


Figure  83.  Predicted  water  level 
elevations  at  Berry's  Creek  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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NO  TUNNEL,  NO  LEVEES  ,  LOW  FLOW -  NO  TUNNEL,  NO  LEVEES  ,  HIGH  FLOW 

WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 

FLOW  FL0W 


PART  VZIX:  STAGE -FREQUENCY  DISTRIBUTIONS  FOR  FUTURE  CONDITIONS 


100.  The  results  of  the  sensitivity  analysis  demonstrate  that  the 
magnitude  and  distribution  of  river /diversion  tunnel  flows  have  little  effect 
(less  than  0.1  ft)  on  tidal  elevations  in  the  lower  parts  of  the  Passaic  and 
Hackensack  Rivers,  and  in  Newark  Bay.  The  effect  of  including  flood 
protection  levees  in  the  vicinity  of  Kearny  Point  is  to  increase  storm  tide 
heights  by  less  than  0.2  foot  over  future  unimproved  conditions  in  this 
vicinity  and  by  less  than  0.4  foot  adjacent  to  Newark  Airport.  Consequently, 
once  stage-frequency  distributions  have  been  developed  for  the  future 
unimproved  condition,  the  effects  of  future  improvements,  such  as  the  tunnel 
and  levee  construction,  will  be  small.  The  effect  of  the  improvements  on 
tidal  elevations  is  within  the  accuracy  of  the  development  of  the  stage- 
frequency  distributions,  and  it  is  standard  practice  to  include  an  appropriate 
factor  of  safety  in  any  design. 

101.  The  33  synthetic  storms  selected  for  a  "modified  bootstrap"  of  the 
underlying  joint  probability  analysis,  were  simulated  with  a  datum  shift  of 
0.89  feet,  which  represents  the  estimated  sea  level  rise  from  1950-2050.  Two 
scenarios  were  simulated:  (1)  future  unimproved,  with  no  levees  or  tunnel 
discharge,  and  one-year  river  flows  in  Passaic  (5,285  cfs)  and  Hackensack 
(1,630  cfs)  rivers;  and  (2)  future  improved,  with  levees,  one-year  river  flows 
in  the  Passaic  River  (1,000  cfs  in  the  river  channel  and  6,000  cfs  in  the 
tunnel,  and  one-year  flow  of  1,630  cfs  in  the  Hackensack  River.  The  resulting 
distributions  at  various  locations  in  the  system  are  shown  on  Figures  85 
through  92. 

102.  It  is  important  to  recognize  that  the  resulting  stage-frequency 
distributions  in  these  figures  do  not  include  the  joint  probabilities  of 
significant  river  flows.  As  the  effect  of  river  flows  in  the  lower  part  of 
the  system  is  small,  the  stage- frequency  distributions  in  these  areas  will  be 
accurate.  However,  in  the  upper  tidal  portions  of  the  Passaic  and  Hackensack 
Rivers,  where  the  effects  of  river  flows  are  significant,  the  stage-frequency 
curves  represent  purely  coastal  processes.  For  the  primary  aim  of  this  study, 
we  are  focussing  on  the  lower  system,  namely  around  Kearny  Point,  and  have 
shown  clearly  that  upstream  flooding  will  not  be  significantly  worsened  by  the 
discharges  from  the  proposed  diversion  tunnel. 
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Figure  84.  Predicted  water  level 
elevations  at  New  Milford  for  year 
2050  conditions  and  synthetic 
hurricane  249 
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NO  TUNNEL,  NO  LEVEES  ,  LOW  FLOW -  NO  TUNNEL,  NO  LEVEES  ,  HIGH  FLOW 

WITH  TUNNEL,  NO  LEVEES,  HIGH  .  WITH  TUNNEL,  WITH  LEVEES,  HIGH 
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1-YR  FLOWS  -  1-YR  FLOWS  WITH 

UNIMPROVED  TUNNEL  AND  LEVEES 


age-frequency  curves  for 
rny  Point)  in  Newark  Bay 
synthetic  storms 
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PROBABILITY  OF  EXCEEDANCE 

1-YR  FLOWS  -  1-YR  FLOWS  WITH 

UNIMPROVED  TUNNEL  AND  LEVEES 


Figure  87.  Stage-frequency  curves  for 
node  173  (East  Newark)  on  the  Passaic 
River  for  33  synthetic  storms 
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1-YR  FLOWS  -  1-YR  FLOWS  WITH 

UNIMPROVED  TUNNEL  AND  LEVEES 
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1-YR  FLOWS  -  1-YR  FLOWS  WITH 

UNIMPROVED  TUNNEL  AND  LEVEES 


Figure  90.  Stage-frequency  curves  for 
ode  182  (Snake  Hill)  on  the  Hackensack 
River  for  33  synthetic  storms 


Figure  91.  Stage-frequency  curves  for 
node  190  (Berry's  Creek)  on  the 
Hackensack  River  for  33  synthetic  storms 
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1-YR  FLOWS  -  1-YR  FLOWS  WITH 

UNIMPROVED  TUNNEL  AND  LEVEES 


Figure  92.  Stage-frequency  curves  for 
node  213  (New  Milford)  on  the  Hackensack 
River  for  33  synthetic  storms 
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PART  ZZt  SUMMARY  AND  CONCLUSIONS 


103.  The  New  York  Corps  is  proposing  to  build  a  20-mile  long,  deep  rock 
tunnel  about  40  feet  in  diameter,  to  divert  flood  waters  from  the  upper 
Passaic  River  Basin  to  Newark  Bay.  One  objective  of  this  study  was  to  examine 
the  effect  that  discharges  from  the  diversion  tunnel  might  have  on  tidal 
elevations  in  Newark  Bay,  and  the  lower  tidal  parts  of  the  Passaic  and 
Hackensack  Rivers.  Another  objective  was  to  develop  stage-frequency  curves  at 
strategic  locations  for  the  design  of  storm  tide  levees.  A  particular  concern 
was  whether  tunnel  flows  might  be  coincident  with  storm  tide  surges  from  New 
York  Bight,  and  possibly  lead  to  increased  flooding  in  the  downstream  areas, 
such  as  the  Meadowlands  and  Newark  Airport. 

104.  To  examine  this  possibility,  a  version  of  the  Dynamic  Estuary  Model 
(DEM),  previously  applied  to  the  system  (CDM  1983a,  1983b,  1988),  was 
revalidated  fcr  this  study  using  recent  hydrographic  data.  The  revalidation, 
having  increased  resolution  built  into  the  model  and  incorporating  nx>re  recent 
data,  showed  better  agreement  than  previous  studies. 

105.  A  coincidence  analysis  was  performed  to  examine  the  distributions 
of  river  flows  and  storm  surges  (the  predicted  tidal  effect  was  removed) . 

This  analysis  showed  that  there  is  a  weak  relationship  between  river  flows  and 
tidal  elevations.  Coincidences  of  large  surges  and  large  river  flows  were 
found  to  be  extremely  rare  occurrences. 

106.  The  statistical  procedure  used  in  the  Hackensack  River  Surge 
Barrier  Study  (CDM  1988)  was  based  on  a  "modified  bootstrap"  method.  This 
method  was  applied  to  the  underlying  joint  probability  method  analysis  (CDM 
1983a,  1983b).  From  the  more  than  70,000  surge/tide/lag  combinations,  24 
events  were  selected  to  provide  a  reasonable  approximation  of  the  stage- 
frequency  distribution.  For  this  study,  we  included  an  additional  nine 
storms,  for  a  total  of  33  storm  events.  These  events  were  selected  to  fill 
out  the  frequency  distribution  for  return  periods  of  less  than  10  years  and 
greater  than  100  years. 

107.  As  the  future  project  conditions  were  developed  for  the  year  2050, 
we  examined  the  difference  between  mean  sea  level  (MSL)  and  NGVD,  and 
concluded  that  the  previous  studies  (CDM  1983b,  1988)  were  based  on  a 
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correction  appropriate  for  about  1950.  Using  an  estimate  of  2.7  mn/yr  (0.0089 
ft/yr)  for  sea  level  rise  (NOAA  1988),  we  estimated  a  sea  level  rise  of  0.89 
feet  between  1950  and  2050.  This  value  was  used  to  adjust  the  specified  datum 
for  future  condition  simulations  in  this  study. 

108.  To  reinforce  the  conclusions  of  the  coincidence  analysis,  a 
sensitivity  analysis  was  performed  in  which  three  tidal  events  ((1)  a  typical 
tide,  (2)  a  small  hurricane,  and  (3)  a  large  hurricane)  were  simulated  for 
four  future  conditions  ((1)  unimproved  with  a  small  river  flow,  (2)  unimproved 
with  a  100-year  river  flow,  (3)  improved  with  100-year  river  and  tunnel  flowB, 
and  (4)  improved  with  flood  protection  levees  and  100-year  river  and  tunnel 
flows).  The  results  showed  that  the  magnitude  and  distribution  of  river  and 
tunnel  flows  had  a  very  small  effect  (less  than  0.2  ft)  on  tidal  elevations  in 
the  lower  parts  of  the  Passaic  and  Hackensack  Rivers,  and  in  Newark  Bay. 

109.  Having  concluded  that  the  magnitude  and  distribution  of  river  and 
tunnel  flows  had  little  effect  on  tidal  elevations  in  the  lower  system,  stage- 
frequency  distributions  for  future  unimproved  conditions  were  developed  at 
various  locations  using  the  33  selected  synthetic  storm  events  and  the  one- 
year  river  flows  in  the  Passaic  and  Hackensack  Rivers.  The  effect  on  these 
stage-frequency  distributions  due  to  large  river  and/or  tunnel  flows  will  be 
small,  generally  on  the  order  of  a  few  tenths  of  a  foot,  which  is  within  the 
accuracy  of  the  model  used  in  this  study. 

Conclusions 

110.  The  revalidated  DEM  model  more  accurately  simulated  tidal 
elevations  in  the  Passaic/Hackensack/Newark  Bay  system  than  previous  studies. 
This  was  due  to  the  updating  and  inclusion  of  more  recent  system  geometry 
data.  The  model  is  therefore  a  valuable  tool  for  simulating  tidal  and  river 
flow  scenarios  in  the  system. 

111.  The  coincidence  analysis  confirmed  the  findings  of  the  Hackensack 
River  Surge  Barrier  Study  (COM  1988)  that  large  river  flows  are  unlikely  to 
coincide  with  large  storm  surges  from  New  York  Bight.  Time-history  plots  of 
several  individual  events  suggest  that  there  is  a  relationship  between  storm 
tides  and  river  flows,  with  the  river  flows  lagged  by  several  days,  but  they 
also  serve  to  further  demonstrate  that  the  relationship  is  weak.  In  general, 
for  a  relative  small  surge,  a  wide  range  of  river  flows  is  possible.  For  a 
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relatively  small  (less  than  one-year)  river  flow,  a  wide  range  of  storm  surges 
is  possible.  However,  large  storm  tides  with  large  Passaic  River  flows 
(lagged  by  0  to  2  days)  are  extremely  rare  events. 

112.  The  data  were  further  analyzed  to  determine  the  average  surge 
condition  that  might  occur  within  72  hours  of  specific  significant  river 
flows.  The  results  suggested  that  over  the  range  of  river  flows  observed,  an 
average  storm  surge  is  on  the  order  of  1.5  to  2  feet.  Therefore,  we  conclude 
that,  if  a  Passaic  River  flow  model  is  to  be  developed  to  simulate  significant 
flow  events,  the  downstream  boundary  condition  could  be  the  superposition  of 
an  average  tide  plus  a  storm  surge  of  1.5  to  2  feet. 

113.  In  a  similar  fashion,  the  data  were  also  analyzed  to  determine  the 
average  Passaic  River  flow  that  might  occur  within  72  hours  of  storm  surges  of 
different  elevations.  The  results  showed  that  for  all  storm  surge  elevations 
observed,  an  average  Passaic  River  flow  of  1000-2000  cfs  would  be  found.  From 
this  finding,  we  recommended  (and  used  for  this  study)  that  the  one-year  river 
flows  in  the  major  tributaries  of  the  system  should  be  used  for  storm  tide 
simulations.  In  the  Passaic  River,  the  one-year  flow  is  in  excess  of  5000 
cfs,  or  about  three  times  the  average  flow  coincident  with  a  storm  surge. 

114.  NOS  concludes  that  sea  level  is  rising  at  an  approximate  rate  of 
2.7  mm/yr  ( 0.0089  ft/yr) .  Based  on  the  previous  model  studies,  we  estimated  a 
sea  level  rise  of  0.89  feet  in  the  100  years  from  1950  to  2050. 

115.  The  sensitivity  analysis  performed  for  future  conditions  (including 
sea  level  rise)  for  three  tide  events  and  four  river  flow  conditions, 
demonstrated  that  the  distribution  and  magnitude  of  river  and  tunnel 
discharges  has  a  very  small  effect  (less  than  0.1  ft)  on  tidal  elevations  in 
the  lower  part  of  the  system,  including  the  lower  tidal  parts  of  the  Passaic 
and  Hackensack  Rivers  and  Newark  Bay.  The  placement  of  flood  protection 
levees  around  Kearny  Point  also  has  a  small  effect  (less  than  0.2  ft)  on  tidal 
elevations.  Upstream  in  the  tidal  Passaic  River,  the  effect  of  the  Diversion 
Tunnel  is  to  lower  water  elevations  during  significant  river  flow  events.  For 
example,  even  though  the  peak  fresh  water  discharge  to  the  lower  Passaic  may 
increase  (from  33,645  cfs  to  48,544  cfs)  for  the  100-year  flow,  flood 
elevations  will  be  lowered  by  several  feet  as  a  result  of  the  tunnel  bypassing 
the  upper  part  of  the  tidal  river.  The  flood  characteristics  of  the 
Hackensack  remain  essentially  the  same. 
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116.  Finally,  stage- frequency  curves  for  various  location  in  the  system 
were  developed  using  33  synthetic  storm  events  and  the  one-year  river  flows  in 
the  Passaic  and  Hackensack.  The  distributions  are  similar  to  the 
distributions  developed  previously  (CDM  1983a,  1988)  with  an  increase  in 
overall  elevations  due  to  sea  level  rise.  These  distributions  permit  the 
selection  of  appropriate  flood  protection  levees  around  Kearny  Point  to 
protect  against  a  given  return  period  storm  surge  event. 

Summary  Conclusions 

117.  From  the  analyses  in  this  study,  the  following  conclusions  are 
drawn: 

1.  The  coincidence  of  large  storm  tides  and  large  river  flows  is  an 
extremely  unlikely  combination,  and  can  be  treated  as  essentially 
independent  events.  However,  the  likelihood  that  somewhat 
elevated  river  flows  occur  following  storm  tides  (by  0-2  days)  is 
quite  strong. 

2.  The  tidal  boundary  condition  for  a  river  flow  model  of  the  Passaic 
River  system  could  be  the  superposition  of  an  average  tide  and  a 
storm  surge  of  1.5  to  2  feet. 

3.  From  the  coincidence  analysis,  we  concluded  that  it  is  reasonable 
to  assume,  as  a  conservative  condition,  that  the  one-year  river 
flows  in  the  Passaic  and  Hackensack,  can  occur  with  significant 
surge  events.  This  conclusion  is  consistent  with  CDM  (1988). 

4.  The  coincidence  analysis  was  used  to  develop  a  storm  tide  stage- 
frequency  distribution  of  events  in  which  the  Passaic  River  flow 
exceeded  6000  cfs.  This  analysis  can  be  used  to  determine 
appropriate  tidal  heights  for  design  of  tunnel  outlet  works. 

5.  Sea  level  rise  is  taken  to  be  2.7  mm/yr  (0.0089  ft/yr),  or  0.89 
feet  between  1950  and  2050. 

6.  The  magnitude  and  distribution  of  significant  river  and  tunnel 
discharges  have  a  very  small  effect  (less  than  0.2  ft)  on  tidal 
elevation  in  the  lower  Passaic  River,  Hackensack  River,  and  Newark 
Bay. 

7.  The  effect  of  tunnel  discharges  to  Newark  Bay  will  have  little 
effect  (less  than  0.1  ft)  on  the  existing  flood  characteristics  of 
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the  Hackensack  River,  including  the  flood-prone  areas  in  the 
vicinity  of  the  Meadow lands. 

8.  The  effect  of  tunnel  discharges  to  Newark  Bay  will  reduce  river 
flood  elevations  in  the  upper  tidal  Passaic  River  by  several  feet. 

9.  The  stage-frequency  distributions  developed  for  the  vicinity  of 
Kearny  Point  will  assist  in  the  selection  of  flood  protection 
levee  heights.  The  levees  will  be  designed  to  protect  against 
tidal  surges  with  a  given  return  period.  The  distributions  in  the 
upper  tidal  Passaic  and  Hackensack  Rivers  do  not  include  the 
effect  of  significant  river  flows,  and  therefore  reflect  only 
tidal  conditions. 

10.  The  inclusion  of  future  project  conditions,  including  flood 
control  levees  in  the  vicinity  of  Kearny  Point,  changes  the  future 
unimproved  stage-frequency  distributions  by  less  than  0.2  feet  in 
the  lower  part  of  the  Passaic  and  Hackensack  Rivers,  and  in  Newark 
Bay. 

11.  Dredging  operation  in  the  Kills  and  Newark  Bay  will  have  little 
effect  on  circulation  and  water  levels  in  the  study  area.  Water 
levels  may  change  by  less  than  0.15  feet,  and  the  tidal  phase  may 
be  altered  by  less  than  1/2  hour. 
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APPENDIX  A:  USERS  GUIDE  FOR  NEWARK  BAY  VERSION  OP 
THE  DYNAMIC  ESTUARY  MODEL  —  DEMNEWARK 

P«ar  Considerations 

Al.  An  early  version  of  the  Dynamic  Estuary  Model  (DEM)  was  developed 
by  Water  Resources  Engineers  (1977).  It  was  the  first  version  to  incorporate 
wetting  and  drying.  The  code  was  substantially  modified  for  a  hurricane  surge 
study  of  the  New  York  Bight  area  (COM  1983a)  and  was  used  in  this  form  for  the 
Passaic  River  surge  study  (CDM  1983b).  For  the  Hackensack  River  surge  barrier 
study  (CDM  1988),  the  code  was  further  modified  to  include  culverts  and  broad- 
crested  weirs  (to  simulate  barrier  overtopping).  This  version  was  called 
DEMHACK.  For  the  Passaic  River  Flood  Control  Study,  described  in  this  report, 
the  Hackensack  River  study  version  was  again  modified  somewhat,  to  be 
compatible  with  Ebasco's  graphics  system.  However,  the  input  formats  are 
unchanged.  This  version  is  called  DEMNEWARK. 

A2.  The  program  is  written  in  FORTRAN  and  compiled  using  a  FORTRAN-77 
compiler.  However,  the  code  uses  no  conventions  for  specific  computers  and 
should  run  on  any  computer  using  FORTRAN-77  with  minor  changes. 

A3.  The  memory  required  by  the  program  should  be  tailored  to  the 
specific  study  being  performed  using  the  PARAMETER  statements  located  in  the 
main  program  and  each  subroutine.  The  maximum  number  of  nodes  or  junctions 
(MXJ)  and  links  or  channels  (MXC)  and  the  maximum  number  of  channels  entering 
a  junction  (MCH)  are  set  in  these  statements  and  used  to  dimension  the  various 
arrays . 

A4.  The  memory  use  and  run  time  requirements  vary  with  the  size  of  the 
network  being  simulated  and  the  computer  used.  To  provide  guidance  for  the 
user,  the  small-scale  grid  of  the  Newark  Bay  area  has  176  nodes  and  263  links. 
Using  a  time  step  of  60  seconds,  one  hour  of  simulated  time  required  41 
seconds  on  an  AST  386/33  with  a  math  coprocessor  (or  about  0.7  seconds  per 
time  step) .  The  source  and  executable  codes  require  approximately  150  KB  and 
200  KB  of  memory  each. 
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Data  Input  Description 


A5.  The  program  inputs  are  supplied  on  14  card  groups,  A  through  N,  as 
summarized  on  Table  Al.  The  actual  input  descriptions  are  given  in  Tables  A2 
and  A3.  In  this  table,  the  variable  types  A,  F  and  I  are  used  for 
alphanumeric,  floating  point  real,  and  integer,  respectively. 

A6.  There  are  three  input  files  and  two  output  files  (see  Figure  Al). 
The  three  input  files  are: 

1.  The  basic  input  information  (as  described  below  except  card  groups 
J  and  K) ; 

2.  The  grid  file  of  node  and  link  information  (card  groups  J  and  K) ; 
and 

3.  Hurricane  or  northeaster  boundary  condition  data  from  the  Tetra 
Tech  model  of  New  York  Bight  (if  XSWCH(3)  *  2  or  3,  respectively]. 

A7.  The  first  two  files  are  described  in  Tables  Al  and  A2.  The  third 
file  (the  hurricane  or  northeaster  file)  is  a  binary  file  created  from  the 
Tetra  Tech  model  of  New  York  Bight  and  is  not  described  here.  Example  input 
files  are  listed  in  Appendix  B. 

Output  Description 

A8.  The  program  creates  two  output  files  (see  Figure  Al).  The  first 
is  a  tabular  output  file  in  which  input  data  is  echoed,  if  the  appropriate 
options  are  specified,  and  velocities/flows  and  heads  are  printed  as  time 
histories  at  each  link  or  node  for  each  day  of  the  simulation.  The  second 
file  is  a  "heads"  output  file  which  contains  heads  at  specified  nodes  every 
half  hour  for  the  simulation  period.  This  file  is  the  basis  of  all  the  time- 
history  plots  in  this  report,  and  is  created  by  editing  the  file  and  reading 
it  into  an  appropriate  graphics  package.  The  output  for  the  input  file  listed 
in  Appendix  B  is  listed  in  Appendix  C. 
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Table  A1 


DEMNEWARK  Input  Card  Groups 


Card  Group 


Brief  Description 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

H 

N 


Alphanumeric  run  description 

Job  control  switches 

Grid  printout  control 

Job  control  information 

Rain  input  cards 

Head  plot  junction  numbers 

Tidal  boundary  cards 

Hind  speed  and  direction  cards 

Map  scale  and  input  datum  reference 

Junction  cards  (nodes) 

Channel  cards  (links) 

Junctions  for  plotting  elevations 

Inflow  hydrographs 

Historical  tidal  boundary  inputs 
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Figure  Al .  Flow  chart  of  DEMNEWARK 
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TABLE  A2 

BASIC  INPUT  DATA  DESCRIPTION 


Card  Group  Card 
Card  Number  Columns 


Variable 
ime  Type 


Variable  Descriptior 


As  ALPHANUMERIC  RUM  DESCRIPTION 


A  1 

1-80 

ALPHA 

A 

A  2 

1-80 

ALPHA 

A 

A  3 

1-80 

TITLE 

A 

A  4 

1-80 

TITLE 

A 

User  Identification  Information 
User  Identification  Information 
Simulation  Title 
Simulation  Title 


Bs  JOB  CONTROL  SNITCHES 


B  1 

1-5 

ISWCH(l) 

I 

1  -  Tidal  system 

2  -  Non-tidal  system  (e.g.,  lake) 

6-10 

ISWCH(2) 

I 

0  -  No  action 

1  -  Suppresses  printout  of  channel 
and  nodal  parameters 

11-15 

ISWCH(3) 

I 

0  -  No  time  series  elevation 

boundary  condition  data  used 
1  -  Read  elevation  boundary 

condition  data  from  this  file 


2  -  Read  hurricane  data  from 
separate  file 


3  -  Read  northeaster  data  from 
separate  file 


16-20 

ISWCH(4) 

I 

0  -  Not  used 

21-25 

ISWCH ( 5 ) 

I 

0  -  Print  channel  velocities  and 
flows 

1  -  Skip  channel  printout 

26-30 

ISWCH(6) 

I 

0  -  Not  used 

31-35 

ISWCH ( 7 ) 

I 

0  -  Not  used 

36-40 

ISWCH (8) 

I 

0  -  No  action 

1  -  Prints  first  wind  and  pressure 
information  and  terminates 
program  (used  for  diagnostic 
checking) 

41-45  ISWCH(9)  I  0  -  Wind  forcing  not  imposed 

1  -  Wind  forcing  included 
46-50  ISWCH(IO)  I  0  -  No  action 

1  -  Surge  elevations  saved  at  1/2 
hour  intervals 

-1  -  Maximum  surge  elevation  and 
time  of  occurrence  saved 
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TABLE  A2 ,  continued 
BASIC  INPUT  DATA  DESCRIPTION 

Card  Group  Card  Variable 

Card  Number  Columns  Name  Type  _ Variable  Description 


C: 

GRID  PRINTOUT 

CONTROL 

C  1 

1-5 

I  CHAR 

I 

0  -  No  action 

1  -  Produces  line  printer 
"map"  of  link-node  grid 
network 

6-10 

XMIN 

F 

Minimum  X  for  grid 

11-15 

XMAX 

F 

Maximum  X  for  grid 

16-20 

YMIN 

F 

Minimum  Y  for  grid 

21-25 

YMAX 

F 

Maximum  Y  for  grid 

D: 

JOB  CONTROL  INFORMATION 

D  1 

1-5 

NDAYS 

I 

Number  of  printout  days  (or 
cycles)  simulated 

6-10 

NPRNT 

I 

Number  of  printout  timesteps 
per  day  (<  26).  Usually  equal 
to  25. 

11-15 

PSTP 

F 

Printout  timestep  (hours). 
Usually  equal  to  one  hour. 

16-20 

DELT 

F 

Hydr au lie  t iroes t ep  ( seconds ) . 

21-25 

TZERO 

F 

Initial  time  of  simulation 
(hours) . 

26-30 

NPLT 

I 

Number  of  stage  plots  (<51) 
(Card  group  F) 

31-35 

EVAP 

F 

Evaporation  rate  (inches  per 
month) . 

36-40 

AVGHEA 

F 

Average  head  used  when  finding 
high  and  low  tides.  Set  to  0. 

46-50 

NPSWRT 

I 

Day  cycle  when  printed  output 
will  start. 

51-55 

NQSWRT 

I 

Number  of  day  cycle  when 
switching  from  curve  fit  to 
measured  tidal  stages. 

56-60 

NJSW 

I 

Number  of  junctions  of 
hydrograph  input  from  cards 
(<21)  (Card  group  M) 

61-65 

INRAIN 

I 

Number  of  points  of  rain 
information  (<51)  (Card  group 
E) 

Number  of  tidal  boundary  nodes 
(<51)  (Card  group  G) 

66-70 

NBOUND 

I 
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TABLE  A2,  continued 
BASIC  INPUT  DATA  DESCRIPTION 


Card  Group 


Card 

Columr 


Variable 
ame  Type 


Variable  Description 


B:  RAIN  INPOT  CARDS 

If  INRAIN  -  0  Skip  Card  Group  E 


INRAIN  pairs  of  values  with  four  per  card 


E  1 

1-10 

RAIN ( 1 ) 

F 

Rate  of  precipitation  (inches 

per  hour) 

11-20 

INTIME (1) 

F 

Time  from  start  of  simulation 
(minutes) 

61-70  RAIN (4)  F  Etc.,  up  to  INRAIN  pairs 

71-80  INTIME (4)  F 


F:  BEAD  PLOT  FUNCTION  NUMBERS 

If  NPLT  =  0  Skip  Card  Group  F 

(note:  If  ISWCH(IO)  *  0,  plots  are  also  sent  to  unit  41) 

Junction  Numbers  to  be  Plotted,  NPLT  values  with  eight  per  card. 

F  1  1-10  JPLT(l)  I  Junction  number  to  be  plotted 


71-80 


JPLT ( 8 ) 


I 


Junction  number  to  be  plotted 


TABLE  A2,  continued 
BASIC  INPUT  DATA  DESCRIPTION 


Card  Group  Card 
Card  Number  Colurr  is 


Variable 
Name  Type 


Variable  Description 


G:  TIDAL  BOUNDARY  CARDS 

If  ISWCH(l)  -2,  Skip  Card  Group  G 

This  card  group  is  repeated  for  each  of  the  NBOUND  tidal  boundary  nodes.  If 
NQSWRT  =  1,  only  the  boundary  junction  number  (JGW)  need  be  specified  using 


only  the 

first  card. 

Otherwise  this 

card  group 

supplies  the  information 

required 

to  perform  a 

curve  fitting 

through  the 

specified  points  to  be 

subsequently  used  as 

the  boundary  condition. 

G  1 

1-5 

KO 

I 

1  -  Only  two  maximum  and  two 
minimum  tide  elevations 
will  be  used  for  curve 
fitting.  NI  must  equal  4. 

0  -  Additional  points  are  read 
in. 

6-10 

NI 

I 

Number  of  tidal  stage  data 
points  (<50). 

11-15 

MAXIT 

I 

Maximum  number  of  iterations 
for  curve  fitting  (<50). 

16-20 

NCHTID 

I 

0  -  Skips  tidal  curve  fitting 
printout . 

1  -  Provides  printout  of  curve 
fitting  calculations 

21-25 

JGW 

I 

Junction  number  of  boundary 

26-35 

PERDT 

F 

Length  of  tidal  day  (hours) 

G  2 

1-10 

TT  ( 1 ) 

F 

Time  (hours)  from  beginning  of 
simulation  of  tidal  stage, 
first  point. 

11-20 

YY(1) 

F 

Tidal  stage  (feet),  first 
point 

• 

TT(NI) 

F 

Time  (hours  from  beginning  of 
simulation  of  tidal  stage, 
last  point) 

• 

YY(NI) 

F 

Tidal  stage  (feet),  last  point 

(use  additional  cards  if 
needed) 
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TABLE  A2,  continued 
BASIC  INPUT  DATA  DESCRIPTION 


Card  Group  Card 
Card  Number  Columns 


Variable 

Name  Type  _ Variable  Description 


Bs  BIND  SPEED  AND  DIRECTION  CARDS 

If  ISWCH(3)  =  0,  Skip  Card  Group  H 


Three  Values  per  Card  with  up  to  50  Cards 


HI 

1-10 

TIMW(l) 

F 

Time  of  wind  change  (hours). 
TIMW  >  99999  terminates  wind 
input . 

11-20 

TWIND(l) 

F 

Wind  magnitude  (mph). 

21-30 

TWDIR(l) 

F 

Wind  direction,  degrees 
measured  counter-clockwise 

from  north  to  direction  in 
which  wind  is  blowing. 


Is  MAP  SCALE  AND  INPOT  DATUM  REFERENCE 


11 

1-10 

I SCALE 

I 

Scale  of  map  uBed  to  measure 
and  y  locations,  and  channel 

lengths .  Program  converts 
input  in  inches,  on  map,  to 
feet.  If  input  is  in  feet, 

I SCALE  *  12. 

21-30 

WRFACT 

F 

Wind  reduction  factor. 

31-40 

DATUM 

F 

Datum  correction  in  feet 
(positive  upward). 

41-50 

FACTN 

F 

Manning's  n  reduction  factor 

L: 

JUNCTIONS  FOR 

PLOTTING  ELEVATIONS 

If  ISWCH(IO) 

=0,  Skip  Card 

Group  L 

LI 

1-5 

NPTS 

I 

Number  of  plotting  locations 

L2 

1-5 

NODES ( 1 ) 

I 

First  plotting  location 

75-80 

NODES (16) 

I 

Sixteenth  plotting  location 

(use  extra  L2  cards  as  needed) 
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TABLE  A2,  continued 
BASIC  INPUT  DATA  DESCRIPTION 

Card  Group  Card  Variable 

Card  Number  Columns  Name  Type  _ Variable  Description 


Mi  INFLOW  HYDROGRAPHS 


If  NJSW  =  0 ,  Skip  Card  Group  M 

Nodes  where  Inflow  will  occur  -16  values  per  card  up  to  the  maximum  of 
NJSW  or  20  nodes. 


Ml 

1-5 

JSW(l) 

I 

Node  for  inflow 

number  1. 

75-80 

JSW(16) 

I 

Node  for  inflow 

number  16. 

M2: 

INFLOW  TIME  HISTORY 

If 

NJSW  ■  0, 

Skip  Card  Group  M2 

M2* 

inflows. 

1-10 

TEP 

F 

Time  (hours)  of 

following 

11-20 

QE  ( 1 ) 

F 

Inflow  (ft3/ sec) 

for  first 

• 

• 

• 

inflow  point. 

• 

71-80 

• 

• 

QE<7) 

• 

• 

F 

Inflow  (ft3/sec) 
input  point. 

for  seventh 

(use  additional  cards  If  NJSW>7) 


Note:  *M2  and  N1  cards  are  ordered  together  through  time  with  M2  cards  first 
if  TEP  =  TEB. 


A10 


TABLE  A2,  continued 
BASIC  INPUT  DATA  DESCRIPTION 


Card  Group  Card 
Card  Number  Columns 


Variable 

Name  Type  _ Variable  Description 


Ms  HISTORICAL  BOUNDARY  STAGE  DATA 


Data  for  boundary  nodes  are  entered  in  the  same  order  as  in  card  group  G. 

Time  and  heads  (feet)  referenced  to  output  datum  for  NBOUND  boundary  junctions 
with  seven  values  on  first  card  and  eight  values  on  subsequent  cards  as  needed 
if  NBOUND  is  greater  than  seven. 


Nl* 


1-10  TEB  F 


11-20  HEB (1,2)  F 
21-30  HEB (2,2)  F 
31-40  HEB (3,2)  F 


Time  (hours)  of  following  head 
inputs  when  switching  from 
curve  fit  to  historical 
inputs.  TEB  must  be  >  the 
start  of  first  day  as 
specified  by  NQSWRT.  The  time 
of  final  changes  must  be  >  end 
of  simulation. 

Head  for  first  boundary. 

Head  for  second  boundary. 

Head  for  third  boundary. 


71-80 


HEB (7, 2) 


F 


Head  for  seventh  boundary. 


Note:  *M2  and  N1  cards  are  ordered  together  through  time  with  M2  cards  first 

if  TEP  =  TEB. 


All 


TABLE  A3 

GRID  FILE  DESCRIPTION 
Card  Group  Card  Variable 

Card  Number  Columns  Name  Type  _ Variable  Description 


J:  JUNCTION  CARDS 

(Note:  Junction  surface  area  Is  based  on  triple  trapezoid  shown  In  Figure  A2) 


1-5 

J 

I 

Junction  number;  99999 
terminates  junction  input. 
Numbers  need  not  be 
continuous. 

6-10 

HEAD 

F 

Initial  or  starting  water 
surface  elevation  (feet) 
referenced  to  output  datum 
plane,  also  used  as  a  minimum 
elevation  reference  for 
continuous  Manning's  n,  and 
nodal  area  calculations 
(upward  from  datum  plane  is 
positive) . 

11-20 

SURF 

F 

Junction  bottom  surface  area 
( f  t^lO4) 

21-30 

QF1 

F 

Constant  diversion  (flow)  into 
junction  (1000  cfs). 

31-40 

DT 

F 

Depth  of  junction  (feet) 
referenced  to  input  datum. 

41-50 

SHO 

F 

Height  to  first  area 
breakpoint  (feet  -  see  Figure 
A2). 

51-60 

SAN 

F 

Increase  in  surface  area 
(ft^lO4)  per  foot  of  head 
over  first  interval  (see 

Figure  A2 ) . 

61-65 

SHI 

F 

Height  to  second  breakpoint 
( feet  -  see  Figure  A2 ) . 

66-70 

SAN1 

F 

Increase  in  surface  area 
(ft^lO4)  per  foot  of  head 
above  second  breakpoint  (see 
Figure  A2 ) . 

71-75 

XI 

F 

x  -  coordinate  (inches) 
(converted  to  feet  using 
ISCALE) . 

76-80 

Y1 

F 

y  -  coordinate  (inches) 
(converted  to  feet  using 
ISCALE) . 
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TABLE  A3,  continued 
GRID  FILE  DESCRIPTION 


Card  Group  Card 
Card  Number  Columns 


Variable 

Name  Type  _ Variable  Description 


K:  CHANNEL  CARDS 


(Note:  Channel  area  based  on  double  trapezoid  shown  in  Figure  A-3 . ) 


K1  1-5 

6-10 

11-15 

16-20 

21-25 

26-35 

36-45 

46-55 

56-65 

66-70 

71-75 

80 


N  I 

N1  I 

N2  1 

ZZ1  F 

ZZ2  F 

ALEN  F 

WIDTH  F 

DEPTH  F 

COEF  F 

CHCA  F 

VEL  F 

INDIC8  I 


Channel  number  (99999 
terminates  input ) .  Numbers 
need  not  be  continuous. 
Junction  at  one  end  of  channel 
(order  not  important). 

Junction  at  other  end  of 
channel . 

Inverse  side  slope  over  first 
interval  ( see  Figure  A3 ) . 
Inverse  side  slope  over  second 
interval  ( see  Figure  A3 ) . 
Length  of  channel  (feet).  If 
unspecified  or  zero,  program 
will  compute  from  coordinates 
of  end  junctions. 

Bottom  width  of  channel 
(feet).  Diameter  of  circular 
culvert  (INDIC8=3) . 

Depth  from  input  datum  to 
channel  bottom  (feet  - 
positive  downward).  Depth  to 
invert  of  culvert  ( INDIC8=2  or 
3). 

Manning's  coefficient. 

Height  to  cross-section 
breakpoint  (feet  -  see  Figure 
A-3 ) .  Elevation  of  weir  crest 
( INDIC8=1 ) .  Height  of  box 
culvert  ( INDIC8=2 ) . 

Initial  velocity  (ft/sec). 

Flow  (103ft3/s)  of  pumping 
across  weir  structure 
( INDIC8-1 ) . 

0  -  regular  channel 

1  -  broad-crested  weir 

2  -  box  culvert 

3  -  circular  culvert 
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B3 


PASSAIC  RIVER  FLOOO  PROTECTION 

STUDY 

OCTOBER  TIDE  11-13  1976  RAN 

FEBRUARY 

1992 

FORCING  AT 

NOOES  1, 

2.  3,  AND  4. 

APPROXIMATE  ACTUAL  FLOWS  IN  RIVERS 

!  1  o 

0  0 

1  1 

0  0 

0  1 

O  2.3  14.2  7.3  24.9 

3  25 

1.  60. 

0.0  0 

1 

1  3 

4 

1  4 

50  1 

1 

25. 

1  4 

50  1 

2 

25. 

1  4 

50  1 

3 

25. 

1  4 

50  1 

4 

25. 

125000 

0.25 

0. 

1.00 

19 

12  62 

168  175 

221  226 

167  182 

242  195 

200  206  211 

213  43  35 

136  137 

138 

229  214 

209 

0. 

1000. 

100. 

30. 

0.0 

•0.28 

-0.28 

-0.28 

2.67 

40.1 

744. 

74. 

25. 

1.0 

-1.020 

-1.020 

-1.020 

2.62 

2.0 

-1.650 

-1.650 

•1.650 

1.85 

3.0 

-1.92 

-1.920 

-1.920 

0.56 

4.0 

-1.570 

-1.570 

-1.570 

-0.84 

5.0 

-0.510 

-0.510 

-0.510 

-1.97 

6.0 

0.840 

0.840 

0.840 

•2.68 

7.0 

1.970 

1.970 

1.970 

-2.65 

8.0 

2.720 

2.720 

2.720 

-1.63 

9.0 

3.070 

3.070 

3.070 

-0.23 

10.0 

2.890 

2.890 

2.890 

1.44 

11.0 

2.160 

2.160 

2.160 

2.96 

12.0 

0.990 

0.990 

0.990 

3.91 

13.0 

-0.160 

-0.160 

-0.160 

4.04 

14.0 

-1.040 

-1.040 

-1.040 

3.23 

15.0 

-1.490 

-1.490 

-1.490 

1.88 

16.0 

-1.360 

-1.360 

-1.360 

0.29 

17.0 

-0.640 

-0.640 

-0.640 

-1.18 

18.0 

0.400 

0.400 

0.400 

-2.29 

19.0 

1.350 

1.350 

1.350 

-2.73 

20.0 

2.080 

2.080 

2.080 

-2.28 

21.0 

2.510 

2.510 

2.510 

-1.07 

22.0 

2.550 

2.550 

2.550 

0.38 

23.0 

2.120 

2.120 

2.120 

1.90 

24. 

1.300 

1.300 

1.300 

3.11 

25. 

0.420 

0.420 

0.420 

3.71 

26. 

-0.450 

-0.450 

-0.450 

3.46 

27. 

-1.110 

-1.110 

-1.110 

2.41 

28. 

-1.110 

-1.110 

-1.110 

1.05 

29. 

-0.360 

-0.360 

-0.360 

-0.35 

30. 

0.810 

0.810 

0.810 

-1.41 

31. 

1.870 

1.870 

1.870 

-1.96 

32. 

2.660 

2.660 

2.660 

-1.70 

33. 

3.160 

3.160 

3.160 

-0.57 

34. 

3.310 

3.310 

3.310 

0.82 

35. 

2.920 

2.920 

2.920 

2.32 

36. 

2.070 

2.070 

2.070 

3.52 

37. 

1.010 

1.010 

1.010 

4.20 

38. 

0.060 

0.060 

0.060 

4.03 

39. 

-0.690 

-0.690 

-0.690 

3.08 

40. 

-1.000 

-1.000 

-1.000 

1.65 

41. 

-0.680 

-0.680 

-0.680 

0.10 

73. 

300. 

30. 

10. 

42. 

0.240 

0.240 

0.240 

-1.20 

43. 

1.260 

1.260 

1.260 

-2.01 

44. 

1.970 

1.970 

1.970 

-2.15 

45. 

2.420 

2.420 

2.420 

-1.37 

46. 

2.720 

2.720 

2.720 

-0.17 

47. 

2.650 

2.650 

2.650 

1.31 

48. 

2.170 

2.170 

2.170 

2.59 

49.0 

1.300 

1.300 

1.300 

3.55 

50.0 

0.240 

0.240 

0.240 

3.83 

51.0 

-0.590 

-0.590 

-0.590 

3.13 

52.0 

-0.900 

-0.900 

-0.900 

1.94 

53.0 

-0.530 

-0.530 

-0.530 

0.57 

54.0 

0.320 

0.320 

0.320 

-0.76 

55.0 

1.370 

1.370 

1.370 

-1.64 

56.0 

2.270 

2.270 

2.270 

-1.91 
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57.0 

2.880 

2.880 

2.880 

-1.27 

58.0 

3.250 

3.250 

3.250 

-0.07 

59.0 

3.280 

3.280 

3.280 

1.34 

60.0 

2.800 

2.800 

2.800 

2.71 

61.0 

1.910 

1.910 

1.910 

3.74 

62.0 

0.930 

0.930 

0.930 

4.21 

63.0 

0.160 

0.160 

0.160 

3.86 

64.0 

-0.360 

-0.360 

-0.360 

2.72 

65.0 

-0.500 

-0.500 

-0.500 

1.31 

66.0 

-0.270 

-0.270 

-0.270 

-0.08 

67.0 

0.320 

0.320 

0.320 

-1.14 

68.0 

1.060 

1.060 

1.060 

-1.84 

69.0 

1.650 

1.650 

1.650 

-1.73 

70.0 

2.210 

2.210 

2.210 

-0.97 

B5 


B6 


1 

990. 

104.30 

-1.8 

11.48  5.66 

2 

863. 

95.30 

-1.8 

13.78  7.58 

3 

1184. 

59.30 

-1.8 

15.99  9.38 

4 

2548. 

105.65 

-3.15 

28.6627.16 

5 

381. 

U. 30 

-1.8 

8.68  7.60 

6 

457. 

53.30 

-1.8 

11.69  9.90 

7 

496. 

42.30 

-2.30 

34.040 

14.0511.36 

8 

217. 

29.30 

-2.30 

0.620 

12.2111.75 

9 

115. 

27.30 

-1.8 

11.1911.20 

10 

200. 

26.30 

-1.8 

9.8110.75 

11 

446. 

44.30 

-1.8 

8.88  9.10 

12 

225. 

26.29 

-1.79 

7.08  8.61 

13 

302. 

20.29 

-1.79 

6.05  8.39 

14 

21.7 

14.29 

-1.79 

6.99  7.90 

15 

133. 

14.29 

-1.79 

7.00  6.21 

16 

185. 

14.29 

-1.79 

4.96  7.00 

17 

112. 

8.28 

-1.78 

3.79  7.60 

18 

208. 

11.28 

-1.78 

2.37  7.81 

19 

266. 

22.28 

-1.78 

3.94  8.98 

20 

266. 

20.30 

-1.8 

5.1510.13 

21 

284. 

29.30 

-1.8 

6.7710.25 

22 

101. 

17.30 

-2.30 

1.290 

8.7710.87 

23 

56.42 

14.31 

-2.31 

0.440 

8.8911.79 

24 

156. 

20.31 

-2.31 

1.720 

7.5711.68 

25 

59.6 

22.31 

-2.31 

0.640 

6.8311.85 

26 

100. 

29.31 

-1.81 

6.4511.29 

27 

43.4 

17.34 

-2.34 

1.810 

7.0912.52 

28 

68.3 

44.20 

-1.7 

6.2112.25 

29 

99. 

9.20 

-2.20 

3.760 

5.1611.84 

30 

78. 

7.24 

-2.24 

2.400 

4.4011.03 

31 

133. 

12.28 

-2.28 

0.620 

3.2210.21 

32 

146. 

20.32 

-2.32 

0.870 

1.73  9.55 

33 

104. 

20.37 

-1.87 

0.88  8.78 

34 

114. 

7.37 

-1.87 

0.86  7.78 

35 

91. 

14.43 

-1.93 

0.08  9.25 

36 

23.8 

20.44 

-1.94 

0.60  9.98 

37 

30.4 

20.46 

-1.96 

0.3910.82 

38 

14.1 

20.48 

-2.48 

1.450 

1.1511.25 

39 

7.8 

20.49 

-2.49 

0.540 

1.5711.70 
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100. 

20. 

5.5 

.025 

10.0 

359 

198 

201 

900. 

0.0 

600. 

-9.0 

.055 

10.5 

360 

199 

201 

11.1 

630. 

4250. 

300. 

22.8 

.03 

2.5 

361 

199 

200 

6.4 

300. 

4800. 

30. 

15.1 

.035 

8.5 

362 

197 

202 

0.0 

100. 

10000. 

20. 

2.5 

.055 

5.0 

363 

198 

202 

9000. 

1000. 

-6.5 

.055 

364 

202 

203 

2000. 

-6.5 

.055 

365 

198 

203 

1100. 

-4.5 

.055 

366 

201 

203 

12.0 

247. 

5900. 

225. 

20.5 

.03 

1.5 

367 

201 

204  6.67 

0.0 

7400. 

40. 

5.5 

.035 

0.5 

368 

192 

193 

4.0 

0.0 

6500. 

30. 

18.5 

.03 

19.5 

369 

202 

205 

100. 

.055 

370 

203 

206 

12.0 

200. 

7800. 

300. 

18.3 

.035 

2.5 

371 

204 

207 

0.0 

120. 

6000. 

40. 

.035 

7.0 

372 

203 

204 

100. 

-5.5 

.055 

373 

164 

165 

3000. 

20. 

.03 

374 

205 

206 

1000. 

-7.5 

.055 

375 

206 

207  3.22 

233. 

7900. 

200. 

.03 

11.5 

376 

207 

208  6.25 

7.7 

6900. 

200. 

.03 

1.5 

377 

208 

209 

7.0 

37.5 

7100. 

100. 

15.3 

.03 

4.0 

378 

206 

210  13.3 

100. 

8200. 

200. 

12.3 

.035 

4.5 

379 

210 

211 

0.0 

970. 

6800. 

130. 

12.3 

.035 

12.0 

380 

211 

212  25.0 

14.0 

12550. 

100. 

.04 

4.5 

381 

212 

213 

5.0 

6.25 

10000. 

100. 

.04 

4.5 

382 

176 

177  14.7121.9 

165. 

13.6 

.035 

6.5 

383 

177 

215 

11.0131.1 

3600. 

175. 

15.7 

.035 

6.5 

384 

215 

216 

15.0121.8 

200. 

13.6 

.035 

6.0 

385 

216 

217  16.0153.8 

2800. 

165. 

12.3 

.03 

6.5 

386 

217 

218 

18.0115.8 

3140. 

190. 

13.6 

.03 

7.5 

387 

218 

219 

18.0  73.8 

2900. 

150. 

13.9 

.03 

9.0 

388 

219 

220 

12.0  92.5 

3300. 

150. 

15.5 

.03 

9.5 

389 

220 

221 

19.0108.5 

2700. 

150. 

13.6 

.03 

9.5 

390 

221 

222 

15.0122.0 

3100. 

150. 

11.1 

.04 

9.5 

391 

222 

223 

16.0139.5 

2700. 

150. 

14.2 

.04 

10.0 

392 

223 

224 

12.0141.6 

2300. 

140. 

14.2 

.04 

10.0 

393 

224 

225 

9.5111.5 

3100. 

165. 

.04 

9.5 

394 

225 

226 

10.5 

83.5 

2600. 

135. 

.03 

10.0 

395 

226 

227 

0. 

50. 

3050. 

250. 

.03 

10. 

396 

227 

228 

16. 

24. 

2100. 

200. 

.03 

10. 

397 

228 

229 

16. 

24. 

2050. 

200. 

.03 

10. 

398 

41 

230 

0.0 

400. 

6400. 

466. 

14.0 

.04 

0.5 

399 

230 

231 

0.0 

350. 

5417. 

387. 

12.3 

.04 

0.5 

400 

231 

232 

0.0 

250. 

5200. 

365. 

7.3 

.04 

0.5 

401 

232 

233 

0.0 

150. 

6990. 

252. 

4.6 

.04 

0.5 

402 

233 

234 

0.0 

100. 

3199. 

97. 

2.7 

.04 

0.5 

403 

233 

235 

0.0 

100. 

4202. 

97. 

2.7 

.04 

0.5 

404 

235 

236 

0.0 

50. 

4197. 

50. 

2.5 

.04 

0.5 

405 

235 

237 

0.0 

50. 

3000. 

50. 

-1.6 

.04 

2.0 

406 

190 

191 

100. 

-6.5 

.055 

407 

241 

189 

2000. 

-7.5 

.055 

408 

184 

241 

2000. 

-7.5 

.055 

409 

240 

185 

7.8 

220. 

5000. 

10. 

9.5 

.035 

-0.5 

B15 


410 

246 

192 

600. 

411 

239 

197  8.6  47.2 

4000. 

412 

205 

238 

413 

238 

210  350.  900. 

0. 

414 

240 

181  800.  0.1 

0. 

415 

242 

243 

600. 

416 

243 

198 

417 

243 

197 

418 

243 

192 

419 

199 

244  0.0  300. 

3000. 

420 

245 

198  0.0  300. 

4200. 

421 

194 

246  14.3  66.2 

5900, 

422 

165 

166 

423 

144 

164  5. 

4065. 

424 

144 

165  5. 

425 

144 

166  30. 

426 

144 

163  20. 

427 

143 

163 

428 

143 

144 

429 

46 

143 

430 

44 

144  7. 

431 

42 

269 

7260 

432 

43 

270 

2900, 

433 

44 

46 

6200 

434 

45 

47  15. 

4200 

435 

48 

50  10. 

5600 

436 

170 

271  18.6 

437 

183 

272  9.6 

4000 

438 

225 

273 

4880 

439 

273 

274 

4880 

99999 

10. 

5.5 

.025 

9.0 

80. 

.055 

-0.5 

100. 

300. 

-8.5 

.055 

.055 

-0.5 

2000. 

-7.0 

.055 

9.5 

7. 

5.0 

.025 

7.0 

100. 

-9.5 

.055 

200. 

-9.5 

.055 

50. 

-9.5 

.055 

40. 

2.5 

.04 

2.0 

40. 

2.5 

.055 

2.0 

40. 

5.5 

.04 

1.5 

3200. 

14. 

.035 

500. 

35. 

.03 

1500. 

40. 

.03 

1100. 

12. 

.035 

2000. 

30. 

.03 

2200. 

8. 

.035 

2000. 

7. 

.035 

3300. 

6. 

.035 

2700. 

40. 

.03 

900. 

35. 

.03 

900. 

35. 

.03 

2500. 

37. 

.03 

1563. 

40. 

.03 

1563. 

40. 

.03 

240. 

17.2 

.035 

400. 

39.6 

.03 

30. 

-2. 

.04 

53. 

-6. 

.04 

B16 


TIME  NODE  12  NOOE  62  N00E168  NOOE175  NOOE221  NOOE226  N00E167  NOOE182  NCOE242  NOOE195 


NOOE2QO  N0CE206  NOOE211  NOOE213  NOOE  43  NOOE  35  NOOE136  N00E137  NOOE138  NOOE 


.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.50 

-.22 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.04 

.00 

.00 

.00 

.00 

.00 

1.00 

-.74 

-.09 

.00 

.00 

.03 

.14 

-.01 

.00 

.00 

.00 

.00 

.01 

.02 

.05 

-.02 

-.02 

.00 

.00 

.00 

1.50 

-1.00 

-.63 

-.12 

.03 

.08 

.51 

-.24 

-.04 

.00 

.01 

.01 

.01 

.01 

.05 

-.30 

-.32 

-.21 

-.06 

-.01 

2.00 

-1.47 

-.85 

-.64 

-.14 

.13 

1.03 

-.72 

-.49 

-.23 

-.13 

-.17 

-.04 

.00 

.02 

-.78 

-.93 

-.81 

-.44 

-.22 

2.50 

-1.77 

-1.07 

-1.17 

-.62 

-.17 

1.05 

-1.23 

-1.07 

-.84 

-.77 

-.79 

-.49 

-.25 

.00 

-1.29 

-1.39 

-1.31 

-.94 

-.83 

3.00 

•1.88 

-1.51 

-1.67 

-1.13 

-.69 

.87 

-1.67 

-1.65 

-1.51 

-1.44 

-1.47 

-1.21 

-.84 

-.21 

-1.65 

-1.75 

-1.71 

-1.49 

-1.43 

3.50 

-2.03 

-1.80 

-1.99 

-1.64 

-1.21 

.77 

-1.95 

-2.07 

-2.11 

-2.09 

-2.11 

-2.00 

-1.58 

-.55 

-1.95 

•2.04 

-2.04 

-1.92 

-1.90 

4.00 

-1.87 

-2.10 

-2.23 

-1.99 

-1.67 

.78 

-2.23 

-2.28 

-2.40 

-2.45 

-2.43 

•2.48 

-2.23 

-.94 

-2.22 

-2.18 

•2.20 

-2.22 

-2.25 

4.50 

-1.65 

-2.08 

-2.45 

-2.28 

-2.01 

.85 

-2.35 

-2.50 

-2.54 

-2.52 

-2.54 

•2.62 

-2.58 

-1.32 

-2.28 

-2.24 

-2.29 

-2.37 

-2.43 

5.00 

-1.09 

-1.84 

-2.26 

-2.49 

-2.29 

.92 

-2.13 

-2.41 

-2.58 

-2.64 

-2.61 

-2.60 

-2.63 

-1.66 

-2.08 

-1.97 

-2.09 

-2.35 

-2.44 

5.50 

-.39 

-1.44 

-1.81 

-2.20 

-2.44 

1.00 

-1.70 

-1.89 

-2.08 

-2.17 

-2.14 

-2.35 

-2.49 

-1.94 

-1.61 

-1.47 

-1.57 

-1.95 

-2.16 

6.00 

.29 

-.63 

-1.10 

-1.68 

-1.94 

1.01 

-.96 

-1.26 

-1.47 

-1.52 

-1.51 

-1.73 

-1.97 

-2.10 

-.88 

-.79 

-.94 

-1.43 

-1.54 

6.50 

1.06 

.03 

-.21 

-.88 

-1.11 

.79 

-.15 

-.32 

-.58 

-.73 

-.64 

-.93 

-1.25 

-1.80 

-.12 

-.03 

-.13 

-.64 

-.83 

7.00 

1.61 

.86 

.59 

.08 

-.17 

.73 

.66 

.52 

.34 

.26 

.30 

-.01 

-.39 

-1.20 

.70 

.74 

.61 

.14 

.03 

7.50 

2.22 

1.45 

1.40 

.97 

.79 

1.12 

1.41 

1.33 

1.12 

.83 

1.07 

.78 

.43 

-.53 

1.42 

1.55 

1.48 

1.04 

.90 

8.00 

2.60 

2.15 

2.04 

1.82 

1.70 

1.81 

2.08 

1.98 

1.83 

1.23 

1.76 

1.51 

1.19 

.17 

2.12 

2.22 

2.15 

1.87 

1.83 

8.50 

2.98 

2.61 

2.69 

2.52 

2.48 

2.55 

2.70 

2.61 

2.39 

1.62 

2.30 

2.07 

1.84 

is 

2.70 

2.83 

2.83 

2.67 

2.64 

9.00 

3.18 

2.95 

3.18 

3.18 

3.14 

3.17 

3.16 

3.13 

2.94 

2.03 

2.82 

2.58 

2.37 

1.66 

3.16 

3.22 

3.23 

3.26 

3.33 

9.50 

3.27 

3.19 

3.50 

3.59 

3.70 

3.76 

3.45 

3.47 

3.34 

2.48 

3.27 

3.01 

2.86 

2.32 

3.42 

3.52 

3.58 

3.67 

3.70 

10.00 

3.19 

3.30 

3.62 

3.87 

3.96 

4.05 

3.56 

3.65 

3.61 

2.91 

3.57 

3.35 

3.25 

2.91 

3.54 

3.54 

3.60 

3.76 

3.85 

10.50 

2.97 

3.30 

3.59 

3.79 

3.96 

4.10 

3.53 

3.65 

3.70 

3.29 

3.75 

3.61 

3.58 

3.39 

3.47 

3.40 

3.46 

3.58 

3.64 

11.00 

2.55 

3.03 

3.38 

3.57 

3.68 

3.86 

3.28 

3.53 

3.70 

3.57 

3.73 

3.76 

3.79 

3.76 

3.23 

3.03 

3.07 

3.20 

3.24 

11.50 

2.02 

2.66 

2.97 

3.21 

3.36 

3.55 

2.86 

3.18 

3.51 

3.67 

3.62 

3.77 

3.87 

3.96 

2.77 

2.57 

2.61 

2.73 

2.76 

12.00 

1.42 

2.12 

2.42 

2.76 

2.96 

3.22 

2.29 

2.67 

3.10 

3.49 

3.29 

3.55 

3.71 

3.93 

2.22 

2.02 

2.08 

2.24 

2.27 

12.50 

.78 

1.57 

1.78 

2.21 

2.50 

2.82 

1.70 

1.96 

2.44 

3.23 

2.66 

3.16 

3.33 

3.63 

1.65 

1.48 

1.53 

1.69 

1.76 

13.00 

.21 

.93 

1.16 

1.62 

1.93 

2.36 

1.07 

1.33 

1.75 

2.93 

1.79 

2.50 

2.80 

3.25 

1.01 

.93 

.99 

1.18 

1.22 

13.50 

-.37 

.40 

.49 

1.02 

1.36 

1.88 

.44 

.62 

1.03 

2.63 

1.15 

1.76 

2.13 

2.80 

.41 

.37 

.42 

.62 

.70 

14.00 

-.81 

-.14 

-.11 

.39 

.77 

1.46 

-.15 

-.04 

.25 

2.32 

.33 

.94 

1.40 

2.28 

-.17 

-.17 

-.10 

.11 

.16 

14.50 

-1.20 

-.60 

-.66 

-.18 

.19 

1.09 

-.67 

-.62 

-.41 

1.96 

-.35 

.09 

.61 

1.73 

-.67 

-.69 

-.63 

-.42 

-.35 

15.00 

-1.47 

-1.01 

-1.10 

-.72 

-.35 

.84 

-1.10 

-1.13 

-1.03 

1.26 

-1.00 

-.77 

-.17 

1.17 

-1.11 

-1.13 

-1.09 

-.90 

-.85 

15.50 

-1.61 

-1.33 

-1.51 

-1.16 

-.84 

.69 

-1.49 

-1.55 

-1.58 

-1.15 

-1.59 

-1.51 

-.95 

.62 

-1.48 

-1.50 

-1.48 

-1.33 

-1.30 

16.00 

-1.61 

-1.52 

-1.81 

-1.57 

-1.26 

.64 

-1.76 

-1.89 

-1.95 

-1.83 

-1.96 

-1.96 

-1.60 

.09 

-1.73 

-1.78 

-1.78 

-1.70 

-1.68 

16.50 

-1.43 

-1.64 

-1.93 

-1.87 

-1.64 

.64 

-1.89 

-2.00 

-2.10 

-2.20 

-2.13 

-2.19 

-2.02 

-.41 

-1.87 

-1.89 

-1.93 

-1.96 

-1.99 

17.00 

-1.11 

-1.50 

-1.90 

-2.00 

-1.91 

.69 

-1.79 

-1.98 

-2.07 

-2.12 

-2.09 

-2.16 

-2.17 

-.86 

-1.75 

-1.77 

-1.88 

-2.06 

-2.14 

17.50 

-.58 

-1.25 

-1.58 

-1.91 

-2.02 

.73 

-1.51 

-1.66 

-1.83 

-1.86 

-1.87 

-1.96 

-2.04 

-1.24 

-1.46 

-1.34 

-1.45 

-1.81 

-2.02 

18.00 

.00 

-.68 

-1.10 

-1.49 

-1.73 

.73 

-.96 

-1.21 

-1.33 

-1.34 

C4 


-1.35 

-1.49 

-1.63 

-1.50 

-.87 

-.80 

-.93 

-1.34 

-1.47 

18.50 

.60 

-.14 

-.33 

-.88 

-1.04 

.61 

-.28 

-.42 

-.66 

-.78 

-.72 

-.92 

-1.10 

-1.44 

-.25 

-.13 

-.22 

-.66 

-.79 

19.00 

1.09 

.53 

.37 

-.02 

-.25 

.60 

.42 

.29 

.15 

.09 

.12 

-.14 

-.44 

-.93 

.46 

.52 

.43 

.07 

-.02 

19.50 

1.60 

1.04 

1.04 

.74 

.61 

.95 

1.05 

1.00 

.84 

.69 

.81 

.57 

.29 

-.39 

1.06 

1.17 

1.13 

.82 

.74 

20.00 

1.99 

1.62 

1.59 

1.45 

1.39 

1.54 

1.61 

1.56 

1.46 

1.07 

1.44 

1.23 

.98 

.19 

1.63 

1.69 

1.65 

1.50 

1.49 

20.50 

2.33 

2.06 

2.12 

2.04 

2.04 

2.13 

2.11 

2.08 

1.96 

1.41 

1.90 

1.75 

1.56 

.84 

2.11 

2.19 

2.20 

2.12 

2.10 

21.00 

2.58 

2.44 

2.53 

2.55 

2.56 

2.62 

2.52 

2.50 

2.39 

1.75 

2.35 

2.19 

2.04 

1.48 

2.53 

2.56 

2.57 

2.57 

2.61 

21.50 

2.72 

2.68 

2.87 

2.91 

2.97 

3.05 

2.85 

2.84 

2.75 

2.10 

2.71 

2.57 

2.46 

2.05 

2.84 

2.86 

2.89 

2.93 

2.94 

22.00 

2.76 

2.76 

3.07 

3.21 

3.23 

3.31 

3.01 

3.09 

3.05 

2.45 

3.01 

2.88 

2.81 

2.55 

2.99 

3.00 

3.03 

3.10 

3.15 

22.50 

2.66 

2.82 

3.07 

3.25 

3.39 

3.49 

3.02 

3.12 

3.19 

2.77 

3.23 

3.14 

3.11 

2.96 

2.99 

2.99 

3.03 

3.13 

3.17 

23.00 

2.42 

2.69 

2.96 

3.13 

3.25 

3.43 

2.89 

3.07 

3.21 

3.03 

3.23 

3.28 

3.31 

3.29 

2.86 

2.77 

2.82 

2.94 

3.00 

23.50 

2.05 

2.46 

2.71 

2.87 

2.99 

3.17 

2.64 

2.86 

3.09 

3.17 

3.16 

3.30 

3.38 

3.47 

2.58 

2.44 

2.48 

2.60 

2.64 

24.00 

1.60 

2.05 

2.34 

2.55 

2.68 

2.90 

2.23 

2.53 

2.82 

3.08 

2.93 

3.13 

3.25 

3.43 

2.17 

2.01 

2.05 

2.17 

2.20 

24.50 

1.12 

1.68 

1.82 

2.14 

2.33 

2.59 

1.75 

1.97 

2.32 

2.87 

2.45 

2.79 

2.94 

3.18 

1.70 

1.56 

1.59 

1.70 

1.75 

25.00 

.67 

1.21 

1.29 

1.64 

1.90 

2.25 

1.24 

1.40 

1.70 

2.62 

1.73 

2.24 

2.47 

2.86 

1.20 

1.10 

1.15 

1.27 

1.29 

25.50 

.21 

.77 

.77 

1.13 

1.39 

1.84 

.74 

.82 

1.06 

2.34 

1.14 

1.57 

1.88 

2.45 

.73 

.65 

.67 

.80 

.86 
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PASSAIC  RIVER  FLOOD  PROTECTION  PROJECT 
Fact  Sheet 

PURPOSE/ DESCRIPTION:  Congress  authorized  the  Passaic  River 

Flood  Protection  Project  for  construction  in  November  1990.  The 
project  is  to  provide  comprehensive  and  environmentally  sound 
flood  protection  to  40  towns  in  the  Passaic  River  Basin. 

Flooding  has  long  been  a  problem  in  the  Passaic  River 
Basin.  Since  Colonial  times  floods  have  claimed  lives  and 
damaged  property.  With  the  explosive  growth  of  residential  and 
industrial  development  in  recent  years,  the  threat  of  serious 
damage  and  loss  of  life  from  flooding  has  multiplied. 

The  most  severe  flood,  the  "flood  of  record,"  occurred  in 
1903,  and  more  recent  floods  in  1968,  1971,  1972,  two  in  1975, 
and  in  1984  were  sufficiently  devastating  to  warrant  Federal 
flood  declarations.  The  flood  of  April  5,  1984,  caused  such 
severe  flooding  in  northern  New  Jersey  that  three  people  died 
and  there  was  $410  million  in  damages  (October  1991  dollars) . 

The  Passaic  River  Basin,  935  square  miles  drained  by  the 
Passaic  River  and  its  tributaries,  lies  mostly  in  northern  New 
Jersey.  The  basin  hold  s  132  municipalities  in  ten  counties  — 
Bergen,  Essex,  Hudson,  Morris,  Passaic,  Sussex,  Somerset,  and 
Union  in  New  Jersey,  and  Orange  and  Rockland  in  New  York'  State 
—  and  has  some  of  the  most  densely  developed  areas  in  the 
nation.  More  than  2.5  mill  ion  people  live  in  the  basin  (1990 
census),  and  a  bout  30,000  homes  and  businesses  lie  in  the 
basin's  flood  plain. 

The  basin  extends  from  Newark  in  the  east  to  Lake  Hopatcong 
in  the  west,  and  from  Monroe,  N.Y.,  in  the  north  to  Somerset 
County,  N.J.,  in  the  south.  It  is  studded  with  rivers  and 
streams  including  the  Pompton,  Pequannock,  Wanaque,  Ramapo, 
Mahwah,  Saddle,  Whippany,  and  Rockaway  rivers  —  major 
tributaries  of  the  Passaic  River. 

STATUS:  The  U.S.  Army  Corps  of  Engineers  has  been  working 
on  plans  to  reduce  flooding  in  the  basin  since  1939,  but  none 
of  those  plans  was  implemented.  Congress  then  authorized  the 
Corps  to  conduct  a  new  study  of  the  Passaic  River  Basin  for  the 
State  of  New  Jersey  in  the  Water  Resources  Development  Act  of 
1976. 

The  Corps  evaluated  more  than  150  plans,  including  a  full 
range  of  non-structural  alternatives  such  as  buying  out  flood- 
prone  homes.  It  presented  an  array  of  plans  to  the  State  of 
New  Jersey.  After  a  series  of  public  hearings,  the  State 
selected  a  dual  inlet  water  diversion  tunnel  system  as  the 
center  piece  of  its  comprehensive  flood  protection  program  for 
the  basin.  On  May  8,  1987,  the  State  formally  announced  its 
intention  to  act  as  the  Passaic  River  Flood  Protection 
Project's  non-Federal  sponsor,  and  the  Corps  recommended  the 
project  in  its  February  1988  General  Design  Memorandum. 
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In  November  1990,  congress  authorized  construction  of  the 
§1.2  billion  Passaic  River  Flood  Protection  Project  in  the 
Water  Resources  Development  Act  of  1990  (Public  Law  101-640) . 
The  Federal  project  cost  is  $890  million  (October  1989 
dollars) ;  the  State  project  cost  is  of  $310  million. 

i 

The  act  changes  the  plan  originally  selected  by  the  State, 
extending  the  tunnel  to  an  outlet  in  Newark  Bay  and  eliminating 
the  need  for  nine  levee-flood-wall  systems  in  eastern  Essex  and 
southern  Bergen  counties.  Eastern  Essex  and  southern  Bergen 
counties  will  now  receive  fluvial  flood  protection  without 
having  any  surface  structures  along  the  lower  Passaic  River. 

The  act  also  authorizes  a  $6  million  stream-bank  restoration 
project  for  the  City  of  Newark. 

The  major  elements  of  the  project  are  two  underground 
tunnels,  a  20.1-mile-long  main  tunnel  about  40  feet  in  diameter 
and  a  1.2-mile  spur  tunnel  about  22  feet  in  diameter..  The  main 
tunnel  will  carry  flood  waters  from  an  inlet  on  the  upper 
Pompton  River  down  to  an  outlet  in  Newark  Bay.  The  spur  tunnel 
will  convey  Central  Basin  flood  waters  from  an  inlet  just 
downstream  of  Two  Bridges  in  Wayne  to  an  underground  connection 
with  the  main  tunnel. 

About  5.9  miles  of  channel  modifications  will  be  requirea 
to  direct  the  flows  into  the  tunnel  inlets  and  17.6  miles  of 
levees  and  6.2  miles  of  flood  walls  will  increase  flood 
protection.  The  diversion  tunnels  and  associated  surface  works 
are  generally  designed  to  protect  against  floods,  ranging  up  to 
the  100-year  event.  The  project  has  a  benefit-to-cost  ratio  of 
1.1  to  1 . 

The  Water  Resources  Act  also  recognizes  New  Jersey's 
earlier  contribution  of  lands  and  project  works  by  crediting 
these  investments  toward  the  State's  share  of  the  project's 
cost.  It  also  establishes  a  wetlands  bank  to  enable  the  State 
of  New  Jersey  to  purchase  wetlands  and  watershed  lands  that  can 
also  be  credited  toward  its  share  of  the  project's  cost. 

The  Passaic  River  Flood  Protection  Project  includes  the 
acquisition  of  5,350  acres  of  natural  water  storage  areas, 

5,200  acres  of  which  are  wetlands.  Without  the  project,  these 
lands  could  conceivably  be  developed,  worsening  existing  flood 
problems.  The  State  of  New  Jersey  has  an  agreement  with  the 
Corps  in  which  it  promises  to  continue  to  protect  6,300  acres 
that  are  annually  inundated  by  the  river,  thus  avoiding  any 
secondary  development.  About  9,500  acres  of  the  Central  Basin 
are  already , protected  as  designated  parkland,  bringing  the 
total  of  natural  storage  areas  that  will  be  permanently 
protected  with  the  project  to  21,000  acres. 

The  Passaic  River  Flood  Protection  Project  includes 
preservation  measures  for  fish  and  wildlife  and  provision  for 
wetlands  mitigation.  The  project  will  also  expand  the 
opportunities  for  recreation  along  the  river. 


CONTACT:  Mr.  Robert  A.  Pietrowsky,  Chief,  Passaic  River 

Division,  Telephone:  (201)  656-2304. 
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APPENDIX  E: 


MODELING  PROCEDURES 


Coincidence  Analysis 

El.  Files  of  daily  values  of  Passaic  River  flow  and  maximum  surge 
height  at  The  Battery  were  used  to  investigate  the  coincidence  between 
occurrences  of  flow  and  surge  magnitudes  that  exceeded  user-specified 
thresholds.  Daily  river  flow  values  were  from  USGS  data  from  the  Little 
Falls,  New  Jersey  gage.  Daily  maximum  surge  values  were  obtained  by 
subtracting  hourly  predicted  tide  heights  from  observed  water  surface 
elevations  at  The  Battery  and  selecting  the  maximum  on  each  day.  With  some 
gaps,  the  data  records  extended  from  1927  to  1992.  The  dates  of  river  flow 
and  surge  height  values  corresponded  on  a  line-by-line  basis  in  their 
respective  data  files.  The  total  number  of  record  days  was  21,174.  Cross 
correlation  calculations  were  first  performed  in  order  to  determine  if  an 
average  lead/lag  relation  existed  between  river  flow  and  surge  height,  other 
programs  were  then  used  to  count  the  joint  occurrences  of  particular 
combinations  of  river  flow  and  surge  height. 

E2.  The  following  procedure  was  used  to  examine  the  coincidence  between 
river  flow  and  surge  height: 

1.  With  the  data  files  "SURGE. DLY"  and  "RIVER. DLY"  as  input,  the 
program  "COIN-BIG"  was  used  to  tally  coincidences  between  values 
of  river  flow  and  surge  height.  The  program  queries  for  input  and 
output  file  names  and  for  threshold  values.  Threshold  values  can 
be  specified  directly  by  the  user  or  can  be  calculated  by  the 
program  by  multiplying  the  respective  data  series  standard 
deviations  by  a  user-specified  factor.  COIN-BIG  reads  a  surge 
height,  and  examines  river  flows  on  the  same  day  and  the  three 
days  before  and  after.  The  appropriate  indexes (s)  are  incremented 
if  the  surge  height  and  river  flow  on  a  given  day  both  exceed 
their  thresholds.  Thus,  seven  counts  are  maintained.  The  program 
reads  successive  daily  surge  heights  and  accumulates  counts  until 
termination  with  the  third-from-the-last  value.  Output  is  an 
ASCII  file  of  pertinent  parameters  and  counts  for  each  of  the 
seven  lead/lag  days. 

2.  A  second  program, "SURRIV" ,  was  used  to  select  particular  joint 
occurrences  of  surge  height  and  river  flow.  The  program  queries 
for  input  (usually  SURGE. DLY  AND  RIVER. DLY)  and  output  file  names 
and  for  threshold  values.  The  program  then  reads  a  surge  height. 
If  surge  height  exceeds  the  specified  threshold,  the  program 
compares  river  flows  on  the  same  day,  the  first  day  after,  and  the 
second  day  after  against  the  specified  river  flow  threshold.  If 
river  flow  exceeds  threshold  on  more  than  one  of  those  three  days, 
the  program  selects  the  maximum.  If  none  of  those  three  river 
flow  values  exceeds  the  threshold,  the  program  reads  the  next 
surge  height  (next  day)  and  examines  the  appropriate  set  of  three 
river  flow  values.  Output  is  a  file  of  line-by-line  surge  height 
and  river  flow  values  that  jointly  exceed  thresholds  accompanied 
by  their  respective  dates.  On  a  given  line,  dates  may  differ  by 
zero,  one,  or  two  days. 
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Identify  Hew  Storms 


E3.  In  this  study,  nine  new  synthetic  storm  tide  events  were  selected 
to  add  to  the  24  events  used  in  the  Hackensack  River  Surge  Barrier  study  (CDM 
1988),  for  a  total  of  33  events.  Each  synthetic  storm  tide  event  is  the 
combination  of  a  synthetic  storm  surge  (from  252  hurricanes  and  43 
northeasters),  a  tide  (from  five  tidal  marigrams),  and  a  lag  (from  50  1/2-hour 
lags).  This  resulted  in  73,750  possible  synthetic  storm  tide  combinations  to 
choose  from. 

E4.  The  following  procedure  was  used  to  identify  the  nine  additional 
events : 

1.  From  the  stage-frequency  distribution  developed  for  the  mouth  of 
Newark  Bay  (node  44)  during  the  CDM  (1983b)  study  (and  reported  in 
CDM  1988),  determine  what  storm  tide  elevation  each  new  event 
should  have.  For  this  study,  we  linearly  interpolated  between 
pairs  of  events  at  each  end  of  the  distribution. 

2.  Using  the  data  files  "HH . SUM"  and  "NE.SUM",  which  contain  maximum 
storm  tide  heights  at  the  mouth  of  Newark  Bay  (node  44)  for  the 
hurricane  and  northeaster  combinations,  respectively,  run  the 
program  "PICKEM.EXE"  for  each  tidal  elevation  identified  and  a 
small  tolerance  (say,  0.001  feet). 

3.  The  output  file  produced,  "PICKEM.OUT" ,  contains  storm  tide 
combinations  that  produced  maximum  tidal  elevations  within  the 
specified  tolerance  of  the  selected  tidal  elevation.  From  each 
group  of  candidate  events,  we  selected  one  event.  Each  selection 
was  made  to  provide  an  overall  mix  of  hurricanes  and  northeasters, 
and  to  distribute  the  storm  types  over  the  range  of  storms 
available.  The  selected  storm  tide  events  are  listed  in  the  file 
"PICKEM1 .OUT" .  The  output  files  list  the  identification  number  of 
the  synthetic  storm  surge,  the  identification  number  of  the 
coincident  tide,  and  the  lag  as  the  number  of  half-hour  intervals 
(with  "1"  being  zero  lag). 

Develop  Synthetic  Storm  Tide  Events 

E5.  Once  a  synthetic  storm  tide  event  had  been  selected,  the  next  step 
was  to  create  the  storm  tide  boundary  condition  marigrams  at  nodes  1-4  of  the 
link-node  model  of  the  interior  system  (including  New  York  Harbor  and  Newark 
Bay).  From  the  CDM  (1983a  and  1983b)  studies  of  storm  tides  in  the  New  York 
City  area,  boundary  condition  marigrams  for  the  252  synthetic  hurricane  surges 
and  43  synthetic  northeaster  surges,  developed  using  the  offshore  surge  model, 
were  saved  on  21  magnetic  tapes.  Table  El  gives  the  contents  of  each  tape. 
These  tapes  reside  at  the  Coastal  Engineering  Research  Center  (CERC), 

Waterways  Experiment  Station,  in  Vicksburg,  Mississippi. 

E6.  The  following  procedure  was  used  to  create  the  boundary  condition 
storm  tide  marigram  file  for  the  interior  link-node  model,  for  each  storm  tide 
combination  identified: 

1.  A  "combination"  input  file  was  created,  with  the  "tidal  day" 

(equal  to  25.0  hours)  and  the  identified  lag  (in  hours)  on  the 
first  line,  and  the  times  and  tidal  elevations  of  the  two  high 
tides  and  two  low  tides,  for  each  link-node  boundary  junction  1-4, 
on  lines  2-5.  Nodes  1-3  were  based  on  tides  at  Sandy  Hook,  and 
node  4  was  based  on  tides  at  Bridgeport,  Connecticut.  Tides  that 
were  combined  with  hurricane  and  northeaster  surges  are  shown  in 
Tables  E2  and  E3,  respectively. 

2.  The  identified  hurricane  or  northeaster  surge  file  (the  results 
from  the  offshore  model  saved  at  cells  corresponding  to  the 
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interior  model's  boundary  nodes  1-4),  was  downloaded  to  the  VAX 
computer  from  the  appropriate  9-track  magnetic  tape  (see  Table  El 
for  cross-reference) . 

3a.  For  hurricane  surge  tide  combinations,  the  program  "COMBINE.EXE" 
was  run  on  the  VAX  to  produce  a  binary  output  file  on  the  VAX. 
Alternatively,  the  program  "COMB15.EXE"  was  run  on  the  VAX  to 
produce  an  ASCII  output  file,  with  a  time  interval  of  15  minutes 
(to  reduce  data  storage),  that  could  be  transferred  from  the  VAX 
to  a  PC.  All  hurricane  storm  tide  simulations  for  this  study  were 
performed  on  a  386  PC. 

3b.  For  northeaster  storm  tide  combinations,  the  program 

NE_COMBINE2.EXE"  was  run  to  produce  a  binary  output  file  on  the 
VAX.  All  northeaster  simulations  were  performed  on  the  VAX. 

3c.  Each  of  the  combination  programs  required  an  input  control  files 
(step  #1)  and  an  offshore  surge  results  file  (step  #2).  The 
models  each  output  a  binary  tidal  boundary  conditions  file  (except 
an  ASCII  file  from  "COMB15.EXE") ,  and  a  results  output  file  that 
echoes  the  steps  the  program  took  in  performing  the  combination. 

4.  ASCII  tidal  boundary  condition  files  were  then  transferred  to  a  PC 
using  standard  transfer  software,  such  as  "KERMIT". 

Link-Node  Model  simulations 

E7 .  The  input  files  required  for  a  simulation  using  the  link-node  model 
are  described  in  detail  in  the  users  manual  (Appendix  A) ,  and  shown  in  Figure 
Al.  The  output  files  produced  are  also  discussed  in  Appendix  A  and  shown  in 
Figure  Al. 

E8.  The  following  procedure  was  used  to  perform  a  simulation  using  the 
link-node  model: 

1.  The  input  (or  "control")  and  "grid"  files  are  prepared  as 

described  in  the  users  manual  (Appendix  A) .  The  tidal  boundary 
conditions  file  is  either  prepared  as  discussed  in  the  above 
section  or  specified  as  a  blank  (empty)  file  if  tidal  boundary 
data  were  included  in  the  "Control"  (or  "Input")  file  on  unit  N5. 
If  the  tidal  boundary  conditions  file  is  not  empty  (i.e.,  a 
synthetic  hurricane  or  northeaster  storm  tide  is  being  simulated), 
then  the  file  should  be  a  binary  file  if  running  on  the  VAX  or  an 
ASCII  file  if  running  on  a  PC. 

2a.  On  the  VAX,  the  link-node  model  is  executed  interactively  by 

typing  "R  DEMNEWARK" .  The  model  then  requests  the  names  of  the 
input  and  output  files  in  the  order  described  in  step  2c.  (A 
batch  file  can  be  set  up  on  the  VAX.) 

2b.  On  the  PC,  the  link-node  model  is  executed  by  typing  "DEMNEW". 

The  model  then  requests  the  names  of  the  input  an  output  files  in 
the  order  described  in  step  2c. 

2c.  The  link-node  model  will  ask  for  the  following  input  and  output 
file  names: 

•  Control  or  input  file  name 

•  Tabular  output  file  name 

•  Heads  output  file  name 

•  Tidal  boundary  conditions  file  name 

•  Link-node  grid  file  name 

3.  Once  the  model  has  executed,  the  computed  heads  in  the  "heads 
output"  file  can  be  plotted  using  a  standard  software  package, 
such  as  "GRAPHER"  or  "HARVARD  GRAPHICS". 
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Stage-Frequency  Curve* 


E9.  Once  the  33  synthetic  hurricane  and  northeaster  storm  tide 
conditions  have  been  simulated,  stage  frequency  curves  at  any  nodes  whose 
computed  heads  are  output  to  the  "heads  output"  file  can  be  developed  as 
follows: 

1.  The  names  of  the  33  "heads  output"  files  are  listed  in  a  "file 
names"  input  file.  The  order  of  the  file  names  corresponds  to  the 
order  of  the  probabilities  of  occurrence  in  the  program 
"FREQ33.FOR". 

2.  A  second  input  file  is  created  that  contains  the  number  of  nodes 
whose  computed  heads  are  contained  in  the  "heads  output"  files, 
and  the  node  numbers. 

3.  The  program  "FREQ33.EXE"  is  run  using  the  two  input  files 
described  in  steps  #1  and  #2.  An  output  file  is  specified,  that 
contains  the  cumulative  occurrences  (stage-frequency)  of  storm 
tide  elevations  at  each  of  the  prescribed  nodes. 

4.  The  output  file  created  is  edited,  and  individual  stage- frequency 
files  created  for  each  node. 

5.  The  program  "CUMP2.EXE"  is  run  to  rescale  the  x-axia  in  terms  of 
probabilities . 

6.  The  final  stage- frequency  distributions  are  plotted  using  a 
standard  software  package,  such  as  "GRAPHER"  or  "HARVARD 
GRAPHICS". 
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laps.  IB 

A02806 

A02807 

A05886 

AOOOSO 

A05894 

A05898 

A00195 

A00201 

A00245 

A00813 

A00814 

A00815 

A00865 

A00876 

A00897 

A00950 

A02809 

A02810 

A02813 

A00195 

A00247 


♦Note: 


Table  El 

Magnetic  Tapes  and  Their  Hurricane /Northeaster  Surge  Contents 


Contents 


Elevation  histories  for  Northeasters  1-  2 

Elevation  histories  for  Northeaster  3 

Elevation  histories  for  Northeasters  4-  8 

Elevation  histories  for  Northeasters  9-  11 

Elevation  histories  for  Northeasters  12-  15 

Elevation  histories  for  Northeasters  18-  20 

Elevation  histories  for  Northeaster  21 

Elevation  histories  for  Northeaster  22 

Elevation  histories  for  Northeasters  23-  25 

Elevation  histories  for  Northeasters  26-  28 

Elevation  histories  for  Northeasters  29-  32 

Elevation  histories  for  Northeasters  33 

Elevation  histories  for  Northeasters  34,  36* 

Elevation  histories  for  Northeasters  37-  39 

Elevation  histories  for  Northeasters  40-  42 

Elevation  histories  for  Northeaster  43 

Elevation  histories  for  Hurricanes  1-  84 

Elevation  histories  for  Hurricanes  85-168 

Elevation  histories  for  Hurricanes  169-252 

Statistics  for  Northeasters  1-43  and  * 

Hurricanes  1-218 

Statistics  for  Hurricanes  219-252 


The  simulation  for  Northeaster  35  failed.  There  is  no  output  for 
this  event. 


E5 


Table  E2 

Hurricane  Astronomic  Tides 


Sandv  Hook 

Bridgeport 

l.D. 

Time 

Height _ 

Time 

_ Height _  1 

No. 

Date 

H.M. 

MLW 

NGVD1 

H.M. 

MLW 

NGVD2 

1 

Oct  25,  1978 

2:13 

3.7 

1.91 

5:47 

6.2 

3.27 

8:30 

1.0 

-0.79 

12:04 

1.2 

-1.73 

14:25 

4.1 

2.31 

18:07 

6.3 

3.37 

21:13 

0.6 

-1.19 

24:29 

0.8 

-2.13 

2 

Aug  25,  1978 

0:39 

4.3 

2.51 

4:22 

6.1 

3.17 

6:52 

0.5 

-1.29 

10:34 

0.5 

-2.43 

13:07 

4.7 

2.91 

16:46 

6.7 

3.77 

19:45 

0.8 

-0.99 

23:15 

0.4 

-2.53 

3 

Sept  11,  1978 

1:56 

4.2 

2.41 

5:27 

6.3 

3.37 

8:23 

0.5 

-1.29 

11:41 

0.6 

-2.33 

14:37 

5.1 

3.31 

17:55 

7.0 

4.07 

21:14 

0.3 

-1.49 

24:27 

0.2 

-2.73 

4 

Sept  13,  1978 

4:20 

4.6 

2.81 

1:30 

-0.1 

-3.03 

10:33 

0.0 

-1.79 

7:38 

6.8 

3.87 

16:49 

5.5 

3.71 

13:53 

0.0 

-2.93 

23:10 

-0.4 

-2.19 

20:04 

7.4 

4.47 

S 

July  20,  1978 

1:37 

-1.0 

-2.79 

5:42 

-1.4 

-4.33 

7:37 

5.3 

3.51 

11:51 

7.5 

4.57 

13:44 

-0.7 

-2.49 

18:03 

-0.9 

-3.89 

19:59 

6.4 

4.61 

24:12 

7.9 

4.97 

1  Sandy  Hook  -  Mean  low  water  is  1.79  ft  below  NGVD 

2  Bridgeport  -  Mean  low  water  is  2.93  ft  below  NGVD 
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Table  E3 

Northeaster  Astronomic  Tides 


l.D. 

No.  pate 

1  Jan  17,  1978 

2  Mar  2,  1978 

3  Feb  4,  1978 

4  Dec  4,  1978 

5  Jan  8,  1978 


Sandv  Hook 

BridaeDort 

Heioht 

Time 

Heioht - 

H.M. 

MLW 

JIGVD1 

H-.Mj- 

MLW 

NGVD2 

2:09 

4.1 

2.31 

5:49 

6.4 

3.47 

8:43 

0.3 

-1.49 

12:15 

0.1 

-2.83 

14:31 

3.6 

1.81 

18:19 

5.6 

2.67 

20:52 

0.0 

-1.79 

24:31 

0.2 

-2.73 

0:59 

4.6 

2.81 

4:17 

6.6 

3.67 

7:32 

0.1 

-1.69 

10:49 

-0.1 

-3.03 

13:24 

3.9 

2.11 

16:55 

5.9 

2.97 

19:48 

0.1 

-1.69 

23:09 

0.2 

-2.73 

4:38 

5.0 

3.21 

1:37 

-0.3 

-3.23 

11:02 

-0.7 

-2.49 

7:51 

7.0 

4.07 

17:07 

4.2 

2.41 

14:20 

-0.8 

-3.73 

23:18 

-0.7 

-2.49 

20:26 

6.3 

3.37 

4:28 

-0.5 

-2.29 

1:42 

6.9 

3.97 

10:45 

5.3 

3.51 

7:52 

-0.3 

-3.23 

17:09 

-0.8 

-2.59 

14:03 

7.3 

4.37 

23:19 

4.4 

2.61 

20:29 

-0.8 

-3.73 

0:29 

-1.1 

-2.89 

3:47 

-0.9 

-3.83 

6:46 

5.8 

4.01 

10:00 

7.8 

4.87 

13:12 

-1.3 

-3.09 

16:27 

-1.5 

-4.43 

19:  J? 

4.8 

3.01 

22:33 

6.9 

3.97 

1  Sandy  Hook  -  mean  low  water  is  1.79  ft  below  NGVD 

2  Bridgeport  -  mean  low  water  is  2.93  ft  below  NGVD 
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